title: Classification of geometrical objects by integrating currents and functional data analysis. An application to a 3D database of Spanish child population

abstract: This paper focuses on the application of Discriminant Analysis to a set of geometrical objects (bodies) characterized by currents. A current is a relevant mathematical object to model geometrical data, like hypersurfaces, through integration of vector fields along them. As a consequence of the choice of a vector-valued Reproducing Kernel Hilbert Space (RKHS) as a test space to integrate hypersurfaces, it is possible to consider that hypersurfaces are embedded in this Hilbert space. This embedding enables us to consider classification algorithms of geometrical objects. A method to apply Functional Discriminant Analysis in the obtained vector-valued RKHS is given. This method is based on the eigenfunction decomposition of the kernel. So, the novelty of this paper is the reformulation of a size and shape classification problem in Functional Data Analysis terms using the theory of currents and vector-valued RKHS. This approach is applied to a 3D database obtained from an anthropometric survey of the Spanish child population with a potential application to online sales of children's wear.


title: Millimeter and Submillimeter Wave Spectroscopy of Propanal

abstract: The rotational spectra of the two stable conformers syn- and gauche-propanal (CH$_3$CH$_2$CHO) were studied in the millimeter and submillimeter wave regions from 75 to 500 GHz with the Cologne (Sub-)Millimeter wave Spectrometer. Furthermore, the first excited states associated with the aldehyde torsion and with the methyl torsion, respectively, of the $syn$-conformer were analyzed. The newly obtained spectroscopic parameters yield better predictions, thus fulfill sensitivity and resolution requirements in new astronomical observations in order to unambiguously assign pure rotational transitions of propanal. This is demonstrated on a radio astronomical spectrum from the Atacama Large Millimeter/submillimeter Array Protostellar Interferometric Line Survey (ALMA-PILS). In particular, an accurate description of observed splittings, caused by internal rotation of the methyl group in the $syn$-conformer and by tunneling rotation interaction from two stable degenerate $gauche$-conformers, is reported. The rotational spectrum of propanal is of additional interest because of its two large amplitude motions pertaining to the methyl and the aldehyde group, respectively.


title: A characterization of equivalent martingale measures in a renewal risk model with applications to premium calculation principles

abstract: Generalizing earlier work of Delbaen and Haezendonck for given compound renewal process $S$ under a probability measure $P$ we characterize all probability measures $Q$ on the domain of $P$ such that $Q$ and $P$ are progressively equivalent and $S$ remains a compound renewal process under $Q$. As a consequence, we prove that any compound renewal process can be converted into a compound Poisson process through a change of measures and we show how this approach is related to premium calculation principles.


title: Asymptotic behavior of 3-D stochastic primitive equations of large-scale moist atmosphere with additive noise

abstract: Using a new and general method, we prove the existence of random attractor for the three dimensional stochastic primitive equations defined on a manifold $\D\subset\R^3$ improving the existence of weak attractor for the deterministic model. Furthermore, we show the existence of the invariant measure.


title: Poincar\'e profiles of groups and spaces

abstract: We introduce a spectrum of monotone coarse invariants for metric measure spaces called Poincar\'{e} profiles. The two extremes of this spectrum determine the growth of the space, and the separation profile as defined by Benjamini--Schramm--Tim\'{a}r. In this paper we focus on properties of the Poincar\'{e} profiles of groups with polynomial growth, and of hyperbolic spaces, where we deduce a connection between these profiles and conformal dimension. As applications, we use these invariants to show the non-existence of coarse embeddings in a variety of examples.


title: Secure Symmetric Private Information Retrieval from Colluding Databases with Adversaries

abstract: The problem of symmetric private information retrieval (SPIR) from replicated databases with colluding servers and adversaries is studied. Specifically, the database comprises $K$ files, which are replicatively stored among $N$ servers. A user wants to retrieve one file from the database by communicating with the $N$ servers, without revealing the identity of the desired file to any server. Furthermore, the user shall learn nothing about the other $K-1$ files. Any $T$ out of $N$ servers may collude, that is, they may communicate their interactions with the user to guess the identity of the requested file. An adversary in the system can tap in on or even try to corrupt the communication. Three types of adversaries are considered: a Byzantine adversary who can overwrite the transmission of any $B$ servers to the user; a passive eavesdropper who can tap in on the incoming and outgoing transmissions of any $E$ servers; and a combination of both -- an adversary who can tap in on a set of any $E$ nodes, and overwrite the transmission of a set of any $B$ nodes. The problems of SPIR with colluding servers and the three types of adversaries are named T-BSPIR, T-ESPIR and T-BESPIR respectively. The capacity of the problem is defined as the maximum number of information bits of the desired file retrieved per downloaded bit. We show that the information-theoretical capacity of T-BSPIR equals $1-\frac{2B+T}{N}$, if the servers share common randomness (unavailable at the user) with amount at least $\frac{2B+T}{N-2B-T}$ times the file size. Otherwise, the capacity equals zero. The capacity of T-ESPIR is proved to equal $1-\frac{\max(T,E)}{N}$, with common randomness at least $\frac{\max(T,E)}{N-\max(T,E)}$ times the file size. Finally, the capacity of T-BESPIR is proved to be $1-\frac{2B+\max(T,E)}{N}$, with common randomness at least $\frac{2B+\max(T,E)}{N-2B-\max(T,E)}$ times the file size.


title: A Cut Discontinuous Galerkin Method for Coupled Bulk-Surface Problems

abstract: We develop a cut Discontinuous Galerkin method (cutDGM) for a diffusion-reaction equation in a bulk domain which is coupled to a corresponding equation on the boundary of the bulk domain. The bulk domain is embedded into a structured, unfitted background mesh. By adding certain stabilization terms to the discrete variational formulation of the coupled bulk-surface problem, the resulting cutDGM is provably stable and exhibits optimal convergence properties as demon- strated by numerical experiments. We also show both theoretically and numerically that the system matrix is well-conditioned, irrespective of the relative position of the bulk domain in the background mesh.


title: A Hybrid High-Order method for nonlinear elasticity

abstract: In this work we propose and analyze a novel Hybrid High-Order discretization of a class of (linear and) nonlinear elasticity models in the small deformation regime which are of common use in solid mechanics. The proposed method is valid in two and three space dimensions, it supports general meshes including polyhedral elements and nonmatching interfaces, enables arbitrary approximation order, and the resolution cost can be reduced by statically condensing a large subset of the unknowns for linearized versions of the problem. Additionally, the method satisfies a local principle of virtual work inside each mesh element, with interface tractions that obey the law of action and reaction. A complete analysis covering very general stress-strain laws is carried out, and optimal error estimates are proved. Extensive numerical validation on model test problems is also provided on two types of nonlinear models.


title: Q-systems and compact W*-algebra objects

abstract: We show that given a rigid C*-tensor category, there is an equivalence of categories between normalized irreducible Q-systems, also known as connected unitary Frobenius algebra objects, and compact connected W*-algebra objects. Although this result could be proved as a corollary of our previous article on realizations of algebra objects and discrete subfactors, we prove it here directly via categorical methods without passing through subfactor theory.


title: A highly accurate boundary integral equation method for surfactant-laden drops in 3D

abstract: The presence of surfactants alters the dynamics of viscous drops immersed in an ambient viscous fluid. This is specifically true at small scales, such as in applications of droplet based microfluidics, where the interface dynamics become of increased importance. At such small scales, viscous forces dominate and inertial effects are often negligible. Considering Stokes flow, a numerical method based on a boundary integral formulation is presented for simulating 3D drops covered by an insoluble surfactant. The method is able to simulate drops with different viscosities and close interactions, automatically controlling the time step size and maintaining high accuracy also when substantial drop deformation appears. To achieve this, the drop surfaces as well as the surfactant concentration on each surface are represented by spherical harmonics expansions. A novel reparameterization method is introduced to ensure a high-quality representation of the drops also under deformation, specialized quadrature methods for singular and nearly singular integrals that appear in the formulation are evoked and the adaptive time stepping scheme for the coupled drop and surfactant evolution is designed with a preconditioned implicit treatment of the surfactant diffusion.


title: Atmospheric turbulence profiling with unknown power spectral density

abstract: Adaptive optics (AO) is a technology in modern ground-based optical telescopes to compensate the wavefront distortions caused by atmospheric turbulence. One method that allows to retrieve information about the atmosphere from telescope data is so-called SLODAR, where the atmospheric turbulence profile is estimated based on correlation data of Shack--Hartmann wavefront measurements. This approach relies on a layered Kolmogorov turbulence model. In this article, we propose a novel extension of the SLODAR concept by including a general non-Kolmogorov turbulence layer close to the ground with an unknown power spectral density. We prove that the joint estimation problem of the turbulence profile above ground simultaneously with the unknown power spectral density at the ground is ill-posed and propose three numerical reconstruction methods. We demonstrate by numerical simulations that our methods lead to substantial improvements in the turbulence profile reconstruction, compared to standard SLODAR-type approach. Also, our methods can accurately locate local perturbations in non-Kolmogorov power spectral densities.


title: Deep Character-Level Click-Through Rate Prediction for Sponsored Search

abstract: Predicting the click-through rate of an advertisement is a critical component of online advertising platforms. In sponsored search, the click-through rate estimates the probability that a displayed advertisement is clicked by a user after she submits a query to the search engine. Commercial search engines typically rely on machine learning models trained with a large number of features to make such predictions. This is inevitably requires a lot of engineering efforts to define, compute, and select the appropriate features. In this paper, we propose two novel approaches (one working at character level and the other working at word level) that use deep convolutional neural networks to predict the click-through rate of a query-advertisement pair. Specially, the proposed architectures only consider the textual content appearing in a query-advertisement pair as input, and produce as output a click-through rate prediction. By comparing the character-level model with the word-level model, we show that language representation can be learnt from scratch at character level when trained on enough data. Through extensive experiments using billions of query-advertisement pairs of a popular commercial search engine, we demonstrate that both approaches significantly outperform a baseline model built on well-selected text features and a state-of-the-art word2vec-based approach. Finally, by combining the predictions of the deep models introduced in this study with the prediction of the model in production of the same commercial search engine, we significantly improve the accuracy and the calibration of the click-through rate prediction of the production system.


