title: Neutron Beam Tests of Barium Fluoride Crystal for Dark Matter Direct Detection

abstract: In order to test the capabilities of Barium Fluoride (BaF2) Crystal for dark matter direct detection, nuclear recoils are studied with mono-energetic neutron beam. The energy spectra of nuclear recoils, quenching factors for elastic scattering neutrons and discrimination capability between neutron inelastic scattering events and {\gamma} events are obtained for various recoil energies of the F content in BaF2.


title: Powerful and Tunable THz Emitters Based on the Fe/Pt Magnetic Heterostructure

abstract: In this work, we report our study on the THz emission in Fe/Pt magnetic heterostructures. We have carried out a comprehensive investigation of THz emission from Fe/Pt magnetic heterostructures, employing time-domain THz spectroscopy. We reveal that by properly tuning the thickness of Fe or Pt layer, THz emission can be greatly improved in this type of heterostructure. We further demonstrate that the THz field strength emitted from a newly designed multilayer (Pt/Fe/MgO)$_n$ with n=3 can reach a value of ~1.6 kV/cm, which is comparable to the values from the conventional GaAs antenna with a bias of 4 kV/cm, and the nonlinear crystals, e.g., 100 micrometer GaP and 2 mm ZnTe. For the first time, the intensity and spectrum of THz wave is demonstrated to be tunable by the magnetic field applied on the patterned magnetic Fe/Pt heterostructures. These findings thus promise novel approaches to fabricate powerful and tunable THz emitters based on magnetic heterostructure.


title: Efficient Activity Detection in Untrimmed Video with Max-Subgraph Search

abstract: We propose an efficient approach for activity detection in video that unifies activity categorization with space-time localization. The main idea is to pose activity detection as a maximum-weight connected subgraph problem. Offline, we learn a binary classifier for an activity category using positive video exemplars that are "trimmed" in time to the activity of interest. Then, given a novel \emph{untrimmed} video sequence, we decompose it into a 3D array of space-time nodes, which are weighted based on the extent to which their component features support the learned activity model. To perform detection, we then directly localize instances of the activity by solving for the maximum-weight connected subgraph in the test video's space-time graph. We show that this detection strategy permits an efficient branch-and-cut solution for the best-scoring---and possibly non-cubically shaped---portion of the video for a given activity classifier. The upshot is a fast method that can search a broader space of space-time region candidates than was previously practical, which we find often leads to more accurate detection. We demonstrate the proposed algorithm on four datasets, and we show its speed and accuracy advantages over multiple existing search strategies.


title: Development of four kinds of waveguide power divider for S band

abstract: Four kinds of waveguide power dividers with different structures working at 2.856 GHz are developed. By comparing the simulation performances of these four structures, the power divider with matching rod in the middle of its structure has got the best performance and finally be chosen for fabrication. The prototype has got qualified microwave measurement results as well as the high vacuum performance. It also works stable during the whole progress of microwave commissioning.


title: Binary Codes with Locality for Four Erasures

abstract: In this paper, codes with locality for four erasures are considered. An upper bound on the rate of codes with locality with sequential recovery from four erasures is derived. The rate bound derived here is field independent. An optimal construction for binary codes meeting this rate bound is also provided. The construction is based on regular graphs of girth $6$ and employs the sequential approach of locally recovering from multiple erasures. An extension of this construction that generates codes which can sequentially recover from five erasures is also presented.


title: Equation-free analysis of a dynamically evolving multigraph

abstract: In order to illustrate the adaptation of traditional continuum numerical techniques to the study of complex network systems, we use the equation-free framework to analyze a dynamically evolving multigraph. This approach is based on coupling short intervals of direct dynamic network simulation with appropriately-defined lifting and restriction operators, mapping the detailed network description to suitable macroscopic (coarse-grained) variables and back. This enables the acceleration of direct simulations through Coarse Projective Integration (CPI), as well as the identification of coarse stationary states via a Newton-GMRES method. We also demonstrate the use of data-mining, both linear (principal component analysis, PCA) and nonlinear (diffusion maps, DMAPS) to determine good macroscopic variables (observables) through which one can coarse-grain the model. These results suggest methods for decreasing simulation times of dynamic real-world systems such as epidemiological network models. Additionally, the data-mining techniques could be applied to a diverse class of problems to search for a succint, low-dimensional description of the system in a small number of variables.


title: A finite element approach to self-consistent field theory calculations of multiblock polymers

abstract: Self-consistent field theory (SCFT) has proven to be a powerful tool for modeling equilibrium microstructures of soft materials, particularly for multiblock polymers. A very successful approach to numerically solving the SCFT set of equations is based on using a spectral approach. While widely successful, this approach has limitations especially in the context of current technologically relevant applications. These limitations include non-trivial approaches for modeling complex geometries, difficulties in extending to non-periodic domains, as well as non-trivial extensions for spatial adaptivity. As a viable alternative to spectral schemes, we develop a finite element formulation of the SCFT paradigm for calculating equilibrium polymer morphologies. We discuss the formulation and address implementation challenges that ensure accuracy and efficiency. We explore higher order chain contour steppers that are efficiently implemented with Richardson Extrapolation. This approach is highly scalable and suitable for systems with arbitrary shapes. We show spatial and temporal convergence and illustrate scaling on up to 2048 cores. Finally, we illustrate confinement effects for selected complex geometries. This has implications for materials design for nanoscale applications where dimensions are such that equilibrium morphologies dramatically differ from the bulk phases.


title: Quantum thermal transport in stanene

abstract: By way of the non-equilibrium Green's function simulations and analytical expressions, the quantum thermal conductance of stanene is studied. We find that, due to the existence of Dirac fermion in stanene, the ratio of electron thermal conductance and electric conductance becomes a chemical-potential-dependent quantity, violating the Wiedemann-Franz law. This finding is applicable to any two-dimensional (2D) materials that possess massless Dirac fermions. In strong contrast to the negligible electronic contribution in graphene, surprisingly, the electrons and phonons in stanene carry a comparable heat current. The unusual behaviours in stanene widen our knowledge of quantum thermal transport in 2D materials.


title: Discontinuity in the asymptotic behavior of planar orthogonal polynomials under a perturbation of the Gaussian weight

abstract: We consider the orthogonal polynomials, $\{P_n(z)\}_{n=0,1,\cdots}$, with respect to the measure $$|z-a|^{2c} e^{-N|z|^2}dA(z)$$ supported over the whole complex plane, where $a>0$, $N>0$ and $c>-1$. We look at the scaling limit where $n$ and $N$ tend to infinity while keeping their ratio, $n/N$, fixed. The support of the limiting zero distribution is given in terms of certain "limiting potential-theoretic skeleton" of the unit disk. We show that, as we vary $c$, both the skeleton and the asymptotic distribution of the zeros behave discontinuously at $c=0$. The smooth interpolation of the discontinuity is obtained by the further scaling of $c=e^{-\eta N}$ in terms of the parameter $\eta\in[0,\infty).$


title: Capacity Bounds for Networks with Correlated Sources and Characterisation of Distributions by Entropies

abstract: Characterising the capacity region for a network can be extremely difficult. Even with independent sources, determining the capacity region can be as hard as the open problem of characterising all information inequalities. The majority of computable outer bounds in the literature are relaxations of the Linear Programming bound which involves entropy functions of random variables related to the sources and link messages. When sources are not independent, the problem is even more complicated. Extension of Linear Programming bounds to networks with correlated sources is largely open. Source dependence is usually specified via a joint probability distribution, and one of the main challenges in extending linear program bounds is the difficulty (or impossibility) of characterising arbitrary dependencies via entropy functions. This paper tackles the problem by answering the question of how well entropy functions can characterise correlation among sources. We show that by using carefully chosen auxiliary random variables, the characterisation can be fairly "accurate" Using such auxiliary random variables we also give implicit and explicit outer bounds on the capacity of networks with correlated sources. The characterisation of correlation or joint distribution via Shannon entropy functions is also applicable to other information measures such as Renyi entropy and Tsallis entropy.


title: $a_0(980)$ photoproduction in the coupled channel model

abstract: We calculated the coupled channel photoproduction amplitudes of the scalar isovector resononace $a_0(980)$, which include the $\pi\eta$, $K^+K^-$ and $K^0\overline{K^0}$ intermediate states. Based on them we obtained the mass distribution in the $\pi\eta$ channel at photon energy $E_\gamma$=7 GeV, corresponding to energy accissible in new JLab experiments CLAS12 and GlueX. We also analyzed the shape of mass distribution as a function of the phase of the background amplitude.


title: Reaching Agreement in Quantum Hybrid Networks

abstract: We consider a basic quantum hybrid network model consisting of a number of nodes each holding a qubit, for which the aim is to drive the network to a consensus in the sense that all qubits reach a common state. Projective measurements are applied serving as control means, and the measurement results are exchanged among the nodes via classical communication channels. We show how to carry out centralized optimal path planning for this network with all-to-all classical communications, in which case the problem becomes a stochastic optimal control problem with a continuous action space. To overcome the computation and communication obstacles facing the centralized solutions, we also develop a distributed Pairwise Qubit Projection (PQP) algorithm, where pairs of nodes meet at a given time and respectively perform measurements at their geometric average. We show that the qubit states are driven to a consensus almost surely along the proposed PQP algorithm, and that the expected qubit density operators converge to the average of the network's initial values.


title: Higher cohomology operations and R-completion

abstract: Let $R=\mathbb{F}_p$ or a field of characteristic $0$. For each $R$-good topological space $Y$, we define a collection of higher cohomology operations which, together with the cohomology algebra $H^*(Y;R)$ suffice to determine $Y$ up to $R$-completion. We also provide a similar collection of higher cohomology operations which determine when two maps $f_0,f_1: Z\to Y$ between $R$-good spaces(inducing the same algebraic homomorphism $H^*(Y;R)\to H^*(Z;R)$) are $R$-equivalent.


title: Comparison of the Smith-Purcell Radiation Yield for Different Models

abstract: Smith-Purcell radiation is used in several applications including the measurement of the longitudinal profile of electron bunches. A correct reconstruction of such profile requires a good understanding of the underlying model. We have compared the leading models of Smith-Purcell radiation and shown that they are in agreement within the experimental errors.


title: Decay of the Higgs boson to $\tau^- \tau^+$ and non-Hermiticity of the Yukawa interaction

abstract: The issue of Hermiticity of the Higgs boson interaction with fermions is addressed. A model for non-Hermitian Yukawa interaction is proposed and approximation of one fermion generation is considered. Symmetry properties of the corresponding $h f \bar{f}$ Lagrangian with respect to the discrete P, C and T transformations are analyzed, and the modified Dirac equation for the free fermion is studied. Longitudinal polarization of the fermions in the decay $h \to f \bar{f}$, which arises due to non-Hermiticity of the $h f \bar{f}$ interaction, is discussed. It is suggested to study effects of this non-Hermiticity in the decay $h \to \tau^- \tau^+ \to \mu^- {\bar \nu}_\mu \nu_\tau \, \mu^+ \nu_\mu {\bar \nu}_\tau$, for which observables (asymmetries) are constructed which take nonzero values for a non-Hermitian $h \tau^- \tau^+$ interaction. These asymmetries are analyzed for various configurations of the muon energies.


title: Elastic Model for Dinucleosome Structure and Energy

abstract: The equilibrium structure of a Dinucleosome is studied using an elastic model that takes into account the force and torque balance conditions. Using the proper boundary conditions, it is found that the conformational energy of the problem does not depend on the length of the linker DNA. In addition it is shown that the two histone octamers are almost perpendicular to each other and the linker DNA in short lengths is almost straight. These findings could shed some light on the role of DNA elasticity in the chromatin structure.


title: Hypergraph Modelling for Geometric Model Fitting

abstract: In this paper, we propose a novel hypergraph based method (called HF) to fit and segment multi-structural data. The proposed HF formulates the geometric model fitting problem as a hypergraph partition problem based on a novel hypergraph model. In the hypergraph model, vertices represent data points and hyperedges denote model hypotheses. The hypergraph, with large and "data-determined" degrees of hyperedges, can express the complex relationships between model hypotheses and data points. In addition, we develop a robust hypergraph partition algorithm to detect sub-hypergraphs for model fitting. HF can effectively and efficiently estimate the number of, and the parameters of, model instances in multi-structural data heavily corrupted with outliers simultaneously. Experimental results show the advantages of the proposed method over previous methods on both synthetic data and real images.


title: Classical limit of diagonal form factors and HHL correlators

abstract: We propose an expression for the classical limit of diagonal form factors in which we integrate the corresponding observable over the moduli space of classical solutions. In infinite volume the integral has to be regularized by proper subtractions and we present the one, which corresponds to the classical limit of the connected diagonal form factors. In finite volume the integral is finite and can be expressed in terms of the classical infinite volume diagonal form factors and subvolumes of the moduli space. We analyze carefully the periodicity properties of the finite volume moduli space and found a classical analogue of the Bethe-Yang equations. By applying the results to the heavy-heavy-light three point functions we can express their strong coupling limit in terms of the classical limit of the sine-Gordon diagonal form factors.


title: Flow of Power-Law Liquids in a Hele-Shaw Cell Driven by Non-Uniform Electroosmotic Slip in the Case of Strong Depletion

abstract: We analyze flow of non-Newtonian fluids in a Hele-Shaw cell, subjected to spatially non-uniform electroosmotic slip. Motivated by their potential use for increasing the characteristic pressure fields, we specifically focus on power-law fluids with wall depletion properties. We derive a p-Poisson equation governing the pressure field, as well as a set of linearized equations representing its asymptotic approximation for weakly non-Newtonian behavior. To investigate the effect of non-Newtonian properties on the resulting fluidic pressure and velocity, we consider several configurations in one- and two-dimensions, and calculate both exact and approximate solutions. We show that the asymptotic approximation is in good agreement with exact solutions even for fluids with significant non-Newtonian behavior, allowing its use in the analysis and design of microfluidic systems involving electro-kinetic transport of such fluids.


title: Monopole excitation of the Hoyle state and linear-chain state in 12C

abstract: Background: Two new $J^\pi=0^+$ states are recently observed around a few MeV above the Hoyle state (the second $0^+$ state in 12C). The characters of them are only poorly discussed in theory and are still mysterious. Purpose: I give for the first time a comprehensive understanding of the structures of the $0^+$ states by analyzing their wave functions, and discuss relationship with the Hoyle state, similarities, and differences between the states. Method: I extend a microscopic $\alpha$-cluster model called Tohsaki-Horiuchi-Schuck-R\"opke (THSR) wave function so as to incorporate $2\alpha+\alpha$ asymmetric configuration explicitly. The so-called $r^2$-constraint method to effectively eliminate spurious continuum components is also used. Results: The $0_3^+$ state is shown to have a very large squared overlap with a single configuration of the extended THSR wave function in an orthogonal space to the Hoyle state as well as to the ground state. The $0_4^+$ state has a maximal squared overlap with a single extended THSR wave function with an extremely prolately-deformed shape. Conclusions: The $0_3^+$ state appears as a family of the Hoyle state to have a higher nodal structure in the internal motions of the $3\alpha$ clusters, due to the orthogonalization to the Hoyle state. The $0_4^+$ state dominantly has a linear-chain structure, where the $3\alpha$ clusters move freely in a non-localized way, like a one-dimensional gas of the $3\alpha$ clusters.


title: Barycentric Subspace Analysis on Manifolds

abstract: This paper investigates the generalization of Principal Component Analysis (PCA) to Riemannian manifolds. We first propose a new and general type of family of subspaces in manifolds that we call barycentric subspaces. They are implicitly defined as the locus of points which are weighted means of $k+1$ reference points. As this definition relies on points and not on tangent vectors, it can also be extended to geodesic spaces which are not Riemannian. For instance, in stratified spaces, it naturally allows principal subspaces that span several strata, which is impossible in previous generalizations of PCA. We show that barycentric subspaces locally define a submanifold of dimension k which generalizes geodesic subspaces.Second, we rephrase PCA in Euclidean spaces as an optimization on flags of linear subspaces (a hierarchy of properly embedded linear subspaces of increasing dimension). We show that the Euclidean PCA minimizes the Accumulated Unexplained Variances by all the subspaces of the flag (AUV). Barycentric subspaces are naturally nested, allowing the construction of hierarchically nested subspaces. Optimizing the AUV criterion to optimally approximate data points with flags of affine spans in Riemannian manifolds lead to a particularly appealing generalization of PCA on manifolds called Barycentric Subspaces Analysis (BSA).


title: Tight Lower Bounds for Multiplicative Weights Algorithmic Families

abstract: We study the fundamental problem of prediction with expert advice and develop regret lower bounds for a large family of algorithms for this problem. We develop simple adversarial primitives, that lend themselves to various combinations leading to sharp lower bounds for many algorithmic families. We use these primitives to show that the classic Multiplicative Weights Algorithm (MWA) has a regret of $\sqrt{\frac{T \ln k}{2}}$, there by completely closing the gap between upper and lower bounds. We further show a regret lower bound of $\frac{2}{3}\sqrt{\frac{T\ln k}{2}}$ for a much more general family of algorithms than MWA, where the learning rate can be arbitrarily varied over time, or even picked from arbitrary distributions over time. We also use our primitives to construct adversaries in the geometric horizon setting for MWA to precisely characterize the regret at $\frac{0.391}{\sqrt{\delta}}$ for the case of $2$ experts and a lower bound of $\frac{1}{2}\sqrt{\frac{\ln k}{2\delta}}$ for the case of arbitrary number of experts $k$.


title: Form Follows Function -- Do algorithms and applications challenge or drag behind the hardware evolution?

abstract: We summarise some of the key statements made at the workshop Form Follows Function at ISC High Performance 2016. The summary highlights what type of co-design the presented projects experience; often in the absence of an explicit co-design agenda. Their software development picks up hardware trends but it also influences the hardware development. Observations illustrate that this cycle not always is optimal for both sides as it is not proactively steered. Key statements characterise ideas how it might be possible to integrate both hardware and software creation closer to the best of both worlds---again even without classic co-design in mind where new pieces of hardware are created. The workshop finally identified three development idioms that might help to improve software and system design with respect to emerging hardware.


title: Maximum principles for nonlocal parabolic Waldenfels operators

abstract: As a class of L\'evy type Markov generators, nonlocal Waldenfels operators appear naturally in the context of investigating stochastic dynamics under L\'evy fluctuations and constructing Markov processes with boundary conditions (in particular the construction with jumps). This work is devoted to prove the weak and strong maximum principles for `parabolic' equations with nonlocal Waldenfels operators. Applications in stochastic differential equations with $\alpha$-stable L\'evy processes are presented to illustrate the maximum principles.


title: Non-Markovian dynamics in the extended cluster spin-1/2 XX chain

abstract: We study the dynamics of entanglement in the extended cluster spin-1/2 XX chain, equivalent to a 1D spin-1/2 XX model with three-spin interaction (TSI). Selecting the nearest neighbor pair spins as an open quantum system, the rest of the chain plays the role of environment. The two-spin Heisenberg and the TSI interaction are responsible for coupling between system and environment. We show the existence of a critical value in the TSI, where the dynamics of concurrence changes from Markovian to the non-Markovian. In the region with the non-Markovian dynamics, entanglement sudden death in the system is observed. By focusing on the nearest neighbor pair spins of the environment, we have showed that the dynamics of entanglement in the environment is sensitive to the Markovian and non-Markovian regions.


title: Suppression of scattering for small dielectric particles: an anapole mode and invisibility

abstract: We reveal that an isotropic homogeneous subwavelength particle with a high refractive index can produce ultra-weak total scattering due to vanishing contribution of the electric dipole moment. This effect can be explained with the help of the Fano resonance and scattering efficiency associated with the excitation of an anapole mode. The latter is a nonradiative mode emerging from destructive interference of electric and toroidal dipole moments, and it can be employed for a design of highly transparent optical materials.


title: Boosting Magnetic Reconnection by Viscosity and Thermal Conduction