title: Interacting and noninteracting integrable systems

abstract: The distinction between interacting and noninteracting integrable systems is proposed to be characterized by the Onsager matrix. It being zero is the defining property of a noninteracting integrable system. To support our view we discuss various classical and quantum integrable chains.


title: TOPCAT: Desktop Exploration of Tabular Data for Astronomy and Beyond

abstract: TOPCAT, the Tool for OPerations on Catalogues And Tables, is an interactive desktop application for retrieval, analysis and manipulation of tabular data, offering a powerful and flexible range of interactive visualization options amongst other features. Its visualization capabilities focus on enabling interactive exploration of large static local tables - millions of rows and hundreds of columns can easily be handled on a standard desktop or laptop machine, and various options are provided for meaningful graphical representation of such large datasets. TOPCAT has been developed in the context of astronomy, but many of its features are equally applicable to other domains. The software, which is free and open source, is written in Java, and the underlying high-performance visualisation library is suitable for re-use in other applications.


title: Interplay between breathing mode distortion and magnetic order in rare-earth nickelates $R$NiO$_3$ within DFT+$U$

abstract: We present a systematic density functional theory (DFT) plus Hubbard $U$ study of structural trends and the stability of different magnetically ordered states across the rare-earth nickelate series, $R$NiO$_3$, with $R$ from Lu to La. In particular, we investigate how the magnetic order, the change of the rare-earth ion, and the Hubbard interaction $U$ are affecting the bond-length disproportionation between the nickel sites. Our results show that structural parameters can be obtained that are in very good agreement with present experimental data, and that DFT+$U$ is in principle able to capture the most important structural trends across the nickelate series. However, the amplitude of the bond-length disproportionation depends very strongly on the specific value used for the Hubbard $U$ parameter and also on the type of magnetic order imposed in the calculation. Regarding the relative stability of different magnetic orderings, a realistic antiferromagnetic order, consistent with the experimental observations, is favored for small $U$ values, and becomes more and more favorable compared to the ferromagnetic state towards the end of the series (i.e., towards $R$=Pr). Nevertheless, it seems that the stability of the ferromagnetic state is generally overestimated within the DFT+$U$ calculations. Our work provides a profound starting point for more detailed experimental investigations, and also for future studies using more advanced computational techniques such as, e.g., DFT combined with dynamical mean-field theory.


title: Rapid Exponentiation using Discrete Operators: Applications in Optimizing Quantum Controls and Simulating Quantum Dynamics

abstract: Matrix exponentiation (ME) is widely used in various fields of science and engineering. For example, the unitary dynamics of quantum systems is described by exponentiation of Hamiltonian operators. However, despite a significant attention, the numerical evaluation of ME remains computationally expensive, particularly for large dimensions. Often this process becomes a bottleneck in algorithms requiring iterative evaluation of ME. Here we propose a method for approximating ME of a single operator with a bounded coefficient into a product of certain discrete operators. This approach, which we refer to as Rapid Exponentiation using Discrete Operators (REDO), is particularly efficient for iterating ME over large numbers. We describe REDO in the context of a quantum system with a constant as well as a time-dependent Hamiltonian, although in principle, it can be adapted in a more general setting. As concrete examples, we choose two applications. First, we incorporate REDO in optimal quantum control algorithms and report a speed-up of several folds over a wide range of system size. Secondly, we propose REDO for numerical simulations of quantum dynamics. In particular, we study exotic quantum freezing with noisy drive fields.


title: Circular-shift Linear Network Coding

abstract: We study a class of linear network coding (LNC) schemes, called circular-shift LNC, whose encoding operations consist of only circular-shifts and bit-wise additions (XOR). Formulated as a special vector linear code over GF($2$), an $L$-dimensional circular-shift linear code of degree $\delta$ restricts its local encoding kernels to be the summation of at most $\delta$ cyclic permutation matrices of size $L$. We show that on a general network, for a certain block length $L$, every scalar linear solution over GF($2^{L-1}$) can induce an $L$-dimensional circular-shift linear solution with 1-bit redundancy per-edge transmission. Consequently, specific to a multicast network, such a circular-shift linear solution of an arbitrary degree $\delta$ can be efficiently constructed, which has an interesting complexity tradeoff between encoding and decoding with different choices of $\delta$. By further proving that circular-shift LNC is insufficient to achieve the exact capacity of certain multicast networks, we show the optimality of the efficiently constructed circular-shift linear solution in the sense that its 1-bit redundancy is inevitable. Finally, both theoretical and numerical analysis imply that with increasing $L$, a randomly constructed circular-shift linear code has linear solvability behavior comparable to a randomly constructed permutation-based linear code, but has shorter overheads.


title: $\xi/\xi_{2nd}$ ratio as a tool to refine Effective Polyakov Loop models

abstract: Effective Polyakov line actions are a powerful tool to study the finite temperature behaviour of lattice gauge theories. They are much simpler to simulate than the original lattice model and are affected by a milder sign problem, but it is not clear to which extent they really capture the rich spectrum of the original theories. We propose here a simple way to address this issue based on the so called second moment correlation length $\xi_{2nd}$. The ratio $\xi/\xi_{2nd}$ between the exponential correlation length and the second moment one is equal to 1 if only a single mass is present in the spectrum, and it becomes larger and larger as the complexity of the spectrum increases. Since both $\xi$ and $\xi_{2nd}$ are easy to measure on the lattice, this is a cheap and efficient way to keep track of the spectrum of the theory. As an example of the information one can obtain with this tool we study the behaviour of $\xi/\xi_{2nd}$ in the confining phase of the ($D=3+1$) $\mathrm{SU}(2)$ gauge theory and show that it is compatible with 1 near the deconfinement transition, but it increases dramatically as the temperature decreases. We also show that this increase can be well understood in the framework of an effective string description of the Polyakov loop correlator. This non-trivial behaviour should be reproduced by the Polyakov loop effective action; thus, it represents a stringent and challenging test of existing proposals and it may be used to fine-tune the couplings and to identify the range of validity of the approximations involved in their construction.


title: Stellar Coronal Mass Ejections I. Estimating occurrence frequencies and mass-loss rates

abstract: Stellar coronal mass ejections (CMEs) may play an important role in mass- and angular momentum loss of young Sun-like stars. If occurring frequently, they may also have a strong effect on planetary evolution by increasing atmospheric erosion. So far it has not been possible to infer the occurrence frequency of stellar CMEs from observations. Based on their close relation with flares on the Sun, we develop an empirical model combining solar flare-CME relationships with stellar flare rates to estimate the CME activity of young Sun-like and late-type main-sequence stars. By comparison of the obtained CME mass-loss rates with observations of total mass-loss rates, we find that our modeled rates may exceed those from observations by orders of magnitude for the most active stars. This reveals a possible limit to the extrapolation of such models to the youngest stars. We find that the most uncertain component in the model is the flare-CME association rate adopted from the Sun, which does not properly account for the likely stronger coronal confinement in active stars. Simple estimates of this effect reveal a possible suppression of CME rates by several orders of magnitude for young stars, indicating that this issue should be addressed in more detail in the future.


title: A new continuum theory for incompressible swelling materials

abstract: Swelling media (e.g. gels, tumors) are usually described by mechanical constitutive laws (e.g. Hooke or Darcy laws). However, constitutive relations of real swelling media are not well known. Here, we take an opposite route and consider a simple packing heuristics, i.e. the particles can't overlap. We deduce a formula for the equilibrium density under a confining potential. We then consider its evolution when the average particle volume and confining potential depend on time under two additional heuristics: (i) any two particles can't swap their position; (ii) motion should obey some energy minimization principle. These heuristics determine the medium velocity consistently with the continuity equation. In the direction normal to the potential level sets the velocity is related with that of the level sets while in the parallel direction, it is determined by a Laplace-Beltrami operator on these sets. This complex geometrical feature cannot be recovered using a simple Darcy law.


title: Tensor-triangular fields: Ruminations

abstract: We examine the concept of field in tensor-triangular geometry. We gather examples and discuss possible approaches, while highlighting open problems. As the construction of residue tt-fields remains elusive, we instead produce suitable homological tensor-functors to Grothendieck categories.


title: Dynamic dipole polarizabilities of heteronuclear alkali dimers: optical response, trapping and control of ultracold molecules

abstract: In this article we address the general approach for calculating dynamical dipole polarizabilities of small quantum systems, based on a sum-over-states formula involving in principle the entire energy spectrum of the system. We complement this method by a few-parameter model involving a limited number of effective transitions, allowing for a compact and accurate representation of both the isotropic and anisotropic components of the polarizability. We apply the method to the series of ten heteronuclear molecules composed of two of ($^7$Li,$^{23}$Na,$^{39}$K,$^{87}$Rb,$^{133}$Cs) alkali-metal atoms. We rely on both up-to-date spectroscopically-determined potential energy curves for the lowest electronic states, and on our systematic studies of these systems performed during the last decade for higher excited states and for permanent and transition dipole moments. Such a compilation is timely for the continuously growing researches on ultracold polar molecules. Indeed the knowledge of the dynamic dipole polarizabilities is crucial to model the optical response of molecules when trapped in optical lattices, and to determine optimal lattice frequencies ensuring optimal transfer to the absolute ground state of initially weakly-bound molecules. When they exist, we determine the so-called "magic frequencies" where the ac-Stark shift and thus the viewed trap depth, is the same for both weakly-bound and ground-state molecules.


title: Enhanced photothermal cooling of nanowires