abstract: Nonlinear evolution of magnetic reconnection is investigated by means of magnetohydrodynamic simulations including uniform resistivity, uniform viscosity, and anisotropic thermal conduction. When viscosity exceeds resistivity (the magnetic Prandtl number Prm > 1), the viscous dissipation dominates outflow dynamics and leads to the decrease in the plasma density inside a current sheet. The low-density current sheet supports the excitation of the vortex. The thickness of the vortex is broader than that of the current for Prm > 1. The broader vortex flow more efficiently carries the upstream magnetic flux toward the reconnection region, and consequently boosts the reconnection. The reconnection rate increases with viscosity provided that thermal conduction is fast enough to take away the thermal energy increased by the viscous dissipation (the fluid Prandtl number Pr < 1). The result suggests the need to control the Prandtl numbers for the reconnection against the conventional resistive model.


title: A microscopic theory for ultra-near-field radiation

abstract: Using the nonequilibrium Green's function (NEGF) formalism, we propose a microscopic theory for near-field radiative heat transfer between metal plates. Tight-binding models for the electrons are coupled to the electromagnetic field continuum. Our approach differs from the established ones based on fluctuational electrodynamics, in that it describes truly nonequilibrium steady states, and is nonlocal in system's dielectric properties. For a two quantum-dot model a new length scale emerges at which the heat current shows a peak. This length scale is related to the physics of parallel plate capacitors. The three-dimensional model results are consistent with the theory of Polder and van Hove except at very short distances.


title: Bipartisanship breakdown, functional networks and forensic analysis in Spanish 2015 and 2016 national elections

abstract: In this paper we present a social network and forensic analysis of the vote counts of Spanish national elections that took place in December 2015 and their sequel in June 2016. Vote counts are extracted at the level of municipalities, yielding an unusually high resolution dataset with over 8000 samples. We initially consider the phenomenon of Bipartisanship breakdown by analysing spatial distributions of several Bipartisanship indices. We find that such breakdown is more prominent close to cosmopolite and largely populated areas and less important in rural areas where Bipartisanship still prevails, and its evolution mildly consolidates in the 2016 round, with some evidence of Bipartisanship reinforcement which we hypothesize to be due to psychological mechanisms of risk aversion. On a third step we explore to which extent vote data are faithful by applying forensic techniques to vote statistics. We first explore the conformance of first digit distributions to Benford's law for each of the main political parties. The results and interpretations are mixed and vary across different levels of aggregation, finding a general good quantitative agreement at the national scale for both municipalities and precincts but finding systematic nonconformance at the level of individual precincts. As a complementary metric, we further explore the co-occurring statistics of voteshare and turnout, finding a mild tendency in the clusters of the conservative party to smear out towards the area of high turnout and voteshare, what has been previously interpreted as a possible sign of incremental fraud. In every case results are qualitatively similar between 2015 and 2016 elections.


title: Staggered Time Average Algorithm for Stochastic Non-smooth Optimization with O(1/T) Convergence

abstract: Stochastic non-smooth convex optimization constitutes a class of problems in machine learning and operations research. This paper considers minimization of a non-smooth function based on stochastic subgradients. When the function has a locally polyhedral structure, a staggered time average algorithm is proven to have O(1/T) convergence rate. A more general convergence result is proven when the locally polyhedral assumption is removed. In that case, the convergence bound depends on the curvature of the function near the minimum. Finally, the locally polyhedral assumption is shown to improve convergence beyond O(1/T) for a special case of deterministic problems.


title: Connected formulas for amplitudes in standard model

abstract: Witten's twistor string theory has led to new representations of S-matrix in massless QFT as a single object, including Cachazo-He-Yuan formulas in general and connected formulas in four dimensions. As a first step towards more realistic processes of the standard model, we extend the construction to QCD tree amplitudes with massless quarks and those with a Higgs boson. For both cases, we find connected formulas in four dimensions for all multiplicities which are very similar to the one for Yang-Mills amplitudes. The formula for quark-gluon color-ordered amplitudes differs from the pure-gluon case only by a Jacobian factor that depends on flavors and orderings of the quarks. In the formula for Higgs plus multi-parton amplitudes, the massive Higgs boson is effectively described by two additional massless legs which do not appear in the Parke-Taylor factor. The latter also represents the first twistor-string/connected formula for form factors.


title: Hurwitz spaces and liftings to the Valentiner group

abstract: We study the components of the Hurwitz scheme of ramified coverings of $\mathbb{P}^1$ with monodromy given by the alternating group $A_6$ and elements in the conjugacy class of product of two disjoint cycles. In order to detect the connected components of the Hurwitz scheme, inspired by the case of the spin structures studied by Fried for the $3$-cycles, we use as invariant the lifting to the Valentiner group, triple covering of $A_6$. We prove that the Hurwitz scheme has two irreducible components when the genus of the covering is greater than zero, in accordance with the asymptotic solution found by Bogomolov and Kulikov.


title: Quantifying some properties of the stochastic quantum force

abstract: We consider for clarity the simple case of the one dimensional non-relativistic Schr\"{o}dinger equation and regard it as an ensemble mean representation of the stochastic motion of a single particle in a vacuum, subject to an undefined stochastic quantum force. By analyzing the Bohm potential it is found that the imaginary part of the quantum momentum is the root mean square fluctuation of the particle around its mean velocity, where the latter is the real part of the quantum momentum. The local mean of the quantum force is found to be proportional to the third spatial derivative of the probability density function, and its associated pressure to the second spatial derivative. The latter is decomposed from the single particle diluted gas pressure, and this pressure partition allows the interpretation of the quantum Bohm potential as the energy required to put a particle in a bath of fluctuating vacuum at constant entropy and volume. The stochastic force expectation value is zero and is uncorrelated with the particle location, thus does not perform work on average. Nonetheless it is anti-correlated with volume and this anti-correlation leads to a new type of an Heisenberg like relation. We imply the dynamic Gaussian solution to the Schr\"{o}dinger equation as a simple example for exploring the mean properties of this quantum force. Still an interesting remained open task is the identification of the stochastic law that leads to these obtained mean properties.


title: Model reduction of linear time-varying systems with applications for moving loads

abstract: In this paper we consider different model reduction techniques for systems with moving loads. Due to the time-dependency of the input and output matrices, the application of time-varying projection matrices for the reduction offers new degrees of freedom, which also come along with some challenges. This paper deals with both simple methods for the reduction of particular linear time-varying systems, as well as with a more advanced technique considering the emerging time derivatives.


title: Cartan subalgebras of tensor products of free quantum group factors with arbitrary factors

abstract: Let $\mathbb{G}$ be a free (unitary or orthogonal) quantum group. We prove that for any non-amenable subfactor $N\subset L^\infty(\mathbb{G})$, which is an image of a faithful normal conditional expectation, and for any $\sigma$-finite factor $B$, the tensor product $N \otimes B$ has no Cartan subalgebras. This generalizes our previous work that provides the same result when $B$ is finite. In the proof, we establish Ozawa--Popa and Popa--Vaes's weakly compact action on the continuous core of $N \otimes B$ as the one relative to B, by using an operator valued weight to B and the central weak amenability of $\mathbb{G}$.


title: Remarks on endpoint Strichartz estimates for Schr\"odinger equations with the critical inverse-square potential

abstract: The purpose of this paper is to study the validity of global-in-time Strichartz estimates for the Schr\"odinger equation on $\mathbb{R}^n$, $n\ge3$, with the negative inverse-square potential $-\sigma|x|^{-2}$ in the critical case $\sigma=(n-2)^2/4$. It turns out that the situation is different from the subcritical case $\sigma<(n-2)^2/4$ in which the full range of Strichartz estimates is known to be hold. More precisely, splitting the solution into the radial and non-radial parts, we show that (i) the radial part satisfies a weak-type endpoint estimate, which can be regarded as an extension to higher dimensions of the endpoint Strichartz estimate with radial data for the two-dimensional free Schr\"odinger equation; (ii) other endpoint estimates in Lorentz spaces for the radial part fail in general; (iii) the non-radial part satisfies the full range of Strichartz estimates.


title: Affine embeddings of Cantor sets on the line

abstract: Let $s\in (0,1)$, and let $F\subset \mathbb{R}$ be a self similar set such that $0 < \dim_H F \leq s$ . We prove that there exists $\delta= \delta(s) >0$ such that if $F$ admits an affine embedding into a homogeneous self similar set $E$ and $0 \leq \dim_H E - \dim_H F < \delta$ then (under some mild conditions on $E$ and $F$) the contraction ratios of $E$ and $F$ are logarithmically commensurable. This provides more evidence for a Conjecture of Feng, Huang, and Rao, that states that these contraction ratios are logarithmically commensurable whenever $F$ admits an affine embedding into $E$ (under some mild conditions). Our method is a combination of an argument based on the approach of Feng, Huang, and Rao, with a new result by Hochman, which is related to the increase of entropy of measures under convolutions.


title: Josephson Quantum Heat Engine

abstract: The design of a mesoscopic self-oscillating heat engine that works thanks to purely quantum effects is presented. The proposed scheme is amenable to experimental implementation with current state-of-the-art nanotechnology and materials. One of the main features of the structure is its versatility: The engine can deliver work to a generic load without galvanic contact. This makes it a promising building block for low-temperature on-chip energy management applications. The heat engine consists of a circuit featuring a thermoelectric element based on a ferromagnetic insulator-superconductor tunnel junction and a Josephson weak link that realizes a purely quantum DC/AC converter. This enables contactless transfer of work to the load (a generic RL circuit). The performance of the heat engine is investigated as a function of the thermal gradient applied to the thermoelectric junction. Power up to $1$ pW can be delivered to a load $R_L=10\ \Omega$.

categories: cond-mat.mes-hall cond-mat.stat-mech cond-mat.supr-con quant-ph


title: Microscopic theory of electron absorption by plasma-facing surfaces

abstract: We describe a method for calculating the probability with which the wall of a plasma absorbs an electron at low energy. The method, based on an invariant embedding principle, expresses the electron absorption probability as the probability for transmission through the wall's long-range surface potential times the probability to stay inside the wall despite of internal backscattering. To illustrate the approach we apply it to a \SiOTwo\ surface. Besides emission of optical phonons inside the wall we take elastic scattering at imperfections of the plasma-wall interface into account and obtain absorption probabilities significantly less than unity in accordance with available electron-beam scattering data but in disagreement with the widely used perfect absorber model.