abstract: We investigate the optomechanical interaction between light and metallic nanowires through the action of bolometric forces. We show that the response time of the photothermal forces induced on the nanowire is fast and the strength of the interaction can overcome the radiation pressure force. Furthermore, we suggest the photothermal forces can be enhanced by surface plasmon excitation to cool the sub-megahertz vibrational modes of the nanowires close to its quantum limit.


title: Barnett relaxation in non-symmetric grains

abstract: Barnett relaxation, first described by Purcell in 1979, appears to play a major role in the alignment of grains with the interstellar magnetic field. In 1999, Lazarian and Draine proposed that Barnett relaxation and its relative, nuclear relaxation, can induce grains to flip. If this thermal flipping is rapid, then the dynamical effect of torques that are fixed relative to the grain body can be greatly reduced. To date, detailed studies of Barnett relaxation have been confined to grains exhibiting dynamic symmetry. In 2003, Weingartner argued that internal relaxation cannot induce flips in any grains, whether they exhibit dynamic symmetry or not. In this work, we develop approximate expressions for the dissipation rate and diffusion coefficient for Barnett relaxation. We revisit the issue of internally induced thermal flipping, finding that it cannot occur for grains with dynamic symmetry but does occur for grains lacking dynamic symmetry.


title: Interpreting and using CPDAGs with background knowledge

abstract: We develop terminology and methods for working with maximally oriented partially directed acyclic graphs (maximal PDAGs). Maximal PDAGs arise from imposing restrictions on a Markov equivalence class of directed acyclic graphs, or equivalently on its graphical representation as a completed partially directed acyclic graph (CPDAG), for example when adding background knowledge about certain edge orientations. Although maximal PDAGs often arise in practice, causal methods have been mostly developed for CPDAGs. In this paper, we extend such methodology to maximal PDAGs. In particular, we develop methodology to read off possible ancestral relationships, we introduce a graphical criterion for covariate adjustment to estimate total causal effects, and we adapt the IDA and joint-IDA frameworks to estimate multi-sets of possible causal effects. We also present a simulation study that illustrates the gain in identifiability of total causal effects as the background knowledge increases. All methods are implemented in the R package pcalg.


title: Statics and kinematics of frameworks in Euclidean and non-Euclidean geometry

abstract: This is a survey article on the infinitesimal rigidity of frameworks in Euclidean, hyperbolic, and spherical geometry. We discuss the equivalence of the static and kinematic formulations of the infinitesimal rigidity, the projective interpretation of statics (representing forces as bivectors), and the infinitesimal Pogorelov maps that establish correspondence between infinitesimal motions of a framework and of its geodesic image. Also we describe the Maxwell-Cremona correspondence between equilibrium loads and polyhedral lifts, both for Euclidean and for non-Euclidean frameworks.


title: Self-organized criticality in type I X-ray bursts

abstract: Type I X-ray bursts in a low-mass X-ray binary (LMXB) are caused by unstable nuclear burning of accreted materials. Semi-analytical and numerical studies of unstable nuclear burning have successfully reproduced partial properties of this kind of burst. However, some other properties (e.g. the waiting time) are not well explained. In this paper, we find that the probability distributions of fluence, peak count, rise time, duration and waiting time can be described as power-law-like distributions. This indicates that type I X-ray bursts may be governed by a self-organized criticality (SOC) process. The power-law index of waiting time distribution (WTD) is around $-1$, which is not predicted by any current waiting time model. We propose a physical burst rate model, in which the mean occurrence rate is inversely proportional to time $\lambda\propto t^{-1}$. In this case, the WTD is well explained by a non-stationary Poisson process within the SOC theory. In this theory, the burst size is also predicted to follow a power-law distribution, which requires that the emission area possesses only part of the neutron star surface. Furthermore, we find that the WTDs of some astrophysical phenomena can also be described by similar occurrence rate models.


title: Methods for finding leader--follower equilibria with multiple followers

abstract: The concept of leader--follower (or Stackelberg) equilibrium plays a central role in a number of real--world applications of game theory. While the case with a single follower has been thoroughly investigated, results with multiple followers are only sporadic and the problem of designing and evaluating computationally tractable equilibrium-finding algorithms is still largely open. In this work, we focus on the fundamental case where multiple followers play a Nash equilibrium once the leader has committed to a strategy---as we illustrate, the corresponding equilibrium finding problem can be easily shown to be $\mathcal{FNP}$--hard and not in Poly--$\mathcal{APX}$ unless $\mathcal{P} = \mathcal{NP}$ and therefore it is one among the hardest problems to solve and approximate. We propose nonconvex mathematical programming formulations and global optimization methods to find both exact and approximate equilibria, as well as a heuristic black box algorithm. All the methods and formulations that we introduce are thoroughly evaluated computationally.


title: Exoplanets and SETI

abstract: The discovery of exoplanets has both focused and expanded the search for extraterrestrial intelligence. The consideration of Earth as an exoplanet, the knowledge of the orbital parameters of individual exoplanets, and our new understanding of the prevalence of exoplanets throughout the galaxy have all altered the search strategies of communication SETI efforts, by inspiring new "Schelling points" (i.e. optimal search strategies for beacons). Future efforts to characterize individual planets photometrically and spectroscopically, with imaging and via transit, will also allow for searches for a variety of technosignatures on their surfaces, in their atmospheres, and in orbit around them. In the near-term, searches for new planetary systems might even turn up free-floating megastructures.


title: Three particle quantization condition in a finite volume: 2. general formalism and the analysis of data

abstract: We derive the three-body quantization condition in a finite volume using an effective field theory in the particle-dimer picture. Moreover, we consider the extraction of physical observables from the lattice spectrum using the quantization condition. To illustrate the general framework, we calculate the volume-dependent three-particle spectrum in a simple model both below and above the three-particle threshold. The relation to existing approaches is discussed in detail.


title: Time-Dependent Ionization in theEnvelopes of Type II Supernovae at the Photospheric Phase

abstract: The importance of allowance for the time-dependent effect in the kinetics at the photospheric phase during a supernova explosion has been confirmed by several independent research groups. The time-dependent effect provides a higher degree of hydrogen ionization in comparison with the steady-state solutions and strengthens the H$\alpha$ line in the resulting simulated spectrum, with the intensity of the effect increasing with time. However, some researchers argue that the time-dependent ionization effect is unimportant. Its allowance leads to an insignificant strengthening of H$\alpha$ in their modeling only in the first days after explosion. We have demonstrated the importance of the time-dependent effect with the models of SN 1999em as an example using the new original LEVELS software package. The role of a number of factors that can weaken the time-dependent effect has been checked. We have confirmed that the intensity of the effect is affected by the abundance of metal admixtures in the envelope, while the addition of extra levels to themodel hydrogen atom weakens the time-dependent effect to a lesser degree and never removes it completely.


title: Applying Parabolic Peterson: Affine Algebras and the Quantum Cohomology of the Grassmannian

abstract: The Peterson isomorphism relates the homology of the affine Grassmannian to the quantum cohomology of any flag variety. In the case of a partial flag, Peterson's map is only a surjection, and one needs to quotient by a suitable ideal on the affine side to map isomorphically onto the quantum cohomology. We provide a detailed exposition of this parabolic Peterson isomorphism in the case of the Grassmannian of m-planes in complex n-space, including an explicit recipe for doing quantum Schubert calculus in terms of the appropriate subset of non-commutative k-Schur functions. As an application, we recast Postnikov's affine approach to the quantum cohomology of the Grassmannian as a consequence of parabolic Peterson by showing that the affine nilTemperley-Lieb algebra arises naturally when forming the requisite quotient of the homology of the affine Grassmannian.


title: Young's moduli of carbon materials investigated by various classical molecular dynamics schemes

abstract: Classical carbon potentials together with classical molecular dynamics are employed to calculate structures and physical properties of such carbon-based materials where quantum mechanical methods fail either due to the excessive size, irregular structure or long-time dynamics. Examples are given by recently synthesized free-standing carbon nanomembranes (CNM) with molecular thickness and macroscopic lateral size as well as by amorphous carbon. Although such potentials, as for instance implemented in LAMMPS, yield reasonably accurate bond lengths and angles for several carbon materials such as graphene, it is not clear how accurate they are in terms of mechanical properties such as Young's moduli. We performed large-scale classical molecular dynamics investigations of three carbon-based materials using the various potentials implemented in LAMMPS as well as the highly sophisticated EDIP potential of Nigel Marks. We demonstrate how the Young's moduli vary with classical potentials and compare to experimental results. Since classical descriptions of carbon are bound to be approximations it is not astonishing that different realizations yield differing results. One should therefore carefully check for which observables a certain potential is suited. We hope to contribute to such a clarification.


title: Signatures of rocky planet engulfment in HAT-P-4. Implications for chemical tagging studies

abstract: Aims. To explore the possible chemical signature of planet formation in the binary system HAT-P-4, by studying abundance vs condensation temperature Tc trends. The star HAT-P-4 hosts a planet detected by transits while its stellar companion does not have any detected planet. We also study the Lithium content, which could shed light on the problem of Li depletion in exoplanet host stars. Conclusions. The exoplanet host star HAT-P-4 is found to be ~0.1 dex more metal rich than its companion, which is one of the highest differences in metallicity observed in similar systems. This could have important implications for chemical tagging studies, disentangling groups of stars with a common origin. We rule out a possible peculiar composition for each star as lambda Boo, delta Scuti or a Blue Straggler. The star HAT-P-4 is enhanced in refractory elements relative to volatile when compared to its stellar companion. Notably, the Lithium abundance in HAT-P-4 is greater than in its companion by ~0.3 dex, which is contrary to the model that explains the Lithium depletion by the presence of planets. We propose a scenario where, at the time of planet formation, the star HAT-P-4 locked the inner refractory material in planetesimals and rocky planets, and formed the outer gas giant planet at a greater distance. The refractories were then accreted onto the star, possibly due to the migration of the giant planet. This explains the higher metallicity, the higher Lithium content, and the negative Tc trend detected. A similar scenario was recently proposed for the solar twin star HIP 68468, which is in some aspects similar to HAT-P-4. We estimate a mass of at least Mrock ~ 10 Mearth locked in refractory material in order to reproduce the observed Tc trends and metallicity.


title: Real eigenvalues in the non-Hermitian Anderson model

abstract: The eigenvalues of the Hatano--Nelson non-Hermitian Anderson matrices, in the spectral regions in which the Lyapunov exponent exceeds the non-Hermiticity parameter, are shown to be real and exponentially close to the Hermitian eigenvalues. This complements previous results, according to which the eigenvalues in the spectral regions in which the non-Hermiticity parameter exceeds the Lyapunov exponent are aligned on curves in the complex plane.


title: A bi-dimensional finite mixture model for longitudinal data subject to dropout

abstract: In longitudinal studies, subjects may be lost to follow-up, or miss some of the planned visits, leading to incomplete response sequences. When the probability of non-response, conditional on the available covariates and the observed responses, still depends on unobserved outcomes, the dropout mechanism is said to be non ignorable. A common objective is to build a reliable association structure to account for dependence between the longitudinal and the dropout processes. Starting from the existing literature, we introduce a random coefficient based dropout model where the association between outcomes is modeled through discrete latent effects. These effects are outcome-specific and account for heterogeneity in the univariate profiles. Dependence between profiles is introduced by using a bi-dimensional representation for the corresponding distribution. In this way, we define a flexible latent class structure which allows to efficiently describe both dependence within the two margins of interest and dependence between them. By using this representation we show that, unlike standard (unidimensional) finite mixture models, the non ignorable dropout model properly nests its ignorable counterpart. We detail the proposed modeling approach by analyzing data from a longitudinal study on the dynamics of cognitive functioning in the elderly. Further, the effects of assumptions about non ignorability of the dropout process on model parameter estimates are (locally) investigated using the index of (local) sensitivity to non-ignorability.


title: Browkin's discriminator conjecture

abstract: Let $q\ge 5$ be a prime and put $q^*=(-1)^{(q-1)/2}\cdot q$. We consider the integer sequence $u_q(1),u_q(2),\ldots,$ with $u_q(j)=(3^j-q^*(-1)^j)/4$. No term in this sequence is repeated and thus for each $n$ there is a smallest integer $m$ such that $u_q(1),\ldots,u_q(n)$ are pairwise incongruent modulo $m$. We write $D_q(n)=m$. The idea of considering the discriminator $D_q(n)$ is due to Browkin (2015) who, in case $3$ is a primitive root modulo $q,$ conjectured that the only values assumed by $D_q(n)$ are powers of $2$ and of $q$. We show that this is true for $n\neq 5$, but false for infinitely many $q$ in case $n=5$. We also determine $D_q(n)$ in case 3 is not a primitive root modulo $q$. Browkin's inspiration for his conjecture came from earlier work of Moree and Zumalac\'arregui (2016), who determined $D_5(n)$ for $n\ge 1$, thus establishing a conjecture of Salajan. For a fixed prime $q$ their approach is easily generalized, but requires some innovations in order to deal with all primes $q\ge 7$ and all $n\ge 1$. Interestingly enough, Fermat and Mirimanoff primes play a special role in this.


title: Measurement of the $\pi K$ atom lifetime and the $\pi K$ scattering length

abstract: After having announced the statistically significant observation (5.6~$\sigma$) of the new exotic $\pi K$ atom, the DIRAC experiment at the CERN proton synchrotron presents the measurement of the corresponding atom lifetime, based on the full $\pi K$ data sample: $\tau = (5.5^{+5.0}_{-2.8}) \cdot 10^{-15}s$. By means of a precise relation ($<1\%$) between atom lifetime and scattering length, the following value for the S-wave isospin-odd $\pi K$ scattering length $a_0^{-}~=~\frac{1}{3}(a_{1/2}-a_{3/2})$ has been derived: $\left|a_0^-\right| = (0.072^{+0.031}_{-0.020}) M_{\pi}^{-1}$.


title: Threshold driven contagion on weighted networks

abstract: Weighted networks capture the structure of complex systems where interaction strength is meaningful. This information is essential to a large number of processes, such as threshold dynamics, where link weights reflect the amount of influence that neighbours have in determining a node's behaviour. Despite describing numerous cascading phenomena, such as neural firing or social contagion, threshold models have never been explicitly addressed on weighted networks. We fill this gap by studying a dynamical threshold model over synthetic and real weighted networks with numerical and analytical tools. We show that the time of cascade emergence depends non-monotonously on weight heterogeneities, which accelerate or decelerate the dynamics, and lead to non-trivial parameter spaces for various networks and weight distributions. Our methodology applies to arbitrary binary state processes and link properties, and may prove instrumental in understanding the role of edge heterogeneities in various natural and social phenomena.

categories: physics.soc-ph cond-mat.stat-mech cs.SI physics.data-an