title: Exact Casimir interaction of perfectly conducting three-spheres in four euclidean dimensions

abstract: Exploiting conformal symmetry, we derive a simple exact formula for the classical electromagnetic Casimir interaction of two perfectly conducting three-spheres, including the sphere-plate geometry as a special case, in four euclidean dimensions. We verify that the short distance expansion of the Casimir energy agrees to leading order with the Proximity Force Approximation (PFA), while the next-to-leading-order is in agreement with a recently proposed derivative expansion of the Casimir energy. At the next-to-next-to-leading order we find a non-analytic correction to PFA, which for a sphere-plate system is of the order of $(d/R)^{3/2} \log(d/R)$, where $d$ is the separation and $R$ the sphere radius.


title: Orbital-dependent charge dynamics in MnP revealed by optical study

abstract: Unconventional superconductivity often emerges at the border of long-range magnetic orders. Understanding the low-energy charge dynamics may provide crucial information on the formation of superconductivity. Here we report the unpolarized/polarized optical conductivity study of high quality MnP single crystals at ambient pressure. Our data reveal two types of charge carriers with very different lifetimes. In combination with the first-principles calculations, we show that the short-lifetime carriers have flat Fermi sheets which become gapped in the helimagnetic phase, causing a dramatic change in the low-frequency optical spectra, while the long-lifetime carriers are anisotropic three-dimensional like which are little affected by the magnetic transitions and provide major contributions to the transport properties. This orbital-dependent charge dynamics originates from the special crystal structure of MnP and may have an influence on the unconventional superconductivity and its interplay with helimagnetism at high pressures.


title: Spin- and angle-resolved photoemission on the topological Kondo Insulator candidate: SmB6

abstract: Topological Kondo insulators are a new class of topological insulators in which metallic surface states protected by topological invariants reside in the bulk band gap at low temperatures. Unlike other three-dimensional topological insulators, a truly insulating bulk state, which is critical for potential applications in next-generation electronic devices, is guaranteed by many-body effects in the topological Kondo insulator. Furthermore, the system has strong electron correlations that can serve as a testbed for interacting topological theories. This topical review focuses on recent advances in the study of SmB6, the most promising candidate for a topological Kondo insulator, from the perspective of spin- and angle-resolved photoemission spectroscopy with highlights of some important transport results.


title: Equidistribution of zeros of random polynomials

abstract: We study the asymptotic distribution of zeros for the random polynomials $P_n(z) = \sum_{k=0}^n A_k B_k(z)$, where $\{A_k\}_{k=0}^{\infty}$ are non-trivial i.i.d. complex random variables. Polynomials $\{B_k\}_{k=0}^{\infty}$ are deterministic, and are selected from a standard basis such as Szeg\H{o}, Bergman, or Faber polynomials associated with a Jordan domain $G$ bounded by an analytic curve. We show that the zero counting measures of $P_n$ converge almost surely to the equilibrium measure on the boundary of $G$ if and only if $\mathbb{E}[\log^+|A_0|]<\infty$.


title: Effect of local environment on moment formation in iron silicides

abstract: The effect of local environment on the formation of magnetic moments on Fe atoms in iron silicides are studied by combination of ab initio and model calculations. The suggested model includes all Fe d- and Si p-orbitals, intra-atomic Coulomb interactions, inter-atomic Fe-Fe exchange and hopping of electrons to nearest and next nearest neighboring atoms and takes into account all symmetries within the Slater-Koster scheme. The parameters of the model are found from the requirement that self-consistent moments on atoms and density of states found from ab initio and model calculations within the Hartree-Fock approximation are close to each other as much as possible. Contrary to the commonly accepted statement that an increase of the Si concentration within nearest environment of Fe atoms in the ordered Fe3Si and FexSi1-x alloys leads to a decrease of Fe magnetic moment we find that a crucial role in the formation of magnetic moments is played by second coordination sphere of Fe atoms. Particularly, the Fe atoms have higher magnetic moments in amorphous films compared to the epitaxial ones due to decrease in the number of iron-atoms in the next nearest environment. Both our model and ab initio calculations confirm existence of known spin crossover with pressure and predict second crossover at higher pressure


title: Classifying Variable-Length Audio Files with All-Convolutional Networks and Masked Global Pooling

abstract: We trained a deep all-convolutional neural network with masked global pooling to perform single-label classification for acoustic scene classification and multi-label classification for domestic audio tagging in the DCASE-2016 contest. Our network achieved an average accuracy of 84.5% on the four-fold cross-validation for acoustic scene recognition, compared to the provided baseline of 72.5%, and an average equal error rate of 0.17 for domestic audio tagging, compared to the baseline of 0.21. The network therefore improves the baselines by a relative amount of 17% and 19%, respectively. The network only consists of convolutional layers to extract features from the short-time Fourier transform and one global pooling layer to combine those features. It particularly possesses neither fully-connected layers, besides the fully-connected output layer, nor dropout layers.


title: Incremental Factorization Machines for Persistently Cold-starting Online Item Recommendation

abstract: Real-world item recommenders commonly suffer from a persistent cold-start problem which is caused by dynamically changing users and items. In order to overcome the problem, several context-aware recommendation techniques have been recently proposed. In terms of both feasibility and performance, factorization machine (FM) is one of the most promising methods as generalization of the conventional matrix factorization techniques. However, since online algorithms are suitable for dynamic data, the static FMs are still inadequate. Thus, this paper proposes incremental FMs (iFMs), a general online factorization framework, and specially extends iFMs into an online item recommender. The proposed framework can be a promising baseline for further development of the production recommender systems. Evaluation is done empirically both on synthetic and real-world unstable datasets.


title: Cubic scaling $GW$: towards fast quasiparticle calculations

abstract: Within the framework of the full potential projector-augmented wave methodology, we present a promising low-scaling $GW$ implementation. It allows for quasiparticle calculations with a scaling that is cubic in the system size and linear in the number of $k$ points used to sample the Brillouin zone. This is achieved by calculating the polarizability and self-energy in the real space and imaginary time domain. The transformation from the imaginary time to the frequency domain is done by an efficient discrete Fourier transformation with only a few nonuniform grid points. Fast Fourier transformations are used to go from real space to reciprocal space and vice versa. The analytic continuation from the imaginary to the real frequency axis is performed by exploiting Thiele's reciprocal difference approach. Finally, the method is applied successfully to predict the quasiparticle energies and spectral functions of typical semiconductors (Si, GaAs, SiC, and ZnO), insulators (C, BN, MgO, and LiF), and metals (Cu and SrVO$_3$). The results are compared with conventional $GW$ calculations. Good agreement is achieved, highlighting the strength of the present method.


title: A simple laser locking system based on a field-programmable gate array

abstract: Frequency stabilization of laser light is crucial in both scientific and industrial applications. Technological developments now allow analog laser stabilization systems to be replaced with digital electronics such as field-programmable gate arrays, which have recently been utilized to develop such locking systems. We have developed a frequency stabilization system based on a field-programmable gate array, with emphasis on hardware simplicity, which offers a user-friendly alternative to commercial and previous home-built solutions. Frequency modulation, lock-in detection and a proportional-integral-derivative controller are programmed on the field-programmable gate array and only minimal additional components are required to frequency stabilize a laser. The locking system is administered from a host-computer which provides comprehensive, long-distance control through a versatile interface. Various measurements were performed to characterize the system. The linewidth of the locked laser was measured to be $0.7 \pm 0.1 \; \rm{MHz}$ with a settling time of $10 \; \rm{ms}$. The system can thus fully match laser systems currently in use for atom trapping and cooling applications.


title: A two-sided academic landscape: portrait of highly-cited documents in Google Scholar (1950-2013)

abstract: The main objective of this paper is to identify the set of highly-cited documents in Google Scholar and to define their core characteristics (document types, language, free availability, source providers, and number of versions), under the hypothesis that the wide coverage of this search engine may provide a different portrait about this document set respect to that offered by the traditional bibliographic databases. To do this, a query per year was carried out from 1950 to 2013 identifying the top 1,000 documents retrieved from Google Scholar and obtaining a final sample of 64,000 documents, of which 40% provided a free full-text link. The results obtained show that the average highly-cited document is a journal article or a book (62% of the top 1% most cited documents of the sample), written in English (92.5% of all documents) and available online in PDF format (86.0% of all documents). Yet, the existence of errors especially when detecting duplicates and linking cites properly must be pointed out. The fact of managing with highly cited papers, however, minimizes the effects of these limitations. Given the high presence of books, and to a lesser extend of other document types (such as proceedings or reports), the research concludes that Google Scholar data offer an original and different vision of the most influential academic documents (measured from the perspective of their citation count), a set composed not only by strictly scientific material (journal articles) but academic in its broad sense


title: A bifurcation analysis for the Lugiato-Lefever equation

abstract: The Lugiato-Lefever equation is a cubic nonlinear Schr\"odinger equation, including damping, detuning and driving, which arises as a model in nonlinear optics. We study the existence of stationary waves which are found as solutions of a four-dimensional reversible dynamical system in which the evolutionary variable is the space variable. Relying upon tools from bifurcation theory and normal forms theory, we discuss the codimension 1 bifurcations. We prove the existence of various types of steady solutions, including spatially localized, periodic, or quasi-periodic solutions.


title: The denominators of harmonic numbers

abstract: The denominators $d_n$ of the harmonic number $1+\frac12+\frac13+\cdots+\frac1n$ do not increase monotonically with~$n$. It is conjectured that $d_n=D_n={\rm LCM}(1,2,\ldots,n)$ infinitely often. For an odd prime $p$, the set $\{n:pd_n|D_n\}$ has a harmonic density and, for $2<p_1<p_2<\cdots<p_k$, there exists $n$ such that $p_1p_2\cdots p_kd_n|D_n$.


title: On the Structure of Equilibrium Strategies in Dynamic Gaussian Signaling Games

abstract: This paper analyzes a finite horizon dynamic signaling game motivated by the well-known strategic information transmission problems in economics. The mathematical model involves information transmission between two agents, a sender who observes two Gaussian processes, state and bias, and a receiver who takes an action based on the received message from the sender. The players incur quadratic instantaneous costs as functions of the state, bias and action variables. Our particular focus is on the Stackelberg equilibrium, which corresponds to information disclosure and Bayesian persuasion problems in economics. Prior work solved the static game, and showed that the Stackelberg equilibrium is achieved by pure strategies that are linear functions of the state and the bias variables. The main focus of this work is on the dynamic (multi-stage) setting, where we show that the existence of a pure strategy Stackelberg equilibrium, within the set of linear strategies, depends on the problem parameters. Surprisingly, for most problem parameters, a pure linear strategy does not achieve the Stackelberg equilibrium which implies the existence of a trade-off between exploiting and revealing information, which was also encountered in several other asymmetric information games.


title: Significance of Natural Scene Statistics in Understanding the Anisotropies of Perceptual Filling-in at the Blind Spot

abstract: Psychophysical experiments reveal our horizontal preference in perceptual filling-in at the blind spot. On the other hand, vertical preference is exhibited in the case of tolerance in filling-in. What causes this anisotropy in our perception? Building upon the general notion, that the functional properties of the early visual system are shaped by the innate specification as well as the statistics of the environment, we reasoned that the anisotropy in filling-in could be understood in terms of anisotropy in orientation distribution inherent to natural scene statistics. We examined this proposition by investigating filling-in of bar stimuli on a Hierarchical Predictive Coding model network. In response to bar stimuli, the model network, trained with natural images, exhibited anisotropic filling-in performance at the blind spot similar to reported in psychophysical experiments i.e. horizontal preference in filling-in and vertical preference in tolerance of filling-in. We suggest that the over-representation of horizontal contours in the natural scene contribute to the observed horizontal superiority while the broader distribution of vertical contours contributes to the observed vertical superiority in tolerance. These results indicate that natural scene statistics plays a significant role in determining the filling-in performance at the blind spot and shaping the associated anisotropies.


title: A quasi-periodic modulation of the iron line centroid energy in the black hole binary H 1743-322