title: Comptage probabiliste sur la fronti\`ere de Furstenberg

abstract: Let $G$ be a real linear semisimple algebraic group without compact factors and $\Gamma$ a Zariski dense subgroup of $G$. In this paper, we use a probabilistic counting in order to study the asymptotic properties of $\Gamma$ acting on the Furstenberg boundary of $G$. First, we show that the $K$ components of the elements of $\Gamma$ in the KAK decomposition of $G$ become asymptotically independent. This result is an analog of a result of Gorodnik-Oh in the context of the Archimedean counting. Then, we give a new proof of a result of Guivarc'h concerning the positivity of the Hausdorff dimension of the unique stationary probability measure on the Furstenberg Boundary of $G$. Finally, we show how these results can be combined to give a probabilistic proof of the Tit's alternative; namely that two independent random walks on $\Gamma$ will eventually generate a free subgroup. This result answered a question of Guivarc'h and was published earlier by the author. Since we're working with the field of real numbers, we give here a more direct proof and a more general statement.


title: Binary strings of length $n$ with $x$ zeros and longest $k$-runs of zeros

abstract: In this paper, we study $F_{n}(x,k)$, the number of binary strings of length $n$ containing $x$ zeros and a longest subword of $k$ zeros. A recurrence relation for $F_{n}(x,k)$ is derived. We expressed few known numbers like Fibonacci, triangular, number of binary strings of length $n$ without $r$-runs of ones and number of compositions of $n+1$ with largest summand $k+1$ in terms of $F_{n}(x,k).$ Similar results and applications were obtained for \^F$_{n}(x,k),$ the number of all palindromic binary strings of length $n$ containing $x$ zeros and longest $k$-runs of zeros.


title: Coherent diversification in corporate technological portfolios

abstract: We study the relationship between firms' performance and their technological portfolios using tools borrowed from the complexity science. In particular, we ask whether the accumulation of knowledge and capabilities related to a coherent set of technologies leads firms to experience advantages in terms of productive efficiency. To this end, we analyzed both the balance sheets and the patenting activity of about 70 thousand firms that have filed at least one patent over the period 2004-2013. From this database it is possible to define a measure of the firms' coherent diversification, based on the network of technological fields, and relate it to the firms' perfomance in terms of labor productivity. Such a measure favors companies with a diversification structure comprising blocks of closely related fields over firms with the same breadth of scope, but a more scattered diversification structure. We find that the coherent diversification of firms is quantitatively related to their economic performance and captures relevant information about their productive structure. In particular, we prove on a statistical basis that a naive definition of technological diversification can explain labor productivity only as a proxy of size and coherent diversification. This approach can be used to investigate possible synergies within firms and to recommend viable partners for merging and acquisitions.


title: Characterization of Chromium Compensated GaAs as an x-ray Sensor Material for Charge-Integrating Pixel Array Detectors

abstract: We studied the properties of chromium compensated GaAs when coupled to charge integrating ASICs as a function of detector temperature, applied bias and x-ray tube energy. The material is a photoresistor and can be biased to collect either electrons or holes by the pixel circuitry. Both are studied here. Previous studies have shown substantial hole trapping. This trapping and other sensor properties give rise to several non-ideal effects which include an extended point spread function, variations in the effective pixel size, and rate dependent offset shifts. The magnitude of these effects varies with temperature and bias, mandating good temperature uniformity in the sensor and very good temperature stabilization, as well as a carefully selected bias voltage.


title: A subexponential parameterized algorithm for Directed Subset Traveling Salesman Problem on planar graphs

abstract: There are numerous examples of the so-called ``square root phenomenon'' in the field of parameterized algorithms: many of the most fundamental graph problems, parameterized by some natural parameter $k$, become significantly simpler when restricted to planar graphs and in particular the best possible running time is exponential in $O(\sqrt{k})$ instead of $O(k)$ (modulo standard complexity assumptions). We consider a classic optimization problem Subset Traveling Salesman, where we are asked to visit all the terminals $T$ by a minimum-weight closed walk. We investigate the parameterized complexity of this problem in planar graphs, where the number $k=|T|$ of terminals is regarded as the parameter. We show that Subset TSP can be solved in time $2^{O(\sqrt{k}\log k)}\cdot n^{O(1)}$ even on edge-weighted directed planar graphs. This improves upon the algorithm of Klein and Marx [SODA 2014] with the same running time that worked only on undirected planar graphs with polynomially large integer weights.


title: Design and Processing of Invertible Orientation Scores of 3D Images for Enhancement of Complex Vasculature

abstract: The enhancement and detection of elongated structures in noisy image data is relevant for many biomedical imaging applications. To handle complex crossing structures in 2D images, 2D orientation scores $U: \mathbb{R} ^ 2\times S ^ 1 \rightarrow \mathbb{C}$ were introduced, which already showed their use in a variety of applications. Here we extend this work to 3D orientation scores $U: \mathbb{R} ^ 3 \times S ^ 2\rightarrow \mathbb{C}$. First, we construct the orientation score from a given dataset, which is achieved by an invertible coherent state type of transform. For this transformation we introduce 3D versions of the 2D cake-wavelets, which are complex wavelets that can simultaneously detect oriented structures and oriented edges. Here we introduce two types of cake-wavelets, the first uses a discrete Fourier transform, the second is designed in the 3D generalized Zernike basis, allowing us to calculate analytical expressions for the spatial filters. Finally, we show two applications of the orientation score transformation. In the first application we propose an extension of crossing-preserving coherence enhancing diffusion via our invertible orientation scores of 3D images which we apply to real medical image data. In the second one we develop a new tubularity measure using 3D orientation scores and apply the tubularity measure to both artificial and real medical data.


title: Ambitwistor formulations of $R^2$ gravity and $(DF)^2$ gauge theories

abstract: We consider $D$-dimensional amplitudes in $R^2$ gravities (conformal gravity in $D=4$) and in the recently introduced $(DF)^2$ gauge theory, from the perspective of the CHY formulae and ambitwistor string theory. These theories are related through the BCJ double-copy construction, and the $(DF)^2$ gauge theory obeys color-kinematics duality. We work out the worldsheet details of these theories and show that they admit a formulation as integrals on the support of the scattering equations, or alternatively, as ambitwistor string theories. For gravity, this generalizes the work done by Berkovits and Witten on conformal gravity to $D$ dimensions. The ambitwistor is also interpreted as a $D$-dimensional generalization of Witten's twistor string (SYM + conformal supergravity). As part of our ambitwistor investigation, we discover another $(DF)^2$ gauge theory containing a photon that couples to Einstein gravity. This theory can provide an alternative KLT description of Einstein gravity compared to the usual Yang-Mills squared.


title: The helium Warm Breeze in IBEX observations as a result of charge exchange collisions in the outer heliosheath

abstract: We simulated the signal due to neutral He atoms, observed by Interstellar Boundary Explorer (IBEX), assuming that charge exchange collisions between neutral He atoms and He+ ions operate everywhere between the heliopause and a distant source region in the local interstellar cloud, where the neutral and charged components are in thermal equilibrium. We simulated several test cases of the plasma flow within the outer heliosheath and investigated the signal generation for plasma flows both in the absence and in the presence of the interstellar magnetic field. We found that a signal in the portion of IBEX data identified as due to the Warm Breeze does not arise when a homogeneous plasma flow in front of the heliopause is assumed, but it appears immediately when any reasonable disturbance in its flow due to the presence of the heliosphere is assumed. We obtained a good qualitative agreement between the data selected for comparison and the simulations for a model flow with the velocity vector of the unperturbed gas and the direction and intensity of magnetic field adopted from recent determinations. We conclude that direct-sampling observations of neutral He atoms at 1~AU from the Sun are a sensitive tool for investigating the flow of interstellar matter in the outer heliosheath, that the Warm Breeze is indeed the secondary population of interstellar helium, as it was hypothesized earlier, and that the WB signal is consistent with the heliosphere distorted from axial symmetry by the interstellar magnetic field.


title: A multi-layer image representation using Regularized Residual Quantization: application to compression and denoising

abstract: A learning-based framework for representation of domain-specific images is proposed where joint compression and denoising can be done using a VQ-based multi-layer network. While it learns to compress the images from a training set, the compression performance is very well generalized on images from a test set. Moreover, when fed with noisy versions of the test set, since it has priors from clean images, the network also efficiently denoises the test images during the reconstruction. The proposed framework is a regularized version of the Residual Quantization (RQ) where at each stage, the quantization error from the previous stage is further quantized. Instead of codebook learning from the k-means which over-trains for high-dimensional vectors, we show that only generating the codewords from a random, but properly regularized distribution suffices to compress the images globally and without the need to resort to patch-based division of images. The experiments are done on the \textit{CroppedYale-B} set of facial images and the method is compared with the JPEG-2000 codec for compression and BM3D for denoising, showing promising results.


title: Quantum interface between photonic and superconducting qubits

abstract: We show that optically active coupled quantum dots embedded in a superconducting microwave cavity can be used to realize a fast quantum interface between photonic and transmon qubits. Single photon absorption by a coupled quantum dot results in generation of a large electric dipole, which in turn ensures efficient coupling to the microwave cavity. Using cavity parameters achieved in prior experiments, we estimate that bi-directional microwave-optics conversion in nanosecond timescales with efficiencies approaching unity is experimentally feasible with current technology. We also outline a protocol for in-principle deterministic quantum state transfer from a time-bin photonic qubit to a transmon qubit. Recent advances in quantum dot based quantum photonics technologies indicate that the scheme we propose could play a central role in connecting quantum nodes incorporating cavity-coupled superconducting qubits.


title: Infinite-server queues with Hawkes input

abstract: In this paper we study the number of customers in infinite-server queues with a self-exciting (Hawkes) arrival process. Initially we assume that service requirements are exponentially distributed and that the Hawkes arrival process is of a Markovian nature. We obtain a system of differential equations that characterizes the joint distribution of the arrival intensity and the number of customers. Moreover, we provide a recursive procedure that explicitly identifies (transient and stationary) moments. Subsequently, we allow for non-Markovian Hawkes arrival processes and non-exponential service times. By viewing the Hawkes process as a branching process, we find that the probability generating function of the number of customers in the system can be expressed in terms of the solution of a fixed-point equation. We also include various asymptotic results: we derive the tail of the distribution of the number of customers for the case that the intensity jumps of the Hawkes process are heavy-tailed, and we consider a heavy-traffic regime. We conclude the paper by discussing how our results can be used computationally and by verifying the numerical results via simulations.


title: Chaotic-Integrable Transition in the Sachdev-Ye-Kitaev Model

abstract: Quantum chaos is one of the distinctive features of the Sachdev-Ye-Kitaev (SYK) model, $N$ Majorana fermions in $0+1$ dimensions with infinite-range two-body interactions, which is attracting a lot of interest as a toy model for holography. Here we show analytically and numerically that a generalized SYK model with an additional one-body infinite-range random interaction, which is a relevant perturbation in the infrared, is still quantum chaotic and retains most of its holographic features for a fixed value of the perturbation and sufficiently high temperature. However a chaotic-integrable transition, characterized by the vanishing of the Lyapunov exponent and spectral correlations given by Poisson statistics, occurs at a temperature that depends on the strength of the perturbation. We speculate about the gravity dual of this transition.


title: Learning Loss Functions for Semi-supervised Learning via Discriminative Adversarial Networks

abstract: We propose discriminative adversarial networks (DAN) for semi-supervised learning and loss function learning. Our DAN approach builds upon generative adversarial networks (GANs) and conditional GANs but includes the key differentiator of using two discriminators instead of a generator and a discriminator. DAN can be seen as a framework to learn loss functions for predictors that also implements semi-supervised learning in a straightforward manner. We propose instantiations of DAN for two different prediction tasks: classification and ranking. Our experimental results on three datasets of different tasks demonstrate that DAN is a promising framework for both semi-supervised learning and learning loss functions for predictors. For all tasks, the semi-supervised capability of DAN can significantly boost the predictor performance for small labeled sets with minor architecture changes across tasks. Moreover, the loss functions automatically learned by DANs are very competitive and usually outperform the standard pairwise and negative log-likelihood loss functions for both semi-supervised and supervised learning.


title: Gr\"obner Bases for Schubert Codes

abstract: We consider the problem of determining Gr\"obner bases of binomial ideals associated with linear error correcting codes. Computation of Gr\"obner bases of linear codes have become a topic of interest to many researchers in coding theory because of its several applications in decoding and error corrections. In this paper, Gr\"obner bases of linear codes associated to Grassmann varieties and Schubert varieties over a binary field have been obtained. We also use them to study the decoding of binary Schubert codes.


title: Few Non-derogatory Directed Graphs from Directed Cycles

abstract: We constructed a few non-derogatory digraphs by adding arcs to a directed cycle and computed their characteristic polynomials and exponents.


title: Learning human behaviors from motion capture by adversarial imitation

abstract: Rapid progress in deep reinforcement learning has made it increasingly feasible to train controllers for high-dimensional humanoid bodies. However, methods that use pure reinforcement learning with simple reward functions tend to produce non-humanlike and overly stereotyped movement behaviors. In this work, we extend generative adversarial imitation learning to enable training of generic neural network policies to produce humanlike movement patterns from limited demonstrations consisting only of partially observed state features, without access to actions, even when the demonstrations come from a body with different and unknown physical parameters. We leverage this approach to build sub-skill policies from motion capture data and show that they can be reused to solve tasks when controlled by a higher level controller.


title: Pandeia: A Multi-mission Exposure Time Calculator for JWST and WFIRST

abstract: Pandeia is the exposure time calculator (ETC) system developed for the James Webb Space Telescope (JWST) that will be used for creating JWST proposals. It includes a simulation-hybrid Python engine that calculates the two-dimensional pixel-by-pixel signal and noise properties of the JWST instruments. This allows for appropriate handling of realistic point spread functions, MULTIACCUM detector readouts, correlated detector readnoise, and multiple photometric and spectral extraction strategies. Pandeia includes support for all the JWST observing modes, including imaging, slitted/slitless spectroscopy, integral field spectroscopy, and coronagraphy. Its highly modular, data-driven design makes it easily adaptable to other observatories. An implementation for use with WFIRST is also available.


title: Non-adiabatic quantum state preparation and quantum state transport in chains of Rydberg atoms

abstract: Motivated by recent progress in the experimental manipulation of cold atoms in optical lattices, we study three different protocols for non-adiabatic quantum state preparation and state transport in chains of Rydberg atoms. The protocols we discuss are based on the blockade mechanism between atoms which, when excited to a Rydberg state, interact through a van der Waals potential, and rely on single-site addressing. Specifically, we discuss protocols for efficient creation of an antiferromagnetic GHZ state, a class of matrix product states including a so-called Rydberg crystal and for the state transport of a single-qubit quantum state between two ends of a chain of atoms. We identify system parameters allowing for the operation of the protocols on timescales shorter than the lifetime of the Rydberg states while yielding high fidelity output states. We discuss the effect of positional disorder on the resulting states and comment on limitations due to other sources of noise such as radiative decay of the Rydberg states. The proposed protocols provide a testbed for benchmarking the performance of quantum information processing platforms based on Rydberg atoms.


title: The structural role of the core literature in history

abstract: The intellectual landscapes of the humanities are mostly uncharted territory. Little is known on the ways published research of humanist scholars defines areas of intellectual activity. An open question relates to the structural role of core literature: highly cited sources, naturally playing a disproportionate role in the definition of intellectual landscapes. We introduce four indicators in order to map the structural role played by core sources into connecting different areas of the intellectual landscape of citing publications (i.e. communities in the bibliographic coupling network). All indicators factor out the influence of degree distributions by internalizing a null configuration model. By considering several datasets focused on history, we show that two distinct structural actions are performed by the core literature: a global one, by connecting otherwise separated communities in the landscape, or a local one, by rising connectivity within communities. In our study, the global action is mainly performed by small sets of scholarly monographs, reference works and primary sources, while the rest of the core, and especially most journal articles, acts mostly locally.


title: A proof of the Flaherty-Keller formula on the effective property of densely packed elastic composites

abstract: We prove in a mathematically rigorous way the asymptotic formula of Flaherty and Keller on the effective property of densely packed periodic elastic composites with hard inclusions. The proof is based on the primal-dual variational principle, where the upper bound is derived by using the Keller-type test functions and the lower bound by singular functions made of nuclei of strain. Singular functions are solutions of the Lam\'e system and capture precisely singular behavior of the stress in the narrow region between two adjacent hard inclusions.


title: Colloidal Gels Tuned by Oscillatory Shear

abstract: We examine microstructural and mechanical changes which occur during oscillatory shear flow and reformation after flow cessation of an intermediate volume fraction colloidal gel using rheometry and Brownian Dynamics (BD) simulations. A model depletion colloid-polymer mixture is used, comprising of a hard sphere colloidal suspension with the addition of non-adsorbing linear polymer chains. Results reveal three distinct regimes depending on the strain amplitude of oscillatory shear. Large shear strain amplitudes fully break the structure which results into a more homogenous and stronger gel after flow cessation. Intermediate strain amplitudes densify the clusters and lead to highly heterogeneous and weak gels. Shearing the gel to even lower strain amplitudes creates a less heterogonous stronger solid. These three regimes of shearing are connected to the microscopic shear-induced structural heterogeneity. A comparison with steady shear flow reveals that the latter does not produce structural heterogeneities as large as oscillatory shear. Therefore oscillatory shear is a much more efficient way of tuning the mechanical properties of colloidal gels. Moreover, colloidal gels presheared at large strain amplitudes exhibit a distinct nonlinear response characterized largely by a single yielding process while in those presheared at lower rates a two step yield process is promoted due to the creation of highly heterogeneous structures.


title: Quantitative characterization of stress concentration in the presence of closely spaced hard inclusions in two-dimensional linear elasticity

abstract: In the region between close-to-touching hard inclusions, the stress may be arbitrarily large as the inclusions get closer. The stress is represented by the gradient of a solution to the Lam\'e system of linear elasticity. We consider the problem of characterizing the gradient blow-up of the solution in the narrow region between two inclusions and estimating its magnitude. We introduce singular functions which are constructed in terms of nuclei of strain and hence are solutions of the Lam\'{e} system, and then show that the singular behavior of the gradient in the narrow region can be precisely captured by singular functions. As a consequence of the characterization, we are able to regain the existing upper bound on the blow-up rate of the gradient, namely, $\epsilon^{-1/2}$ where $\epsilon$ is the distance between two inclusions. We then show that it is in fact an optimal bound by showing that there are cases where $\epsilon^{-1/2}$ is also a lower bound. This work is the first to completely reveal the singular nature of the gradient blow-up in the context of the Lam\'{e} system with hard inclusions. The singular functions introduced in this paper play essential roles to overcome the difficulties in applying the methods of previous works. Main tools of this paper are the layer potential techniques and the variational principle. The variational principle can be applied because the singular functions of this paper are solutions of the Lam\'{e} system.


title: The properties of bordered matrix of symmetric block design

abstract: An incidence structure consists simply of a set P of points and a set B of blocks, with a relation of incidence between points and blocks.A symmetric (v,k,\lambda) block design is the subject of this paper. The symmetric (n^2+n+1, n+1,1) block design is a projective plane of order n. Despite much research no one has uncovered any further necessary conditions for the existence of a symmetric (v,k,\lambda) design apart from the equation (v-1)\lambda=k(k-1), Schutzenberger's Theorem and the Bruck-Ryser-Chowla Theorem. For no (v,k,\lambda) satisfying these requirements has it been shown that a symmetric (v,k,\lambda) design does not exist. Projective planes of order n exist for all prime powers n (aside from PG(2,n) a host of other constructions are known ) but for no other n is a construction known. The first open values are n=10, 12, 15, 18, 20, 24, 26$ and 28. It was proved by a computer search that there does not exist any projective plane of order 10 by Lam, C.W.H., Thiel, L. and Swiercz, S. Whether there exists any projective plane of order 12 is still open.The author introduces the bordered matrix of a (v,k,\lambda) symmetric design and gives some new necessary conditions for the existence of the symmetric (v,k,\lambda) design. As their application it is easy to determine that there does not exist finite projective plane of order n if n is one of the first open values 10, 12, 15, 18, 20, 24, 26 and 28, for which the Bruck-Ryser-Chowla Theorem can not be used. For large n the new method is also valid. Also some symmetric designs are excluded by the new method.


title: Analysis and Control of a Non-Standard Hyperbolic PDE Traffic Flow Model

abstract: The paper provides results for a non-standard, hyperbolic, 1-D, nonlinear traffic flow model on a bounded domain. The model consists of two first-order PDEs with a dynamic boundary condition that involves the time derivative of the velocity. The proposed model has features that are important from a traffic-theoretic point of view: is completely anisotropic and information travels forward exactly at the same speed as traffic. It is shown that, for all physically meaningful initial conditions, the model admits a globally defined, unique, classical solution that remains positive and bounded for all times. Moreover, it is shown that global stabilization can be achieved for arbitrary equilibria by means of an explicit boundary feedback law. The stabilizing feedback law depends only on the inlet velocity and consequently, the measurement requirements for the implementation of the proposed boundary feedback law are minimal. The efficiency of the proposed boundary feedback law is demonstrated by means of a numerical example.


title: On a Construction of Integrally Invertible Graphs and their Spectral Properties

abstract: Godsil (1985) defined a graph to be invertible if it has a non-singular adjacency matrix whose inverse is diagonally similar to a nonnegative integral matrix; the graph defined by the last matrix is then the inverse of the original graph. In this paper we call such graphs positively invertible and introduce a new concept of a negatively invertible graph by replacing the adjective `nonnegative' by `nonpositive in Godsil's definition; the graph defined by the negative of the resulting matrix is then the negative inverse of the original graph. We propose new constructions of integrally invertible graphs (those with non-singular adjacency matrix whose inverse is integral) based on an operation of `bridging' a pair of integrally invertible graphs over subsets of their vertices, with sufficient conditions for their positive and negative invertibility. We also analyze spectral properties of graphs arising from bridging and derive lower bounds for their least positive eigenvalue. As an illustration we present a census of graphs with a unique 1-factor on $m\le 6$ vertices and determine their positive and negative invertibility.


title: The Stellar decomposition: A compact representation for simplicial complexes and beyond

abstract: We introduce the Stellar decomposition, a model for efficient topological data structures over a broad range of simplicial and cell complexes. A Stellar decomposition of a complex is a collection of regions indexing the complex's vertices and cells such that each region has sufficient information to locally reconstruct the star of its vertices, i.e., the cells incident in the region's vertices. Stellar decompositions are general in that they can compactly represent and efficiently traverse arbitrary complexes with a manifold or non-manifold domain. They are scalable to complexes in high dimension and of large size, and they enable users to easily construct tailored application-dependent data structures using a fraction of the memory required by a corresponding global topological data structure on the complex. As a concrete realization of this model for spatially embedded complexes, we introduce the Stellar tree, which combines a nested spatial tree with a simple tuning parameter to control the number of vertices in a region. Stellar trees exploit the complex's spatial locality by reordering vertex and cell indices according to the spatial decomposition and by compressing sequential ranges of indices. Stellar trees are competitive with state-ofthe-art topological data structures for manifold simplicial complexes and offer significant improvements for cell complexes and non-manifold simplicial complexes. We conclude with a highlevel description of several mesh processing and analysis applications that utilize Stellar trees to process large datasets.


title: Radio-flaring Ultracool Dwarf Population Synthesis

abstract: Over a dozen ultracool dwarfs (UCDs), low-mass objects of spectral types $\geq$M7, are known to be sources of radio flares. These typically several-minutes-long radio bursts can be up to 100\% circularly polarized and have high brightness temperatures, consistent with coherent emission via the electron cyclotron maser operating in $\sim$kG magnetic fields. Recently, the statistical properties of the bulk physical parameters that describe these UCDs have become adequately described to permit synthesis of the population of radio-flaring objects. For the first time, I construct a Monte Carlo simulator to model the population of these radio-flaring UCDs. This simulator is powered by Intel Secure Key (ISK)- a new processor technology that uses a local entropy source to improve random number generation that has heretofore been used to improve cryptography. The results from this simulator indicate that only $\sim$5% of radio-flaring UCDs within the local interstellar neighborhood ($<$25 pc away) have been discovered. I discuss a number of scenarios which may explain this radio-flaring fraction, and suggest that the observed behavior is likely a result of several factors. The performance of ISK as compared to other pseudorandom number generators is also evaluated, and its potential utility for other astrophysical codes briefly described.

categories: astro-ph.SR astro-ph.IM physics.comp-ph physics.data-an


title: Photonic Nambu-Goldstone bosons

abstract: We study numerically the spatial dynamics of light in periodic square lattices in the presence of a Kerr term, emphasizing the peculiarities stemming from the nonlinearity. We find that, under rather general circumstances, the phase pattern of the stable ground state depends on the character of the nonlinearity: the phase is spatially uniform if it is defocusing whereas in the focusing case, it presents a chess board pattern, with a difference of $\pi$ between neighboring sites. We show that the lowest lying perturbative excitations can be described as perturbations of the phase and that finite-sized structures can act as tunable metawaveguides for them. The tuning is made by varying the intensity of the light that, because of the nonlinearity, affects the dynamics of the phase fluctuations. We interpret the results using methods of condensed matter physics, based on an effective description of the optical system. This interpretation sheds new light on the phenomena, facilitating the understanding of individual systems and leading to a framework for relating different problems with the same symmetry. In this context, we show that the perturbative excitations of the phase are Nambu-Goldstone bosons of a spontaneously broken $U(1)$ symmetry.


title: Optical trapping of anti-hydrogen towards an atomic anti-clock

abstract: The Anti-Matter Factory at CERN is gearing up, commissioning of the Extra Low ENergy Antiprotons (ELENA) ring is ongoing and the first anti-protons are foreseen to circulate in the decelerator very soon. The unprecedented flux of low energy antiprotons delivered by ELENA will open a new era for precision tests with antimatter including laser and microwave spectroscopy and tests of its gravitational behaviour. Here we propose a scheme to load the ultra cold anti-hydrogen atoms that will be produced by the GBAR experiment in an optical lattice tuned at the magic wavelength of the 1S-2S transition in order to measure this interval at a level comparable or even better than its matter counter part. This will provide the most sensitive test of CPT symmetry parametrised in the framework of the Standard Model Extension.


title: Mendelian randomization with fine-mapped genetic data: choosing from large numbers of correlated instrumental variables

abstract: Mendelian randomization uses genetic variants to make causal inferences about the effect of a risk factor on an outcome. With fine-mapped genetic data, there may be hundreds of genetic variants in a single gene region any of which could be used to assess this causal relationship. However, using too many genetic variants in the analysis can lead to spurious estimates and inflated Type 1 error rates. But if only a few genetic variants are used, then the majority of the data is ignored and estimates are highly sensitive to the particular choice of variants. We propose an approach based on summarized data only (genetic association and correlation estimates) that uses principal components analysis to form instruments. This approach has desirable theoretical properties: it takes the totality of data into account and does not suffer from numerical instabilities. It also has good properties in simulation studies: it is not particularly sensitive to varying the genetic variants included in the analysis or the genetic correlation matrix, and it does not have greatly inflated Type 1 error rates. Overall, the method gives estimates that are not so precise as those from variable selection approaches (such as using a conditional analysis or pruning approach to select variants), but are more robust to seemingly arbitrary choices in the variable selection step. Methods are illustrated by an example using genetic associations with testosterone for 320 genetic variants to assess the effect of sex hormone-related pathways on coronary artery disease risk, in which variable selection approaches give inconsistent inferences.


title: Generalized Nonsingular Solutions for the Scalar Meson Field of a Point Charge in General Relativity Theory

abstract: It has been shown that in investigating the electromagnetic and meson fields of elementary particles, the effects of gravitational interaction cannot be neglected. The present note is an attempt to find solutions to the equations of the general relativistic gravitational and scalar meson fields of a point nuclear charge. The solution is nonsingular at all points of the gravitational and meson fields. Although it contains several functions whose form is not given, this solution can be used to obtain the potential of the scalar meson field, which is a generalization of the Yukawa potential; furthermore, it makes it possible to calculate the mass and self-energy of the nucleon, which turn out to be finite.


title: General Covariant Equations for Fields of Arbitrary Spin

abstract: To obtain a generalized wave equation for a field in general covariant form, one must define covariant differentiation of a generalized wave function describing particles with arbitrary spin. Gel'fand and Iaglom, Dirac, and Fierz and Pauli have studied the generalized wave equation in the special theory of relativity. In this paper, the Gel'fand-Iaglom field equation was generalized to the curved space based on the veilbein framework. The general covariant field equations for arbitrary spin were obtained. More importantly, in this work, we introduced the gauge covariant derivative with the local Lorentz transformation, which is the foundation for a gauge theory of gravitation.


title: Universal magneto-orbital ordering in the divalent $A$-site quadruple perovskite manganites $A$Mn$_7$O$_{12}$ ($A$ = Ca, Sr, Cd, and Pb)

abstract: Through analysis of variable temperature neutron powder diffraction data, we present solutions for the magnetic structures of SrMn$_7$O$_{12}$, CdMn$_7$O$_{12}$, and PbMn$_7$O$_{12}$ in all long-range ordered phases. The three compounds were found to have magnetic structures analogous to that reported for CaMn$_7$O$_{12}$. They all feature a higher temperature lock-in phase with \emph{commensurate} magneto-orbital coupling, and a delocked, multi-\textbf{k} magnetic ground state where \emph{incommensurate} magneto-orbital coupling gives rise to a constant-moment magnetic helix with modulated spin helicity. CdMn$_7$O$_{12}$ represents a special case in which the orbital modulation is commensurate with the crystal lattice and involves stacking of fully and partially polarized orbital states. Our results provide a robust confirmation of the phenomenological model for magneto-orbital coupling previously presented for CaMn$_7$O$_{12}$. Furthermore, we show that the model is universal to the $A^{2+}$ quadruple perovskite manganites synthesised to date, and that it is tunable by selection of the $A$-site ionic radius.


title: On the growth behaviour of Hironaka quotients

abstract: We consider a finite analytic morphism $\phi = (f,g) : (X,p)\to (\C^2,0)$ where $(X,p)$ is a complex analytic normal surface germ and $f$ and $g$ are complex analytic function germs. Let $\pi : (Y,E_{Y})\to (X,p)$ be a good resolution of $\phi$ with exceptional divisor $E_{Y}=\pi ^{-1}(p)$. We denote $G(Y)$ the dual graph of the resolution $\pi $. We study the behaviour of the Hironaka quotients of $(f,g)$ associated to the vertices of $G(Y)$. We show that there exists maximal oriented arcs in $G(Y)$ along which the Hironaka quotients of $(f,g)$ strictly increase and they are constant on the connected components of the closure of the complement of the union of the maximal oriented arcs.


title: Scattered light intensity measurements of plasma treated Polydimethylsiloxane films: A measure to detect surface modification

abstract: Polydimethylsiloxane (PDMS) films possess different chemical and physical properties based on surface modification. The bond structure of pristine PDMS films and plasma treated PDMS films differ in a particular region of silicate bonds. We have studied the surface physical properties of pristine PDMS films and plasma treated PDMS films through optical technique. It is already known that plasma treated PDMS films forms very thin SiO2 layer on its surface. Due to difference in coefficient of thermal expansion of the surface SiO2 and the remaining bulk layer, the SiO2 layer develops cracks. These formations are explored to characterize PDMS surface by observing intensity of scattered light while the films are stretched. Pristine PDMS films do not show such optical scattering. Further the intensity measurements were repeated over period of time to monitor surface properties with time. It was observed that plasma treated PDMS film; the scattered light intensity is linearly dependent on the applied force for stretching. After a substantial time, scattering intensity is reduced to the value which is almost equal to that of pristine PDMS surface. It is proposed that basic light scattering through plasma treated PDMS occurs due to the decrease in the width of the cracks when the PDMS film is stretched. The change in the scattering property of plasma treated PDMS surface over time could be attributed to healing of cracks by the migration of polymer chain molecules from the bulk to the surface. It is also reported that with reference to various bonds present on surface plasma treated PDMS surface regains its original properties as that of pristine PDMS with time. Therefore the optical technique could be employed to study surface characteristics of PDMS surface as an alternative approach to conventional spectroscopic techniques.


title: Geant4 simulations of soft proton scattering in X-ray optics. A tentative validation using laboratory measurements

abstract: Low energy protons (< 300 keV) can enter the field of view of X-ray space telescopes, scatter at small incident angles, and deposit energy on the detector, causing intense background flares at the focal plane or in the most extreme cases, damaging the X-ray detector. A correct modelization of the physics process responsible for the grazing angle scattering processes is mandatory to evaluate the impact of such events on the performance of future X-ray telescopes as the ESA ATHENA mission. For the first time the Remizovich model, in the approximation of no energy losses, is implemented top of the Geant4 release 10.2. Both the new scattering physics and the built-in Coulomb scattering are used to reproduce the latest experimental results on grazing angle proton scattering. At 250 keV multiple scattering delivers large proton angles and it is not consistent with the observation. Among the tested models, the single scattering seems to better reproduce the scattering efficiency at the three energies but energy loss obtained at small scattering angles is significantly lower than the experimental values. In general, the energy losses obtained in the experiment are higher than what obtained by the simulation. The experimental data are not completely representative of the soft proton scattering experienced by current X-ray telescopes because of the lack of measurements at low energies (< 200 keV) and small reflection angles, so we are not able to address any of the tested models as the one that can certainly reproduce the scattering behavior of low energy protons expected for the ATHENA mission. We can, however, discard multiple scattering as the model able to reproduce soft proton funneling, and affirm that Coulomb single scattering can represent, until further measurements, the best approximation of the proton scattered angular distribution at the exit of X-ray optics.


title: Interference Mitigation via Relaying

abstract: This paper studies the effectiveness of relaying for interference mitigation in an interference-limited communication scenario. We are motivated by the observation that in a cellular network, a relay node placed at the cell edge observes a combination of intended signal and inter-cell interference that is correlated with the received signal at a nearby destination, so a relaying link can effectively allow the antennas at the relay and at the destination to be pooled together for both signal enhancement and interference mitigation. We model this scenario by a MIMO Gaussian relay channel with a digital relay-to-destination link of finite capacity, and with correlated noise across the relay and destination antennas. Assuming a compress-and-forward strategy with Gaussian input distribution and quantization noise, we propose a coordinate ascent algorithm for obtaining a stationary point of the non-convex joint optimization of the transmit and quantization covariance matrices. For fixed input distribution, the globally optimum quantization noise covariance matrix can be found in closed-form using a transformation of the relay's observation that simultaneously diagonalizes two conditional covariance matrices by congruence. For fixed quantization, the globally optimum transmit covariance matrix can be found via convex optimization. This paper further shows that such an optimized achievable rate is within a constant additive gap of the MIMO relay channel capacity. The optimal structure of the quantization noise covariance enables a characterization of the slope of the achievable rate as a function of the relay-to-destination link capacity. Moreover, this paper shows that the improvement in spatial degrees of freedom by MIMO relaying in the presence of noise correlation is related to the aforementioned slope via a connection to the deterministic relay channel.


title: Coordinate, Momentum and Dispersion operators in Phase space representation

abstract: The aim of this paper is to present a study on the representations of coordinate, momentum and dispersion operators in the framework of a phase space representation of quantum mechanics that we have introduced and studied in previous works. We begin in the introduction section with a recall about the concept of representation of operators on wave function spaces. Then, we show that in the case of the phase space representation the coordinate and momentum operators can be represented either with differential operators or with matrices. The explicit expressions of both the differential operators and matrices representations are established. Multidimensional generalization of the obtained results are performed and phase space representation of dispersion operators are given.


title: Possible signature of magnetic fields related to quasi-periodic oscillation observed in microquasars

abstract: The study of quasi-periodic oscillations (QPOs) of X-ray flux observed in the stellar-mass black hole binaries can provide a powerful tool for testing of the phenomena occurring in strong gravity regime. The high frequency (HF) QPOs usually come in pairs of double peaks which have a frequency ratio close to $3:2$. In addition to HF QPOs, some sources display simultaneous existence of the low frequency (LF) QPOs in Fourier power spectra. We demonstrate that the explanation of these phenomena can be well related to the epicyclic oscillations of charged particles in accretion disks orbiting Kerr black holes immersed in external large-scale magnetic fields. Magnetized versions of the standard geodesic models of QPOs can explain the observationally fixed data from the three microquasars. We perform a successful fitting of the HF QPOs observed in three microquasars, GRS 1915+105, XTE 1550-564 and GRO 1655-40, containing black holes, for magnetised versions of both epicyclic rezonance and relativistic precession models and discuss the corresponding constraints of parameters of the model, which are the mass and spin of black hole and the parameter related to the external magnetic field. Assuming the main source of synchrotron radiation producing X-rays are the relativistic electrons, we estimate the magnetic field in the vicinity of the black hole in the three sources to be of order $10^{-5}$ Gs which can serve as possible signature of the Galactic magnetic field magnitude. For heavier particles (protons, ions) larger magnetic fields are necessary for fitting the data.


title: Onset of Quantum Criticality in the Topological-to-Nematic Transition in a Two-dimensional Electron Gas at Filling Factor $\nu=5/2$

abstract: Under hydrostatic pressure, the ground state of a two-dimensional electron gas at $\nu=5/2$ changes from a fractional quantum Hall state to the stripe phase. By measuring the energy gap of the fractional quantum Hall state and of the onset temperature of the stripe phase we mapped out a phase diagram of these competing phases in the pressure-temperature plane. Our data highlight the dichotomy of two descriptions of the half-filled Landau level near the quantum critical point: one based on electrons and another on composite fermions.


title: Genuine equivariant operads

abstract: We build new algebraic structures, which we call genuine equivariant operads, which can be thought of as a hybrid between equivariant operads and coefficient systems. We then prove an Elmendorf-Piacenza type theorem stating that equivariant operads, with their graph model structure, are equivalent to genuine equivariant operads, with their projective model structure. As an application, we build explicit models for the $N_{\infty}$-operads of Blumberg and Hill.


title: Coloring Fibonacci-Cayley tree: An application to neural networks

abstract: This paper investigates the coloring problem on Fibonacci-Cayley tree, which is a Cayley graph whose vertex set is the Fibonacci sequence. More precisely, we elucidate the complexity of shifts of finite type defined on Fibonacci-Cayley tree via an invariant called entropy. It comes that computing the entropy of a Fibonacci tree-shift of finite type is equivalent to studying a nonlinear recursive system. After proposing an algorithm for the computation of entropy, we apply the result to neural networks defined on Fibonacci-Cayley tree, which reflect those neural systems with neuronal dysfunction. Aside from demonstrating a surprising phenomenon that there are only two possibilities of entropy for neural networks on Fibonacci-Cayley tree, we reveal the formula of the boundary in the parameter space.


title: Enhanced output entanglement with reservoir engineering

abstract: We study the output entanglement in a three-mode optomechanical system via reservoir engineering by shifting the center frequency of filter function away from resonant frequency. We find the bandwidth of the filter function can suppress the entanglement in the vicinity of resonant frequency of the system, while the entanglement will become prosperous if the center frequency departs from the resonant frequency. We obtain the approximate analytical expressions of the output entanglement, and from which we give the optimal center frequency at which the entanglement takes the maximum. Furthermore, we study the effects of time delay between the two output fields on the output entanglement, and obtain the optimal time delay for the case of large filter bandwidth.


title: A Lower Bound Technique for Communication in BSP

abstract: Communication is a major factor determining the performance of algorithms on current computing systems; it is therefore valuable to provide tight lower bounds on the communication complexity of computations. This paper presents a lower bound technique for the communication complexity in the bulk-synchronous parallel (BSP) model of a given class of DAG computations. The derived bound is expressed in terms of the switching potential of a DAG, that is, the number of permutations that the DAG can realize when viewed as a switching network. The proposed technique yields tight lower bounds for the fast Fourier transform (FFT), and for any sorting and permutation network. A stronger bound is also derived for the periodic balanced sorting network, by applying this technique to suitable subnetworks. Finally, we demonstrate that the switching potential captures communication requirements even in computational models different from BSP, such as the I/O model and the LPRAM.


title: Computational Models of Tutor Feedback in Language Acquisition

abstract: This paper investigates the role of tutor feedback in language learning using computational models. We compare two dominant paradigms in language learning: interactive learning and cross-situational learning - which differ primarily in the role of social feedback such as gaze or pointing. We analyze the relationship between these two paradigms and propose a new mixed paradigm that combines the two paradigms and allows to test algorithms in experiments that combine no feedback and social feedback. To deal with mixed feedback experiments, we develop new algorithms and show how they perform with respect to traditional knn and prototype approaches.


title: Black holes in massive gravity as heat engines

abstract: The paper at hand studies the heat engine provided by black holes in the presence of massive gravity. The main motivation is to investigate the effects of massive gravity on different properties of the heat engine. It will be shown that massive gravity parameters and graviton's mass modify the efficiency of engine on a significant level. Furthermore, it will be shown that it is possible to have the heat engine for non-spherical black holes in massive gravity and we study the effects of topological factor on properties of the heat engine. Surprisingly, it will be shown that the highest efficiency for the heat engine belongs to black holes with hyperbolic horizon, while the lowest one belongs to spherical black holes.


title: The structure and evolution of semiconducting buffer graphene grown on SiC(0001)

abstract: Using highly controlled coverages of graphene on SiC(0001), we have studied the structure of the first graphene layer that grows on the SiC interface. This layer, known as the buffer layer, is semiconducting. Using x-ray reflectivity and x-ray standing waves analysis we have performed a comparative study of the buffer layer structure with and without an additional monolayer graphene layer above it. We show that no more than 26\% of the buffer carbon is covalently bonded to Si in the SiC interface. We also show that the top SiC bilayer is Si depleted and is the likely the cause of the incommensuration previously observed in this system. When a monolayer graphene layer forms above the buffer, the buffer layer becomes less corrugated with signs of a change in the bonding geometry with the SiC interface. At the same time, the entire SiC interface becomes more disordered, presumably due to entropy associated with the higher growth temperature.


title: The Impact of Model Assumptions in Scalar-on-Image Regression

abstract: Complex statistical models such as scalar-on-image regression often require strong assumptions to overcome the issue of non-identifiability. While in theory it is well understood that model assumptions can strongly influence the results, this seems to be underappreciated, or played down, in practice. The article gives a systematic overview of the main approaches for scalar-on-image regression with a special focus on their assumptions. We categorize the assumptions and develop measures to quantify the degree to which they are met. The impact of model assumptions and the practical usage of the proposed measures are illustrated in a simulation study and in an application to neuroimaging data. The results show that different assumptions indeed lead to quite different estimates with similar predictive ability, raising the question of their interpretability. We give recommendations for making modeling and interpretation decisions in practice, based on the new measures and simulations using hypothetic coefficient images and the observed data.


title: How to observe the giant thermal effect in the Casimir force for graphene systems

abstract: A differential measurement scheme is proposed which allows for a clear observation of the giant thermal effect for the Casimir force, that was recently predicted to occur in graphene systems at short separation distances. The difference among the Casimir forces acting between a metal-coated sphere and the two halves of a dielectric plate, one uncoated and the other coated with graphene, is calculated in the framework of the Dirac model using the rigorous formalism of the polarization tensor. It is shown that in the proposed configuration both the difference among the Casimir forces and its thermal contributioncan be easily measured using already existing experimental setups. An observation of the giant thermal effect should open opportunities for modulation and control of dispersion forces in micromechanical systems based on graphene and other novel 2D-materials.


title: Measured Multiseries and Integration

abstract: A paper by Bruno Salvy and the author introduced measured multiseries and gave an algorithm to compute these for a large class of elementary functions, modulo a zero-equivalence method for constants. This gave a theoretical background for the implementation that Salvy was developing at that time. The main result of the present article is an algorithm to calculate measured multiseries for integrals of functions of the form h*sin G, where h and G belong to a Hardy field. The process can reiterated with the resulting algebra, and also applied to solutions of a second order differential equation of a particular form.