abstract: Accreting stellar-mass black holes often show a `Type-C' quasi-periodic oscillation (QPO) in their X-ray flux, and an iron emission line in their X-ray spectrum. The iron line is generated through continuum photons reflecting off the accretion disk, and its shape is distorted by relativistic motion of the orbiting plasma and the gravitational pull of the black hole. The physical origin of the QPO has long been debated, but is often attributed to Lense-Thirring precession, a General Relativistic effect causing the inner flow to precess as the spinning black hole twists up the surrounding space-time. This predicts a characteristic rocking of the iron line between red and blue shift as the receding and approaching sides of the disk are respectively illuminated. Here we report on XMM-Newton and NuSTAR observations of the black hole binary H 1743-322 in which the line energy varies systematically over the ~4 s QPO cycle (3.70 sigma significance), as predicted. This provides strong evidence that the QPO is produced by Lense-Thirring precession, constituting the first detection of this effect in the strong gravitation regime. There are however elements of our results harder to explain, with one section of data behaving differently to all the others. Our result enables the future application of tomographic techniques to map the inner regions of black hole accretion disks.


title: Initial Experiments on Learning-Based Randomized Bin-Picking Allowing Finger Contact with Neighboring Objects

abstract: This paper proposes a novel method for randomized bin-picking based on learning. When a two-fingered gripper tries to pick an object from the pile, a finger often contacts a neighboring object. Even if a finger contacts a neighboring object, the target object will be successfully picked depending on the configuration of neighboring objects. In our proposed method, we use the visual information on neighboring objects to train the discriminator. Corresponding to a grasping posture of an object, the discriminator predicts whether or not the pick will be successful even if a finger contacts a neighboring object. We examine two learning algorithms, the linear support vector machine (SVM) and the random forest (RF) approaches. By using both methods, we demonstrate that the picking success rate is significantly higher than with conventional methods without learning.


title: Proceedings of the CAS - CERN Accelerator School: Power Converters, Baden, Switzerland, 7-14 May 2014

abstract: These proceedings collate lectures given at the twenty-eighth specialized course organised by the CERN Accelerator School (CAS). The course was held at the Hotel du Parc, Baden, Switzerland from 7 - 14 May 2014, in collaboration with the Paul Scherrer Institute. Following introductory lectures on accelerators and the requirements on power converters, the course covered components and topologies of the different types of power converters needed for particle accelerators. Issues of design, control and exploitation in a sometimes-hostile environment were addressed. Site visits to ABB and PSI provided an insight into state-of-the-art power converter production and operation, while topical seminars completed the programme.


title: Some Features of Scattering Problem in a $\kappa$-Deformed Minkowski Spacetime

abstract: The doubly special relativity (DSR) theories are suggested in order to incorporate an observer-independent length scale in special theory of relativity. The Magueijo-Smolin proposal of DSR is realizable through a particular form of the noncommutative (NC) spacetime (known as $\kappa$-Minkowski spacetime) in which the Lorentz symmetry is preserved. In this framework, the NC parameter $\kappa$ provides the origin of natural cutoff energy scale. Using a nonlinear deformed relativistic dispersion relation along with the Lorentz transformations, we investigate some phenomenological facets of two-body collision problem (without creation of new particles) in a $\kappa$-Minkowski spacetime. By treating an elastic scattering problem, we study effects of the Planck scale energy cutoff on some relativistic kinematical properties of this scattering problem. The results are challenging in the sense that as soon as one turns on the $\kappa$-spacetime extension, the nature of the two-body collision alters from elastic to inelastic one. It is shown also that a significant kinematical variable involving in heavy ion collisions, the rapidity, is not essentially an additive quantity under a sequence of the nonlinear representation of the Lorentz transformations.


title: Geomagnetically trapped, albedo and solar energetic particles: trajectory analysis and flux reconstruction with PAMELA

abstract: The PAMELA satellite experiment is providing comprehensive observations of the interplanetary and magnetospheric radiation in the near-Earth environment. Thanks to its identification capabilities and the semi-polar orbit, PAMELA is able to precisely measure the energetic spectra and the angular distributions of the different cosmic-ray populations over a wide latitude region, including geomagnetically trapped and albedo particles. Its observations comprise the solar energetic particle events between solar cycles 23 and 24, and the geomagnetic cutoff variations during magnetospheric storms. PAMELA's measurements are supported by an accurate analysis of particle trajectories in the Earth's magnetosphere based on a realistic geomagnetic field modeling, which allows the classification of particle populations of different origin and the investigation of the asymptotic directions of arrival.


title: Dirac Delta Function of Matrix Argument

abstract: Dirac delta function of matrix argument is employed frequently in the development of diverse fields such as Random Matrix Theory, Quantum Information Theory, etc. The purpose of the article is pedagogical, it begins by recalling detailed knowledge about Heaviside unit step function and Dirac delta function. Then its extensions of Dirac delta function to vector spaces and matrix spaces are discussed systematically, respectively. The detailed and elementary proofs of these results are provided. Though we have not seen these results formulated in the literature, there certainly are predecessors. Applications are also mentioned.


title: Two rate periodic protocol with dynamics driven through many cycles

abstract: We study the long time dynamics in closed quantum systems periodically driven via time dependent parameters with two frequencies $\omega_1$ and $\omega_2=r \omega_1$. Tuning of the ratio $r$ there can unleash plenty of dynamical phenomena to occur. Our study includes integrable models like Ising and XY models in $d=1$ and Kitaev model in $d=1$ and $2$ and can also be extended to Dirac fermions in graphene. We witness the wave-function overlap or dynamic freezing to occur within some small/ intermediate frequency regimes in the $(\omega_1 ,r)$ plane (with $r\ne0$) when the ground state is evolved through single cycle of driving. However, evolved states soon become steady with long driving and the freezing scenario gets rarer. We extend the formalism of adiabatic-impulse approximation for many cycle driving within our two-rate protocol and show the near-exact comparisons at small frequencies. An extension of the rotating wave approximation is also developed to gather an analytical framework of the dynamics at high frequencies. Finally we compute the entanglement entropy in the stroboscopically evolved states within the gapped phases of the system and observe how it gets tuned with the ratio $r$ in our protocol. The minimally entangled states are found to fall within the regime of dynamical freezing. In general, the results indicate that the entanglement entropy in our driven short-ranged integrable systems follow genuine non-area law of scaling and show a convergence (with a $r$ dependent pace) towards volume scaling behavior as the driving is continued for long time.


title: Deformations of Legendrian Curves

abstract: We construct versal and equimultiple versal deformations of the parametrization of a Legendrian curve.


title: Mapping Reactive Flow Patterns in Monolithic Nanoporous Catalysts

abstract: The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nano-scale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domain. To address this problem, we developed a multi-scale mass transport computational framework based on the Lattice Boltzmann Method (LBM) that allows us to account for catalytic reactions at the gas-solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane.

categories: physics.flu-dyn cond-mat.mtrl-sci physics.chem-ph physics.comp-ph


title: Radio-Optical Galaxy Shape Correlations in the COSMOS Field

abstract: We investigate the correlations in galaxy shapes between optical and radio wavelengths using archival observations of the COSMOS field. Cross-correlation studies between different wavebands will become increasingly important for precision cosmology as future large surveys may be dominated by systematic rather than statistical errors. In the case of weak lensing, galaxy shapes must be measured to extraordinary accuracy (shear systematics of $< 0.01\%$) in order to achieve good constraints on dark energy parameters. By using shape information from overlapping surveys in optical and radio bands, robustness to systematics may be significantly improved without loss of constraining power. Here we use HST-ACS optical data, VLA radio data, and extensive simulations to investigate both our ability to make precision measurements of source shapes from realistic radio data, and to constrain the intrinsic astrophysical scatter between the shapes of galaxies as measured in the optical and radio wavebands. By producing a new image from the VLA-COSMOS L-band radio visibility data that is well suited to galaxy shape measurements, we are able to extract precise measurements of galaxy position angles. Comparing to corresponding measurements from the HST optical image, we set a lower limit on the intrinsic astrophysical scatter in position angles, between the optical and radio bands, of $\sigma_\alpha > 0.212\pi$ radians (or $38.2^{\circ}$) at a $95\%$ confidence level.


title: Studying the evolution of AGB stars in the Gaia epoch

abstract: We present asymptotic giant branch (AGB) models of solar metallicity, to allow the interpretation of observations of Galactic AGB stars, whose distances should be soon available after the first release of the Gaia catalogue. We find an abrupt change in the AGB physical and chemical properties, occurring at the threshold mass to ignite hot bottom burning,i.e. $3.5M_{\odot}$. Stars with mass below $3.5 M_{\odot}$ reach the C-star stage and eject into the interstellar medium gas enriched in carbon , nitrogen and $^{17}O$. The higher mass counterparts evolve at large luminosities, between $3\times 10^4 L_{\odot}$ and $10^5 L_{\odot}$. The mass expelled from the massive AGB stars shows the imprinting of proton-capture nucleosynthesis, with considerable production of nitrogen and sodium and destruction of $^{12}C$ and $^{18}O$. The comparison with the most recent results from other research groups are discussed, to evaluate the robustness of the present findings. Finally, we compare the models with recent observations of galactic AGB stars, outlining the possibility offered by Gaia to shed new light on the evolution properties of this class of objects.


title: Dynamical symmetry enhancement near N=2, D=4 gauged supergravity horizons

abstract: We show that all smooth Killing horizons with compact horizon sections of 4-dimensional gauged N=2 supergravity coupled to any number of vector multiplets preserve $2 c_1({\cal K})+4 \ell$ supersymmetries, where ${\cal K}$ is a pull-back of the Hodge bundle of the special K\"ahler manifold on the horizon spatial section. We also demonstrate that all such horizons with $c_1({\cal K})=0$ exhibit an SL(2,R) symmetry and preserve either 4 or 8 supersymmetries. If the orbits of the SL(2,R) symmetry are 2-dimensional, the horizons are warped products of AdS2 with the horizon spatial section. Otherwise, the horizon section admits an isometry which preserves all the fields. The proof of these results is centered on the use of index theorem in conjunction with an appropriate generalization of the Lichnerowicz theorem for horizons that preserve at least one supersymmetry. In all $c_1({\cal K})=0$ cases, we specify the local geometry of spatial horizon sections and demonstrate that the solutions are determined by first order non-linear ordinary differential equations on some of the fields.


title: A duality theory for non-convex problems in the Calculus of Variations

abstract: We present a new duality theory for non-convex variational problems, under possibly mixed Dirichlet and Neumann boundary conditions. The dual problem reads nicely as a linear programming problem, and our main result states that there is no duality gap. Further, we provide necessary and sufficient optimality conditions, and we show that our duality principle can be reformulated as a min-max result which is quite useful for numerical implementations. As an example, we illustrate the application of our method to a celebrated free boundary problem. The results were announced in \cite{BoFr}.


title: A mathematical evaluation of vote transfer systems

abstract: The paper builds a general model of vote transfer systems on the basis of the current Hungarian electoral rules. It combines single-seat districts and list mandates with three possible compensation rule for 'wasted' votes in constituencies: no compensation (direct vote transfer, DVT), compensation for votes cast for losing party candidates (positive vote transfer, PVT) and compensation for all votes that are not necessary to win the district (negative vote transfer, NVT). The model is studied in the case of two parties. When the number of votes for the majority party follows a uniform distribution in each district, DVT results in the greatest expected seat share, however, application of PVT, and, especially, NVT increases the probability of winning the election. The trade-off between vote transfer formulas and the number of list mandates reveals that the majority party should use an appropriately calibrated NVT system if it focuses on these two variables.


title: In Cohen generic extension, every countable OD set of reals belongs to the ground model

abstract: It is true in the Cohen generic extension of L, the constructible universe, that every countable ordinal-definable set of reals belongs to L.


title: X-ray spectral and optical properties of a ULX in NGC 4258 (M106)

abstract: We study the X-ray and optical properties of the ultraluminous X-ray source (ULX) X-6 in the nearby galaxy NGC 4258 (M106) based on the archival {\it XMM-Newton}, {\it Chandra}, {\it Swift}, and {\it Hubble Space Telescope} ({\it HST}) observations. The source has a peak luminosity of $L_{\mathrm{X}} \sim 2 \times 10^{39}$ erg s$^{-1}$ in the {\it XMM-Newton} observation of 2004 June. Consideration of the hardness ratios and spectral model parameters shows that the source seems to exhibit possible spectral variations throughout the X-ray observations. In the images from the {\it HST}/Advanced Camera for Surveys (ACS), three optical sources have been identified as counterpart candidates within the 1$\sigma$ error radius of 0$\arcsec$.3. The brightest one has an absolute magnitude of $M_{V} \approx$ $-$7.0 and shows extended structure. The remaining two sources have absolute magnitudes of $M_{V} \approx$ $-$5.8 and $-$5.3 mag. The possible spectral types of the candidates from brightest to dimmest were determined as B6$-$A5, B0$-$A7, and B2$-$A3, respectively. The counterparts of the X-ray source possibly belong to a young star cluster. Neither the standard disk model nor the slim disk model provides firm evidence to determine the spectral characteristics of ULX X-6. We argue that the mass of the compact object lies in the range $10-15M_{\sun}$ indicating that the compact source is most likely a stellar-mass black hole.


title: Existence of torsional solitons in a beam model of suspension bridge

abstract: This paper studies the existence of solitons, namely stable solitary waves, in a suspension bridge. The bridge is modeled as a degenerate plate, that is, a central beam with cross sections, and displays two degrees of freedom: the vertical displacement of the beam and the torsional angles of the cross sections. Under fairy general assumptions, we prove the existence of solitons. Under the additional assumption of large tension in the sustaining cables, we prove that these solitons have a nontrivial torsional component. This appears relevant for the security since several suspension bridges collapsed due to torsional oscillations.


title: Non-Concave Penalization in Linear Mixed-Effects Models and Regularized Selection of Fixed Effects

abstract: Mixed-effect models are very popular for analyzing data with a hierarchical structure, e.g. repeated observations within subjects in a longitudinal design, patients nested within centers in a multicenter design. However, recently, due to the medical advances, the number of fixed effect covariates collected from each patient can be quite large, e.g. data on gene expressions of each patient, and all of these variables are not necessarily important for the outcome. So, it is very important to choose the relevant covariates correctly for obtaining the optimal inference for the overall study. On the other hand, the relevant random effects will often be low-dimensional and pre-specified. In this paper, we consider regularized selection of important fixed effect variables in linear mixed-effects models along with maximum penalized likelihood estimation of both fixed and random effect parameters based on general non-concave penalties. Asymptotic and variable selection consistency with oracle properties are proved for low-dimensional cases as well as for high-dimensionality of non-polynomial order of sample size (number of parameters is much larger than sample size). We also provide a suitable computationally efficient algorithm for implementation. Additionally, all the theoretical results are proved for a general non-convex optimization problem that applies to several important situations well beyond the mixed model set-up (like finite mixture of regressions etc.) illustrating the huge range of applicability of our proposal.


title: RemoveYoung: A tool for the removal of the young stellar component in galaxies within an adjustable age cutoff

abstract: The optical morphology of galaxies holds the cumulative record of their assembly history, and techniques for its quantitative characterization offer a promising avenue toward understanding galaxy formation and evolution. However, the morphology of star-forming galaxies is generally dictated by the youngest stellar component, which can readily overshine faint structural/morphological features in the older underlying stellar background (e.g., relics from recent minor mergers) that could hold important insights into the galaxy build-up process. Stripping off galaxy images from the emission from stellar populations younger than an adjustable age cutoff tcut can, therefore, provide a valuable tool in extragalactic research. RemoveYoung (RY), a publicly available tool that is presented here, exploits the combined power of integral field spectroscopy (IFS) and spectral population synthesis (SPS) toward this goal. Two-dimensional (2D) post-processing of SPS models to IFS data cubes with RY permits computation of the spectral energy, surface brightness, and stellar surface density distribution of stellar populations older than a user-defined tcut. This suggests a variety of applications of star-forming galaxies, such as interacting or merging galaxy pairs and lower mass starburst galaxies near and far; these include blue compact and tidal dwarf galaxies.


title: The Wulff construction for convex integrands

abstract: For any given Wulff shape $\mathcal{W}$, we can define the unique continuous function $S^{n}\to \mathbb{R}_{+}$ called convex integrand, denoted by $\gamma_{{}_{\mathcal{W}}}$. In this paper, we show that, for any Wulff shapes $\mathcal{W}_{1}$ and $\mathcal{W}_{2}$, the equality $d(\gamma_{{}_{\mathcal{W}_{1}}}, \gamma_{{}_{\mathcal{W}_{2}}})= h(\mathcal{W}_{1}, \mathcal{W}_{2})$ holds, where $d$ is the maximum distance of the function space consisting of convex integrands and $h$ is the Pompeiu-Hausdorff distance of the space consisting of Wulff shapes. Moreover, applications of this result are given.


title: Reply to Comment on "Tip induced unconventional superconductivity on Weyl semimetal TaAs"

abstract: Recently, we presented a paper about the tip-induced superconducting phase on TaAs single crystal (arXiv:1607.00513). A conductance plateau with double conductance peaks at finite bias, sharp double dips and a zero bias conductance peak were observed in the point contact spectra (PCS). These three features in one PCS suggest the possibility of p-wave like superconductivity and Majorana zero modes in the surface. Soon after our paper posted on arXiv, Gayen el al. (arXiv:1607.01405) questioned the underlying mechanism of our observations based on their experimental measurements in conventional superconductors(Pb and Nb) and conclusively indicated that the superconductivity observed in our TaAs crystal is conventional. We disagree with their quick conclusion without any specific study on the TaAs materials. We will explain in the following that for the clean (Z=0) point contact (PC) in ballistic limit, the superconducting transition feature in temperature dependence of PC resistance is expected. And the critical current effect claimed by Gayen el al. cannot explain the observed PCS in our paper.


title: On the growth of the Kronecker coefficients

abstract: We study the rate of growth experienced by the Kronecker coefficients as we add cells to the rows and columns indexing partitions. We do this by moving to the setting of the reduced Kronecker coefficients.


title: Coverings: variations on a result of Rogers and on the Epsilon-net theorem of Haussler and Welzl

abstract: We consider four problems. Rogers proved that for any convex body $K$, we can cover ${\mathbb R}^d$ by translates of $K$ of density very roughly $d\ln d$. First, we extend this result by showing that, if we are given a family of positive homothets of $K$ of infinite total volume, then we can find appropriate translation vectors for each given homothet to cover ${\mathbb R}^d$ with the same (or, in certain cases, smaller) density. Second, we extend Rogers' result to multiple coverings of space by translates of a convex body: we give a non-trivial upper bound on the density of the most economical covering where each point is covered by at least a certain number of translates. Third, we show that for any sufficiently large $n$, the sphere ${\mathbb S}^2$ can be covered by $n$ strips of width $20n/\ln n$, where no point is covered too many times. Finally, we give another proof of the previous result based on a combinatorial observation: an extension of the Epsilon-net Theorem of Haussler and Welzl. We show that for a hypergraph of bounded Vapnik--Chervonenkis dimension, in which each edge is of a certain measure, there is a not-too large transversal set which does not intersect any edge too many times.


title: Nonlinear embeddings: Applications to analysis, fractals and polynomial root finding

abstract: We introduce $\mathcal{B}_{\kappa}$-embeddings, nonlinear mathematical structures that connect, through smooth paths parameterized by $\kappa$, a finite or denumerable set of objects at $\kappa=0$ (e.g. numbers, functions, vectors, coefficients of a generating function...) to their ordinary sum at $\kappa \to \infty$. We show that $\mathcal{B}_{\kappa}$-embeddings can be used to design nonlinear irreversible processes through this connection. A number of examples of increasing complexity are worked out to illustrate the possibilities uncovered by this concept. These include not only smooth functions but also fractals on the real line and on the complex plane. As an application, we use $\mathcal{B}_{\kappa}$-embeddings to formulate a robust method for finding all roots of a univariate polynomial without factorizing or deflating the polynomial. We illustrate this method by finding all roots of a polynomial of 19th degree.


title: Note on three-generation models in heterotic string and F-theory on elliptic Calabi-Yau manifolds over Hirzebruch varieties

abstract: We give a complete list of a class of three-generation models in E8 x E8 heterotic string theory and its dual F-theory on an elliptic Calabi-Yau over a (generalized) Hirzebruch variety in which the divisors of the relevant line bundles needed for a smooth Weierstrass model are effective. The most stringent constraint on the bound of the eta class comes from the effectiveness of the divisor of the bundle corresponding to the highest Casimir in Looijenga's weighted projective space, as well as from the compactness of the toric variety. Comparison is also made with the list obtained in the literature.


title: Effects of an induced electric field on \pi^{-}/\pi^{+} ratio in heavy-ion collisions

abstract: Using an isospin- and momentum-dependent transport model, we examine the effects of an electric field induced by a variable magnetic field on the \pi^{-}/\pi^{+} ratio in central to peripheral heavy-ion collisions at beam energies of 400 and 1500MeV/nucleon. It is shown that while the induced electric field does not affect the total multiplicities of both $\pi^{-}$ and $\pi^{+}$ mesons at both the lower beam energy of 400MeV/nucleon and the higher beam energy of 1500MeV/nucleon, it reduces (enhances) the emission of $\pi^{-}$ ($\pi^{+}$) mesons in midrapidity, but enhances (reduces) the emission of $\pi^{-}$ ($\pi^{+}$) mesons in forward and backward rapidities especially for the more peripheral collisions at the lower beam energy because of the rapidly transient variable magnetic field at more peripheral collisions and longer reaction duration time at the lower beam energy. These findings indicate that the total \pi^{-}/\pi^{+} ratio is still a precisely reliable probe of symmetry energy at both the lower and higher beam energies, but one should consider the induced electric field when using the differential \pi^{-}/\pi^{+} ratio to probe the symmetry energy especially for the lower beam energy and more peripheral collisions. Finally, the relative suppression factor based on the ratio of \pi^{-}/\pi^{+} in different rapidities is proposed to be an effective probe of the induced electric field generated in heavy-ion collisions due to its maximizing effects of induced electric fields on the differential \pi^{-}/\pi^{+} ratio but minimizing effects of some uncertainty factors in heavy-ion collisions.


title: Filter-feeding, near-field flows, and the morphologies of colonial choanoflagellates

abstract: Efficient uptake of nutrients from the environment is an important component in the fitness of all microorganisms, and its dependence on size may reveal clues to the origins of evolutionary transitions to multicellularity. Because potential benefits in uptake rates must be viewed in the context of other costs and benefits of size, such as varying predation rates and the increased metabolic costs associated with larger and more complex body plans, the uptake rate itself is not necessarily that which is optimized by evolution. Uptake rates can be strongly dependent on local organism geometry and its swimming speed, providing selective pressure for particular arrangements. Here we examine these issues for choanoflagellates, filter-feeding microorganisms that are the closest relatives of the animals. We explore the different morphological variations of the choanoflagellete $Salpingoeca~rosetta$, which can exist as a swimming cell, a sessile thecate cell, and as colonies of cells in various shapes. In the absence of other requirements and in a homogeneously nutritious environment, we find that the optimal strategy to maximize filter feeding by the collar of microvilli is to swim fast, which favours swimming unicells. In large external flows, the sessile thecate cell becomes advantageous. Effects of prey diffusion are discussed and also found to be to the advantage of the swimming unicell.


title: Deterministic Annealing Optimization for Witsenhausen's and Related Decentralized Stochastic Control Problems

abstract: This note studies the global optimization of controller mappings in discrete-time stochastic control problems including Witsenhausen's celebrated 1968 counter-example. We propose a generally applicable non-convex numerical optimization method based on the concept of deterministic annealing-which is derived from information-theoretic principles and was successfully employed in several problems including vector quantization, classification, and regression. We present comparative numerical results for two test problems that show the strict superiority of the proposed method over prior approaches in the literature.


title: A new MC-based method to evaluate the fission fraction uncertainty at reactor neutrino experiment

abstract: Uncertainties of fission fraction is an important uncertainty source for the antineutrino flux prediction in a reactor antineutrino experiment. A new MC-based method of evaluating the covariance coefficients between isotopes was proposed. It was found that the covariance coefficients will varying with reactor burnup and which may change from positive to negative because of fissioning balance effect, for example, the covariance coefficient between $^{235}$U and $^{239}$Pu changes from 0.15 to -0.13. Using the equation between fission fraction and atomic density, the consistent of uncertainty of fission fraction and the covariance matrix were obtained. The antineutrino flux uncertainty is 0.55\% which does not vary with reactor burnup, and the new value is about 8.3\% smaller.


title: Lagrangian Decomposition based Multi Agent Model Predictive Control for Electric Vehicles Charging integrating Real Time Pricing

abstract: This paper presents a real time distributed control strategy for electric vehicles charging covering both drivers and grid players' needs. Computation of the charging load curve is performed by agents working at the level of each single vehicle, with the information exchanged with grid players being restricted to the chosen load curve and energy price feedback from the market, elaborated according to the charging infrastructure congestion. The distributed control mechanism is based on model predictive control methodology and Lagrangian decomposition of the optimization control problem at its basis. The simulation results show the effectiveness of the proposed distributed approach and the mutual coherence between the computed charging load curves and the resulting energy price over the time.


title: Conjugacy properties of time-evolving Dirichlet and gamma random measures

abstract: We extend classic characterisations of posterior distributions under Dirichlet process and gamma random measures priors to a dynamic framework. We consider the problem of learning, from indirect observations, two families of time-dependent processes of interest in Bayesian nonparametrics: the first is a dependent Dirichlet process driven by a Fleming-Viot model, and the data are random samples from the process state at discrete times; the second is a collection of dependent gamma random measures driven by a Dawson-Watanabe model, and the data are collected according to a Poisson point process with intensity given by the process state at discrete times. Both driving processes are diffusions taking values in the space of discrete measures whose support varies with time, and are stationary and reversible with respect to Dirichlet and gamma priors respectively. A common methodology is developed to obtain in closed form the time-marginal posteriors given past and present data. These are shown to belong to classes of finite mixtures of Dirichlet processes and gamma random measures for the two models respectively, yielding conjugacy of these classes to the type of data we consider. We provide explicit results on the parameters of the mixture components and on the mixing weights, which are time-varying and drive the mixtures towards the respective priors in absence of further data. Explicit algorithms are provided to recursively compute the parameters of the mixtures. Our results are based on the projective properties of the signals and on certain duality properties of their projections.


title: Continuous measurements of real-life bidirectional pedestrian flows on a wide walkway

abstract: Employing partially overlapping overhead \kinectTMS sensors and automatic pedestrian tracking algorithms we recorded the crowd traffic in a rectilinear section of the main walkway of Eindhoven train station on a 24/7 basis. Beside giving access to the train platforms (it passes underneath the railways), the walkway plays an important connection role in the city. Several crowding scenarios occur during the day, including high- and low-density dynamics in uni- and bi-directional regimes. In this paper we discuss our recording technique and we illustrate preliminary data analyses. Via fundamental diagrams-like representations we report pedestrian velocities and fluxes vs. pedestrian density. Considering the density range $0$ - $1.1\,$ped/m$^2$, we find that at densities lower than $0.8\,$ped/m$^2$ pedestrians in unidirectional flows walk faster than in bidirectional regimes. On the opposite, velocities and fluxes for even bidirectional flows are higher above $0.8\,$ped/m$^2$.


title: Phase dependent advection-diffusion in drift wave - zonal flow turbulence

abstract: In plasma turbulence theory, due to the complexity of the system with many non-linearly interacting waves, the dynamics of the phases is often disregarded and the so-called random-phase approximation (RPA) is used assuming the existence of a Chirikov-like criterion for the onset of wave stochasticity. The dynamical amplitudes are represented as complex numbers, $\psi = \psi_r + i\psi_i = ae^{i\theta}$, with the amplitudes slowly varying whereas the phases are rapidly varying and, in particular, distributed uniformly over the interval $[0;2\pi)$. However, one could expect that the phase dynamics can play a role in the self-organisation and the formation of coherent structures. In the same manner it is also expected that the RPA falls short to take coherent interaction between phases into account. In this work therefore, we studied the role of phase dynamics and the coupling of phases between different modes on the characteristic time evolution of the turbulent. We assume a simple turbulent system where the so-called stochastic oscillator model can be employed. The idea of interpreting turbulence by stochastic oscillators. The stochastic oscillator models can be derived from radical simplifications of the nonlinear terms in the Navier-Stokes or Gyro-Kinetic equations. In this particular case we adopt the basic equation for the stochastic oscillator model with passive advection and random forcing from Ref.


title: Collapse of the mean curvature flow for certain kind of invariant hypersurfaces in a Hilbert space

abstract: In this paper, we investigate the regularized mean curvature flow starting from an invariant hypersurface in a Hilbert space equipped with an isometric and almost free action of a Hilbert Lie group whose orbits are regularized minimal. We prove that, if the invariant hypersurface satisfies a certain kind of horizontally convexity condition, then it collapses to an orbit of the Hilbert Lie group action along the regularized mean curvature flow.


title: Reaction and fusion cross sections for the near-symmetric system $^{129}Xe+^{nat}Sn$ from $8$ to $35$ $AMeV$

abstract: \item[Background]Heavy-ion reactions from barrier up to Fermi energy. \item[Purpose]Reaction and fusion cross sections determination. Fusion reactions induced by $^{129}Xe$ projectiles on $^{nat}Sn$ targets for energies ranging from $8$ A.MeV to $35$ A.MeV were measured with the INDRA $4\pi$-array.\\ The evaluation of the fusion/incomplete fusion cross sections for the incident energies from 8 to 35 A.MeV is the main purpose of this paper. \item[Method] The reaction cross sections are evaluated for each beam energy thanks to INDRA $4\pi$-array. The events are also sorted in order to focus the study on a selected sample of events, in such a way that the fusion/fusion incomplete cross section is estimated. \item[Results] The excitation function of reaction and fusion cross sections were measured for the heavy and nearly symmetric system $^{129}Xe + ^{nat}Sn$ from 8 to 35 A.MeV. \item[Conclusions] The fusion-like cross-sections evaluated show a good agrement with a recent systematics for beam energies greater than 20 A.MeV. For low beam energies the cross-section values are lower than the expected ones. A probable reason for these low values is in the fusion hindrance at energies above/close the barrier.


title: The electronic structure of BSCCO in the presence of a super-current: Flux-flow, Doppler shift and quasiparticle pockets

abstract: There are several ways to turn a superconductor into a normal conductor: increase the temperature, apply a high magnetic field, or run a large current. High-T$_c$ cuprate superconductors are unusual in the sense that experiments suggest that destroying superconductivity by heating the sample to temperatures above T$_c$ or by applying a high magnetic field result in different 'normal' states. Spectroscopic probes show that above T$_c$, in the pseudogap regime, the Fermi surface is partly gapped and there are no well-defined quasiparticles. Transport measurements, on the contrary, reveal quantum oscillations in high magnetic fields and at low temperatures, suggesting a more usual Fermi liquid state. Studying the electronic structure while suppressing superconductivity by using current, will hopefully shed new light on this problem. In type II superconductors, such as the cuprates, the resistive state created by the current is a result of vortex motion. The large dissipation in the flux-flow regime makes spectroscopic measurements in the presence of current challenging. We performed angle-resolved photoemission experiments in thin films of Bi$_2$Sr$_2$CaCu$_2$O$_y$ while running high-density current through the samples. Clear evidence was found for non-uniform flux flow, leaving most of the sample volume free of mobile vortices and dissipation. The super-current changes the electronic spectrum, creating quasiparticle and quasihole pockets. The size of these pockets as a function of the current is found to be doping dependent; it depends both on the superfluid stiffness and on the strength of interactions.


title: sk_p: a neural program corrector for MOOCs

abstract: We present a novel technique for automatic program correction in MOOCs, capable of fixing both syntactic and semantic errors without manual, problem specific correction strategies. Given an incorrect student program, it generates candidate programs from a distribution of likely corrections, and checks each candidate for correctness against a test suite. The key observation is that in MOOCs many programs share similar code fragments, and the seq2seq neural network model, used in the natural-language processing task of machine translation, can be modified and trained to recover these fragments. Experiment shows our scheme can correct 29% of all incorrect submissions and out-performs state of the art approach which requires manual, problem specific correction strategies.


title: Magnetorotational dynamo chimeras. The missing link to turbulent accretion disk dynamo models?

abstract: In Keplerian accretion disks, turbulence and magnetic fields may be jointly excited through a subcritical dynamo process involving the magnetorotational instability (MRI). High-resolution simulations exhibit a tendency towards statistical self-organization of MRI dynamo turbulence into large-scale cyclic dynamics. Understanding the physical origin of these structures, and whether they can be sustained and transport angular momentum efficiently in astrophysical conditions, represents a significant theoretical challenge. The discovery of simple periodic nonlinear MRI dynamo solutions has recently proven useful in this respect, and has notably served to highlight the role of turbulent magnetic diffusion in the seeming decay of the dynamics at low magnetic Prandtl number Pm (magnetic diffusivity larger than viscosity), a common regime in accretion disks. The connection between these simple structures and the statistical organization reported in turbulent simulations remained elusive, though. Here, we report the numerical discovery in moderate aspect ratio Keplerian shearing boxes of new periodic, incompressible, three-dimensional nonlinear MRI dynamo solutions with a larger dynamical complexity reminiscent of such simulations. These "chimera" cycles are characterized by multiple MRI-unstable dynamical stages, but their basic physical principles of self-sustainment are nevertheless identical to those of simpler cycles found in azimuthally elongated boxes. In particular, we find that they are not sustained at low Pm either due to subcritical turbulent magnetic diffusion. These solutions offer a new perspective into the transition from laminar to turbulent instability-driven dynamos, and may prove useful to devise improved statistical models of turbulent accretion disk dynamos.


title: The Vectorization of the Tersoff Multi-Body Potential: An Exercise in Performance Portability

abstract: Molecular dynamics simulations, an indispensable research tool in computational chemistry and materials science, consume a significant portion of the supercomputing cycles around the world. We focus on multi-body potentials and aim at achieving performance portability. Compared with well-studied pair potentials, multibody potentials deliver increased simulation accuracy but are too complex for effective compiler optimization. Because of this, achieving cross-platform performance remains an open question. By abstracting from target architecture and computing precision, we develop a vectorization scheme applicable to both CPUs and accelerators. We present results for the Tersoff potential within the molecular dynamics code LAMMPS on several architectures, demonstrating efficiency gains not only for computational kernels, but also for large-scale simulations. On a cluster of Intel Xeon Phi's, our optimized solver is between 3 and 5 times faster than the pure MPI reference.


title: High pressure studies on properties of FeGa3: role of on-site coulomb correlation

abstract: High pressure X-ray diffraction measurements have been carried out on the intermetallic semiconductor FeGa$_3$ and the equation of state for FeGa$_3$ has been determined. First principles based DFT calculations within the GGA approximation indicate that although the unit cell volume matches well with the experimentally obtained value at ambient pressure, it is significantly underestimated at high pressures and the difference between them increases as pressure increases. GGA + U calculations with increasing values of U$_{Fe(3d)}$ (on-site Coulomb repulsion between the Fe 3d electrons) at high pressures, correct this discrepancy. Further, the GGA+U calculations also show that along with U$_{Fe(3d)}$, the Fe 3d band width also increases with pressure and around a pressure of 4 GPa, a small density of states appear at the Fermi level. High pressure resistance measurements carried out on FeGa$_3$ also clearly show a signature of an electronic transition. Beyond the pressure of 19.7 GPa, the diffraction peaks reduce in intensity and are not observable beyond $\sim$ 26 GPa, leading to an amorphous state.


title: On the amplification of acoustic phonons in carbon nanotube

abstract: We present a theoretical study of acoustic phonons amplification in Carbon Nanotubes (CNT). The phenomenon is via Cerenkov emission (CE) of acoustic phonons using intraband transitions proposed by Mensah et. al.,~\cite{1} in Semiconductor Superlattices (SSL) and confirmed in ~\cite{2}. From this, an asymmetric graph of $\Gamma^{CNT}$ on $\frac{V_d}{V_s}$ and $\Omega\tau$ were obtained where amplification ($\Gamma_{amp}^{CNT}$) $>>$ absorption ($\Gamma_{abs}^{CNT}$). The ratio, $\frac{\vert \Gamma_{amp}^{CNT}\vert}{\vert\Gamma_{abs}^{CNT}\vert}\approx 3.5$, at $V_d = 1.02V_s$, $\omega_q = 3.0\ \mathrm{THz}$ and $T = 85\ K$ for scattering angle $\theta > 0$ . A threshold field at which $\Gamma_{abs}^{CNT}$ switches over to $\Gamma_{amp}^{CNT}$ was calculated to be $E_{z}^{dc} = 6.2\times 10^3\ \mathrm{V/m}$. This field is far less than that deduced using Bloch-Type Oscillation (BTO)~\cite{3} which is $E_{BTO}^{dc} = 3.0\times 10^5\ \mathrm{V/m}$. The obtained $\Gamma_{amp}^{CNT}$ would enable the use of CNT for the production of SASER.


title: On contact and symplectic Lie algebroids

abstract: In this paper, we will study compatible triples on Lie algebroids. Using a suitable decomposition for a Lie algebroid, we construct an integrable generalized distribution on the base manifold. As a result, the symplectic form on the Lie algebroid induces a symplectic form on each integral submanifold of the distribution. The induced Poisson structure on the base manifold can be represented by means of the induced Poisson structures on the integral submanifolds. Moreover, for any compatible triple with invariant metric and admissible almost complex structure, we show that the bracket annihilates on the kernel of the anchor map.