title: Error Control for Exponential Integration of the Master Equation

abstract: Error estimates for the numerical solution of the master equation are presented. Estimates are based on adjoint methods. We find that a good estimate can often be computed without spending computational effort on a dual problem. Estimates are applicable to both settings with time-independent, and time-dependent propensity functions. The Finite State Projection algorithm reduces the dimensionality of the problem and time propagation is based on an Arnoldi exponential integrator, which in the time-dependent setting is combined with a Magnus method. Local error estimates are devised for the truncation of both the Magnus expansion and the Krylov subspace in the Arnoldi algorithm. An issue with existing methods is that error estimates for truncation of the state space depend on measuring a loss of probability mass in a way that is not usually compatible with the approximation of the exponential. We suggest an alternative error estimate that is compatible with a Krylov approximation to the matrix exponential. Finally, we apply the new error estimates to develop an adaptive simulation algorithm. Numerical examples demonstrate the benefits of the approach.


title: Bounds for radii of starlikeness and convexity of some special functions

abstract: In this paper we consider some normalized Bessel, Struve and Lommel functions of the first kind, and by using the Euler-Rayleigh inequalities for the first positive zeros of combination of special functions we obtain tight lower and upper bounds for the radii of starlikeness of these functions. By considering two different normalization of Bessel and Struve functions we give some inequalities for the radii of convexity of the same functions. On the other hand, we show that the radii of univalence of some normalized Struve and Lommel functions are exactly the radii of starlikeness of the same functions. In addition, by using some ideas of Ismail and Muldoon we present some new lower and upper bounds for the zeros of derivatives of some normalized Struve and Lommel functions. The Laguerre-P\'olya class of real entire functions plays an important role in our study.


title: An approach to calculation of mass spectra and two-photon decays of $c\overline{c}$ mesons in the framework of Bethe-Salpeter Equation

abstract: In this work we calculate the mass spectrum of charmonium for $1P,...,3P$ states of $0^{++}$ , $1^{++}$, as well as for $1S,...,4S$ states of $0^{-+}$, and $1S, ...,5D$ states of $1^{--}$ along with the two-photon decay widths of ground and first excited state of $0^{++}$ quarkonia for the process, $O^{++}\rightarrow \gamma\gamma$ in the framework of a QCD motivated Bethe-Salpeter Equation. In this $4\times 4$ BSE framework, the coupled Salpeter equations are first shown to decouple for the confining part of interaction, under heavy-quark approximation, and analyically solved, and later the one-gluon-exchange interaction is perturbatively incorporated leading to mass spectral equations for various quarkonia. The analytic forms of wave functions obtained are used for calculation of two-photon decay widths of $\chi_{c0}$. Our results are in reasonable agreement with data (where ever available) and other models.


title: Model-Free Closed-Loop Stability Analysis: A Linear Functional Approach

abstract: Performing a stability analysis during the design of any electronic circuit is critical to guarantee its correct operation. A closed-loop stability analysis can be performed by analysing the impedance presented by the circuit at a well-chosen node without internal access to the simulator. If any of the poles of this impedance lie in the complex right half-plane, the circuit is unstable. The classic way to detect unstable poles is to fit a rational model on the impedance. In this paper, a projection-based method is proposed which splits the impedance into a stable and an unstable part by projecting on an orthogonal basis of stable and unstable functions. When the unstable part lies significantly above the interpolation error of the method, the circuit is considered unstable. Working with a projection provides one, at small cost, with a first appraisal of the unstable part of the system. Both small-signal and large-signal stability analysis can be performed with this projection-based method. In the small-signal case, a low-order rational approximation can be fitted on the unstable part to find the location of the unstable poles.


title: Measuring Structural Parameters Through Stacking Galaxy Images

abstract: It remains challenging to detect the low surface brightness structures of faint high-z galaxies, which is key to understanding the structural evolution of galaxies. The technique of image stacking allows us to measure the averaged light profile beneath the detection limit and probe the extended structure of a group of galaxies. We carry out simulations to examine the recovery of the averaged surface brightness profile through stacking model HST/ACS images of a set of galaxies as functions of Sersic index (n), effective radius (Re) and axis ratio (AR). The Sersic profile best fitting the radial profile of the stacked image is taken as the recovered profile, in comparison with the intrinsic mean profile of the model galaxies. Our results show that, in general, the structural parameters of the mean profile can be properly determined through stacking, although systematic biases need to be corrected when spreads of Re and AR are counted. We find that Sersic index is slightly overestimated and Re is underestimated at AR < 0.5 as the stacked image appears to be more compact due to the presence of inclined galaxies; the spread of Re biases the stacked profile to have a higher Sersic index. We stress that the measurements of structural parameters through stacking should take these biases into account. We estimate the biases in the recovered structural parameters from stacks of galaxies when the samples have distributions of Re, AR and n seen in local galaxies.


title: Uniqueness of scatterer in inverse acoustic obstacle scattering with a single incident plane wave

abstract: In this paper, we give a positive answer to a challenging open problem for recovering unknown obstacle (which is usually referred to as a scatterer) by acoustic wave probe associated to the Helmholtz equation. We show that the acoustic scattering amplitude $A(\beta, {\alpha}_0, k_0)$, known for all ${\beta}\in {\mathbb S}^2$, where ${\mathbb {S}}^2$ is the unit sphere in ${\mathbb{R}}^3$, ${{\alpha}}_0\in {\mathbb{S}}^2$ is fixed, $k_0>0$ is fixed, determines the obstacle $D$ and the boundary condition on $\partial D$ uniquely (The boundary condition on $\partial D$ is either the Dirichlet, or Neumann, or the impedance one).


title: Elastohydrodynamic Synchronization of Adjacent Beating Flagella

abstract: It is now well established that nearby beating pairs of eukaryotic flagella or cilia typically synchronize in phase. A substantial body of evidence supports the hypothesis that hydrodynamic coupling between the active filaments, combined with waveform compliance, provides a robust mechanism for synchrony. This elastohydrodynamic mechanism has been incorporated into `bead-spring' models in which flagella are represented by microspheres tethered by radial springs as internal forces drive them about orbits. While these low-dimensional models reproduce the phenomenon of synchrony, their parameters are not readily relatable to those of flagella. More realistic models which reflect the elasticity of the axonemes and active force generation take the form of fourth-order nonlinear PDEs. While computational studies have shown synchrony, the effects of hydrodynamic coupling between nearby filaments governed by such models have been theoretically examined only in the regime of interflagellar distances $d$ large compared to flagellar length $L$. Yet, in many biological situations $d/L \ll 1$. Here, we first present an asymptotic analysis of the hydrodynamic coupling between two filaments in the regime $d/L \ll 1$, and find that the form of the coupling is independent of the details of the internal forces that govern the motion of the filaments. The analysis is like the localized induction approximation for vortex filament motion, extended to the case of mutual induction. To understand how the coupling mechanism leads to synchrony of extended objects, we introduce a heuristic model of flagellar beating, a single fourth-order nonlinear PDE whose form is derived from symmetry considerations, the physics of elasticity, and the overdamped nature of the dynamics. Analytical and numerical studies of this model illustrate how synchrony between two filaments is achieved through the asymptotic coupling.


title: Quantum frequency conversion between infrared and ultraviolet

abstract: We report on the implementation of quantum frequency conversion (QFC) between infrared (IR) and ultraviolet (UV) wavelengths by using single-stage upconversion in a periodically poled KTP waveguide. Due to the monolithic waveguide design, we manage to transfer a telecommunication band input photon to the wavelength of the ionic dipole transition of Yb${}^{+}$ at 369.5 nm. The external (internal) conversion efficiency is around 5\% (10\%). The high energy pump used in this converter introduces a spontaneous parametric downconversion (SPDC) process, which is a cause for noise in the UV mode. Using this SPDC process, we show that the converter preserves non-classical correlations in the upconversion process, rendering this miniaturized interface a source for quantum states of light in the UV.


title: Determining pitch-angle diffusion coefficients from test particle simulations

abstract: Transport and acceleration of charged particles in turbulent media is a topic of great interest in space physics and interstellar astrophysics. These processes are dominated by the scattering of particles off magnetic irregularities. The scattering process itself is usually described by small-angle scattering with the pitch-angle coefficient $D_{\mu\mu}$ playing a major role. Since the diffusion coefficient $D_{\mu\mu}$ can be determined analytically only for the approximation of quasi-linear theory, the determination of this coefficient from numerical simulations has, therefore, become more important. So far these simulations yield particle tracks for small-scale scattering, which can then be interpreted using the running diffusion coefficients. This method has a limited range of validity. This paper presents two new methods that allow for the calculation of the pitch-angle diffusion coefficient from numerical simulations. These methods no longer analyse particle trajectories, but the change of particle distribution functions. It is shown that they provide better resolved results and allow for the analysis of strong turbulence. The application of these methods to Monte Carlo simulations of particle scattering and hybrid MHD-particle simulations is presented. Both analysis methods are able to recover the diffusion coefficients used as input for the Monte Carlo simulations and provide better results in MHD simulations especially for stronger turbulence.


title: Braids, orderings and minimal volume cusped hyperbolic 3-manifolds

abstract: It is well-known that there is a faithful representation of braid groups on automorphism groups of free groups, and it is also well-known that free groups are bi-orderable. We investigate which n-strand braids give rise to automorphisms which preserve some bi-ordering of the free group rank n. As a consequence of our work we find that of the two minimal volume hyperbolic 2-cusped orientable 3-manifolds, one has bi-orderable fundamental group whereas the other does not. We prove a similar result for the 1-cusped case, and have further results for more cusps. In addition, we study pseudo-Anosov braids and find that typically those with minimal dilatation are not order-preserving.


title: Solitons of a vector model on the honeycomb lattice

abstract: We study a simple nonlinear vector model defined on the honeycomb lattice. We propose a bilinearization scheme for the field equations and demonstrate that the resulting system is closely related to the well-studied integrable models, such as the Hirota bilinear difference equation and the Ablowitz-Ladik system. This result is used to derive the N-soliton solutions.


title: Evidence of AGB pollution in Galactic globular clusters from the Mg-Al anticorrelations observed by the APOGEE survey

abstract: We study the formation of multiple populations in globular clusters (GC), under the hypothesis that stars in the second generation formed from the winds of intermediate-mass stars, ejected during the asymptotic giant branch (AGB) phase, possibly diluted with pristine gas, sharing the same chemical composition of first-generation stars. To this aim, we use the recent APOGEE data, which provide the surface chemistry of a large sample of giant stars, belonging to clusters that span a wide metallicity range. The APOGEE data set is particularly suitable to discriminate among the various pollution scenarios proposed so far, as it provides the surface abundances of Mg and Al, the two elements involved in a nuclear channel extremely sensitive to the temperature, hence to the metallicity of the polluters. The present analysis shows a remarkable agreement between the observations and the theoretical yields from massive AGB stars. In particular, the observed extension of the depletion of Mg and O and the increase in Al is well reproduced by the models and the trend with the metallicity is also fully accounted for. This study further supports the idea that AGB stars were the key players in the pollution of the intra-cluster medium, from which additional generations of stars formed in GC.


title: Partonic structure of proton in the resonance region

abstract: We separate the contributions of parton distributions from higher twist corrections to the deeply inelastic lepton-proton scattering in the resonance region using the Jefferson Lab data at low $Q^2$. The study indicates that the concept of the valence quarks and their distributions are indispensable even at $Q^2<1GeV^2$. The quark-hadron duality is also discussed.


title: A Flat and Scalable Data Center Network Topology Based on De Bruijn Graphs

abstract: Due to the requirement of hosting tens of thousands of hosts in today's data centers, data center networks strive for scalability and high throughput on the one hand. On the other hand, the cost for networking hardware should be minimized. Consequently, the number and complexity (e.g. TCAM size) of switches has to be minimized. These requirements led to network topologies like Clos and Leaf-Spine networks only requiring a shallow hierarchy of switches---two levels for Leaf-Spine networks. The drawback of these topologies is that switches at higher levels like Spine switches need a high port density and, thus, are expensive and limit the scalability of the network. In this paper, we propose a data center network topology based on De Bruijn graphs completely avoiding a switch hierarchy and implementing a flat network topology of top-of-rack switches instead. This topology guarantees logarithmic (short) path length. We show that the required routing logic can be implemented by standard prefix matching operations in hardware (TCAM) allowing for using commodity switches without any modification. Moreover, forwarding requires only a very small number of forwarding table entries, saving costly and energy-intensive TCAM.


title: Toward a new instances of NELL

abstract: We are developing the method to start new instances of NELL in various languages and develop then NELL multilingualism. We base our method on our experience on NELL Portuguese and NELL French. This reports explain our method and develops some research perspectives.


title: Strongly interacting matter from holographic QCD model

abstract: We introduce the 5-dimension dynamical holographic QCD model, which is constructed in the graviton-dilaton-scalar framework with the dilaton background field $\Phi$ and the scalar field $X$ responsible for the gluodynamics and chiral dynamics, respectively. We review our results on the hadron spectra including the glueball and light meson spectra, QCD phase transitions and transport properties in the framework of the dynamical holographic QCD model.


title: Quantum state transfer with ultracold atoms in optical lattices

abstract: Ultracold atoms can be used to perform quantum simulations of a variety of condensed matter systems, including spin systems. These progresses point to the implementation of the manipulation of quantum states and to observe and exploit the effect of quantum correlations. A natural direction along this line is provided by the possibility to perform quantum state transfer (QST). After presenting a brief discussion of the simulation of quantum spin chains with ultracold gases and reminding the basic facts of QST, we discuss how to potentially use the tools of present-day ultracold technology to implement the QST between two regions of the atomic system (the sender and the receiver). The fidelity and the typical timescale of the QST are discussed, together with possible limitations and applications of the presented results.


title: Runaway-electron formation and electron slide-away in an ITER post-disruption scenario

abstract: Mitigation of runaway electrons is one of the outstanding issues for the reliable operation of ITER and other large tokamaks, and accurate estimates for the expected runaway-electron energies and current are needed. Previously, linearized tools (which assume the runaway population to be small) have been used to study the runaway dynamics, but these tools are not valid in the cases of most interest, i.e. when the runaway population becomes substantial. We study runaway-electron formation in a post-disruption ITER plasma using the newly developed non-linear code NORSE, and describe a feedback mechanism by which a transition to electron slide-away can be induced at field strengths significantly lower than previously expected. If the electric field is actively imposed using the control system, the entire electron population is quickly converted to runaways in the scenario considered. We find the time until the feedback mechanism sets in to be highly dependent on the details of the mechanisms removing heat from the thermal electron population.


title: Searches for the FCNC couplings from top-Higgs associated production signal with $h\to \gamma\gamma$ at the LHC

abstract: In this paper, we study the observability of the top-Higgs flavor changing neutral current (FCNC) $tqh$ coupling through the process $pp\to qg\to t(\to \ell^{+}b\nu)h(\to \gamma\gamma)$ at the Large Hadron Collider~(LHC), where $\ell=e, \mu$. Our numerical results show that, in some parameter regions, the LHC may observe the above signals at the $5\sigma$ level. Otherwise, the branching ratios $Br(t\to uh)$ and $Br(t\to ch)$ can be respectively probed to $0.036\%$ and $0.13\%$ at $3\sigma$ level at 14 TeV LHC with the high integrated luminosity of 3000 fb$^{-1}$. On the other hand, studying the charge ratio for the lepton in top quark decay can be not only used to discriminate between signal and backgrounds, but also used to discriminate between $tuh$ and $tch$ couplings, for which anomalous single top production comes from the up initiated channel and charm initiated channel.


title: Supersymmetry in Open Superstring Field Theory

abstract: We realize the 16 unbroken supersymmetries on a BPS D-brane as invariances of the action of the corresponding open superstring field theory. We work in the small Hilbert space approach, where a symmetry of the action translates into a symmetry of the associated cyclic $A_\infty$ structure. We compute the supersymmetry algebra, being careful to disentangle the components which produce a translation, a gauge transformation, and a symmetry transformation which vanishes on-shell. Via the minimal model theorem, we illustrate how supersymmetry of the action implies supersymmetry of the tree level open string scattering amplitudes.


title: Resonance poles and threshold energies for hadron physical problems by a model-independent universal algorithm

abstract: We show how complex resonance poles and threshold energies for systems in hadron physics can be accurately obtained by using a method based on the Pad\'{e}-approximant which was recently developed for the calculation of resonance poles for atomic and molecular auto-ionization systems. The main advantage of this method is the ability to calculate the resonance poles and threshold energies from \emph{real} spectral data. In order to demonstrate the capabilities of this method we apply it here to an analytical model as well as to experimental data for the squared modulus of the vector pion form factor, the S0 partial wave amplitude for $\pi\pi$ scattering and the cross section ratio $R(s)$ for $e^+e^-$ collisions. The extracted values for the resonance poles of the $\rho(770)$ and the $f_0(500)$ or $\sigma$ meson are in very good agreement with the literature. When the data are noisy the prediction of decay thresholds proves to be less accurate but feasible.


title: A quantum spectrum analyzer enhanced by a nuclear spin memory

abstract: We realize a two-qubit sensor designed for achieving high spectral resolution in quantum sensing experiments. Our sensor consists of an active "sensing qubit" and a long-lived "memory qubit", implemented by the electronic and the nitrogen-15 nuclear spins of a nitrogen-vacancy center in diamond, respectively. Using state storage times of up to 45 ms, we demonstrate spectroscopy of external ac signals with a line width of 19 Hz (~2.9 ppm) and of carbon-13 nuclear magnetic resonance (NMR) signals with a line width of 190 Hz (~ 74 ppm). This represents an up to 100-fold improvement in spectral resolution compared to measurements without nuclear memory.


title: Fast sausage modes in magnetic tubes with continuous transverse profiles: effects of a finite plasma beta

abstract: While standing fast sausage modes in flare loops are often invoked to interpret quasi-periodic pulsations (QPPs) in solar flares, it is unclear as to how they are influenced by the combined effects of a continuous transverse structuring and a finite internal plasma beta ($\beta_{\rm i}$). We derive a generic dispersion relation (DR) governing linear sausage waves in straight magnetic tubes for which plasma pressure is not negligible and the density and temperature inhomogeneities of essentially arbitrary form take place in a layer of arbitrary width. Focusing on fast modes, we find that $\beta_{\rm i}$ only weakly influences $k_{\rm c}$, the critical longitudinal wavenumber separating the leaky from trapped modes. Likewise, for both trapped and leaky modes, the periods $P$ in units of the transverse fast time depend only weakly on $\beta_{\rm i}$, which is compatible with the fact that the effective wavevectors of fast sausage modes are largely perpendicular to the background magnetic field. However, a weak $\beta_{\rm i}$ dependence of the damping times $\tau$ is seen only when the length-to-radius ratio $L/R$ is $\sim 50\%$ larger than some critical value $\pi/(k_{\rm c} R)$, which itself rather sensitively depends on the density contrast, profile steepness as well as on how the transverse structuring is described. In the context of QPPs, we conclude that the much simpler zero-beta theory can be employed for trapped modes, as long as one sees the deduced internal Alfv\'en speed as actually being the fast speed. In contrast, effects due to a finite beta in flare loops should be considered when leaky modes are exploited.


title: Finiteness of gibbs measures on noncompact manifolds with pinched negative curvature

abstract: We characterize the finiteness of Gibbs measures for geodesic flows on negatively curved manifolds by several criteria, analogous to those proposed by Sarig for symbolic dynamical systems over an infinite alphabet. As an application, we recover Dal'bo-Otal-Peign{\'e} criterion of finiteness for the Bowen-Margulis measure on geometrically finite hyperbolic manifolds, as well as Peign{\'e}' examples of gemetrically infinite manifolds having a finite Bowen-Margulis measure. These criteria should be useful in the future to find more examples with finite Gibbs measures.


title: GMM-Free Flat Start Sequence-Discriminative DNN Training

abstract: Recently, attempts have been made to remove Gaussian mixture models (GMM) from the training process of deep neural network-based hidden Markov models (HMM/DNN). For the GMM-free training of a HMM/DNN hybrid we have to solve two problems, namely the initial alignment of the frame-level state labels and the creation of context-dependent states. Although flat-start training via iteratively realigning and retraining the DNN using a frame-level error function is viable, it is quite cumbersome. Here, we propose to use a sequence-discriminative training criterion for flat start. While sequence-discriminative training is routinely applied only in the final phase of model training, we show that with proper caution it is also suitable for getting an alignment of context-independent DNN models. For the construction of tied states we apply a recently proposed KL-divergence-based state clustering method, hence our whole training process is GMM-free. In the experimental evaluation we found that the sequence-discriminative flat start training method is not only significantly faster than the straightforward approach of iterative retraining and realignment, but the word error rates attained are slightly better as well.


title: Catastrophic thinning of dielectric elastomers

abstract: We provide a clear energetic insight into the catastrophic nature of the so-called creasing and pull-in instabilities in soft electro-active elastomers. These phenomena are ubiquitous for thin electro-elastic plates and are a major obstacle to the development of giant actuators; yet they are not completely understood nor modelled accurately. Here, in complete agreement with experiments, we give a simple formula to predict the voltage thresholds for these instability patterns and model their shape, and show that equilibrium is impossible beyond their onset. Our analysis is fully analytical, does not require finite element simulations, and can be extended to include pre-stretch and to encompass any material behaviour.


title: CTA in the Context of Searches for Particle Dark Matter - a glimpse

abstract: In this contribution, CTAs potential role in detection of particle dark matter in the context of other detection approaches is briefly discussed for an audience of gamma-ray astronomers. In particular searches for new particles at the large hadron collider and detection of dark matter particles in deep underground detectors are considered. We will focus on Weakly Interacting Massive Particles (WIMP). Approaches will be compared in terms of (a) robustness of sensitivity predictions, (b) timeline and (c) reach. The estimate of the reach will be model-dependent. Given our ignorance about the nature of dark matter, and the complementarity of detection techniques even within a given framework (e.g. Supersymmetry), the trivial conclusion is that we might need all approaches and the most sensitive experiments. Our discussion will be somewhat more restrictive in order to be able to be more concrete. With the caveat of incompleteness, under the assumption that the WIMP paradigm describes nature, CTA is more likely to discover multi-TeV WIMP dark matter, whereas for lower masses direct detection and LHC has significantly better prospects. We will illustrate this conclusion with examples from foremost Supersymmetry, but mention effective field theory or simplified models. We will comment on a few models predicting high mass WIMPs, in particular 1 TeV higgsino and wino WIMPs, as well as Minimal Dark Matter and point out the relevance of updated measurements of the anomalous magnetic moment of the muon for CTAs role in searches for Supersymmetry.


title: Epidemics of Liquidity Shortages in Interbank Markets

abstract: Financial contagion from liquidity shocks has being recently ascribed as a prominent driver of systemic risk in interbank lending markets. Building on standard compartment models used in epidemics, in this work we develop an EDB (Exposed-Distressed-Bankrupted) model for the dynamics of liquidity shocks reverberation between banks, and validate it on electronic market for interbank deposits data. We show that the interbank network was highly susceptible to liquidity contagion at the beginning of the 2007/2008 global financial crisis, and that the subsequent micro-prudential and liquidity hoarding policies adopted by banks increased the network resilience to systemic risk---yet with the undesired side effect of drying out liquidity from the market. We finally show that the individual riskiness of a bank is better captured by its network centrality than by its participation to the market, along with the currently debated concept of "too interconnected to fail".


title: Channel Training for Analog FDD Repeaters: Optimal Estimators and Cram\'er-Rao Bounds

abstract: For frequency division duplex channels, a simple pilot loop-back procedure has been proposed that allows the estimation of the UL & DL channels at an antenna array without relying on any digital signal processing at the terminal side. For this scheme, we derive the maximum likelihood (ML) estimators for the UL & DL channel subspaces, formulate the corresponding Cram\'er-Rao bounds and show the asymptotic efficiency of both (SVD-based) estimators by means of Monte Carlo simulations. In addition, we illustrate how to compute the underlying (rank-1) SVD with quadratic time complexity by employing the power iteration method. To enable power control for the data transmission, knowledge of the channel gains is needed. Assuming that the UL & DL channels have on average the same gain, we formulate the ML estimator for the channel norm, and illustrate its robustness against strong noise by means of simulations.


title: Propagation of chaos for aggregation equations with no-flux boundary conditions and sharp sensing zones

abstract: We consider a $N$-particle interacting particle system with the vision geometrical constraints and reflected noises, proposed as a model for collective behavior of individuals. We rigorously derive a continuity-type of mean-field equation with discontinuous kernels and the normal reflecting boundary conditions from that stochastic particle system as the number of particles $N$ goes to infinity. More precisely, we provide a quantitative estimate of the convergence in law of the empirical measure associated to the particle system to a probability measure which possesses a density which is a weak solution to the continuity equation. This extends previous results on an interacting particle system with bounded and Lipschitz continuous drift terms and normal reflecting boundary conditions by Sznitman[J. Funct. Anal., 56, (1984), 311--336] to that one with discontinuous kernels.


title: Model Reduction for Systems with Inhomogeneous Initial Conditions

abstract: We consider the model reduction problem for linear time-invariant dynamical systems having nonzero (but otherwise indeterminate) initial conditions. Building upon the observation that the full system response is decomposable as a superposition of the response map for an unforced system having nontrivial initial conditions and the response map for a forced system having null initial conditions, we develop a new approach that involves reducing these component responses independently and then combining the reduced responses into an aggregate reduced system response. This approach allows greater flexibility and offers better approximation properties than other comparable methods.


title: Error Asymmetry in Causal and Anticausal Regression

abstract: It is generally difficult to make any statements about the expected prediction error in an univariate setting without further knowledge about how the data were generated. Recent work showed that knowledge about the real underlying causal structure of a data generation process has implications for various machine learning settings. Assuming an additive noise and an independence between data generating mechanism and its input, we draw a novel connection between the intrinsic causal relationship of two variables and the expected prediction error. We formulate the theorem that the expected error of the true data generating function as prediction model is generally smaller when the effect is predicted from its cause and, on the contrary, greater when the cause is predicted from its effect. The theorem implies an asymmetry in the error depending on the prediction direction. This is further corroborated with empirical evaluations in artificial and real-world data sets.


title: Observations of Binary Systems with the H.E.S.S. Telescopes

abstract: Observations of binary systems obtained recently with the High Energy Stereoscopic System (H.E.S.S) of Cherenkov telescopes are reported. The outcomes of a detailed observation campaign on PSR B1259-63 during its periastron passage in 2014 will be presented. This system was observed for the first time with H.E.S.S. II, providing spectra and light curves down to 200 GeV, which will be compared with observations conducted during previous periastron passages and with results from an analysis of contemporaneously taken Fermi-LAT data. Also long-term observations of LS 5039 with H.E.S.S in phase I and phase II are reported. This source was monitored at very high energies (VHEs) in a period of time spanning more than ten years. Its spectral energy distribution measured with H.E.S.S. II extends down to 120 GeV. Spectral results from the Fermi-LAT observations are shown as well, and the compatibility with H.E.S.S. results in the overlapping energy range is discussed. The identification of the new gamma-ray binary candidate HESS J1832-093 will also be presented. Furthermore, the search for VHE emission from the microquasars GRS 1915+105, Circinus X-1 and V4641 Sgr based on data from H.E.S.S. observations conducted contemporaneously with the RXTE satellite experiment will be reported. These data provide constraints on the integral gamma-ray flux at different X-ray states of the three sources.


title: Lower bounds on the size of general Schr\"odinger-cat states from experimental data

abstract: Experimental progress with meso- and macroscopic quantum states (i.e., general Schrodinger-cat states) was recently accompanied by theoretical proposals on how to measure the merit of these efforts. So far, experiment and theory were disconnected as theoretical analysis of actual experimental data was missing. Here, we consider a proposal for macroscopic quantum states that measures the extend of quantum coherence present in the system. For this, the quantum Fisher information is used. We calculate lower bounds from real experimental data. The results are expressed as an "effective size", that is, relative to "classical" reference states. We find remarkable numbers of up to 70 in photonic and atomic systems.


title: On the Optimality of Tape Merge of Two Lists with Similar Size

abstract: The problem of merging sorted lists in the least number of pairwise comparisons has been solved completely only for a few special cases. Graham and Karp \cite{taocp} independently discovered that the tape merge algorithm is optimal in the worst case when the two lists have the same size. In the seminal papers, Stockmeyer and Yao\cite{yao}, Murphy and Paull\cite{3k3}, and Christen\cite{christen1978optimality} independently showed when the lists to be merged are of size $m$ and $n$ satisfying $m\leq n\leq\lfloor\frac{3}{2}m\rfloor+1$, the tape merge algorithm is optimal in the worst case. This paper extends this result by showing that the tape merge algorithm is optimal in the worst case whenever the size of one list is no larger than 1.52 times the size of the other. The main tool we used to prove lower bounds is Knuth's adversary methods \cite{taocp}. In addition, we show that the lower bound cannot be improved to 1.8 via Knuth's adversary methods. We also develop a new inequality about Knuth's adversary methods, which might be interesting in its own right. Moreover, we design a simple procedure to achieve constant improvement of the upper bounds for $2m-2\leq n\leq 3m $.


title: Mixed Hodge structures and formality

abstract: We prove a formality theorem for algebraic objects internal to smooth complex varieties that are not compact but whose mixed Hodge structure has a certain purity property.


title: Ratchet Baryogenesis with an Analogy to the Forced Pendulum

abstract: A new scenario of baryogenesis via the ratchet mechanism is proposed based on an analogy to the forced pendulum. The oscillation of the inflaton field during the reheating epoch after inflation plays the role of the driving force, while the phase $\theta$ of a scalar baryon field (a complex scalar field with baryon number) plays the role of the angle of the pendulum. When the inflaton is coupled to the scalar baryon, the behavior of the phase $\theta$ can be analogous to that of the angle of the forced pendulum. If the oscillation of the driving force is adjusted to the pendulum's motion, a directed rotation of the pendulum is obtained with a non-vanishing value of $\dot{\theta}$, which models successful baryogenesis since $\dot{\theta}$ is proportional to the baryon number density. Similar ratchet models which lead to directed motion have been used in the study of molecular motors in biology. There, the driving force is supplied by chemical reactions, while in our scenario this role is played by the inflaton during the reheating epoch. Keywords: Baryogenesis; Reheating; Inflaton; Ratchet Model; Forced Rotating Pendulum


title: Separate chemical freeze-outs of strange and non-strange hadrons and problem of residual chemical non-equilibrium of strangeness in relativistic heavy ion collisions

abstract: We present an elaborate version of the hadron resonance gas model with the combined treatment of separate chemical freeze-outs for strange and non-strange hadrons and with an additional $\gamma_{s}$ factor which accounts for the remaining strange particle non-equilibration. Within suggested approach the parameters of two chemical freeze-outs are connected by the conservation laws of entropy, baryonic charge, third isospin projection and strangeness. The developed model enables us to perform a high-quality fit of the hadron multiplicity ratios measured at AGS, SPS and RHIC with $\chi^2/dof \simeq 0.93$. A special attention is paid to a successful description of the Strangeness Horn. The well-known problem of selective suppression of $\bar \Lambda $ and $\bar \Xi$ hyperons is also discussed. The main result is that for all collision energies the $\gamma_{s}$ factor is about 1 within the error bars, except for the center of mass collision energy 7.6 GeV at which we find about 20\% enhancement of strangeness. Also we confirm an existence of strong jumps in pressure, temperature and effective number of degrees of freedom at the stage of strange particle chemical freeze-out, when the center of mass collision energy changes from 4.3 to 4.9 GeV. We argue that these irregularities may signal about the quark-hadron phase transition.


title: Schauder estimates at the boundary for sub-laplacians in Carnot groups

abstract: In this paper we prove Schauder estimates at the boundary for sub-Laplacian type operators in Carnot groups. While internal Schauder estimates have been deeply studied, up to now subriemannian estimates at the boundary are known only in the Heisenberg groups. The proof of these estimates in the Heisenberg setting, due to Jerison, is based on the Fourier transform technique and can not be repeated in general Lie groups. After the result of Jerison no new contribution to the boundary problem has been provided. In this paper we introduce a new approach, which allows to built a Poisson kernel starting from the fundamental solution, from which we deduce the Schauder estimates at non characteristic boundary points.


title: A highly accurate {\it ab initio} potential energy surface for methane

abstract: A new nine-dimensional potential energy surface (PES) for methane has been generated using state-of-the-art \textit{ab initio} theory. The PES is based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set limit and incorporates a range of higher-level additive energy corrections. These include: core-valence electron correlation, higher-order coupled cluster terms beyond perturbative triples, scalar relativistic effects and the diagonal Born-Oppenheimer correction. Sub-wavenumber accuracy is achieved for the majority of experimentally known vibrational energy levels with the four fundamentals of $^{12}$CH$_4$ reproduced with a root-mean-square error of $0.70{\,}$cm$^{-1}$. The computed \textit{ab initio} equilibrium C{--}H bond length is in excellent agreement with previous values despite pure rotational energies displaying minor systematic errors as $J$ (rotational excitation) increases. It is shown that these errors can be significantly reduced by adjusting the equilibrium geometry. The PES represents the most accurate \textit{ab initio} surface to date and will serve as a good starting point for empirical refinement.


title: Quantum coherence in the Dynamical Casimir Effect

abstract: We propose to use quantum coherence as ultimate proof of quantum nature in the radiation that appears by means of the Dynamical Casimir Effect in experiments with superconducting microwave waveguides. We show that, unlike previously considered measurements of quantum correlations such as entanglement and discord, quantum coherence does not require a threshold value of the external pump amplitude and is highly robust to thermal noise.


title: Fabrication of sharp atomic force microscope probes using in-situ local electric field induced deposition under ambient conditions

abstract: We demonstrate a simple method to significantly improve the sharpness of standard silicon probes for an atomic force microscope, or to repair a damaged probe. The method is based on creating and maintaining a strong, spatially localized electric field in the air gap between the probe tip and the surface of conductive sample. Under these conditions, nanostructure growth takes place on both the sample and the tip. The most likely mechanism is the decomposition of atmospheric adsorbate with subsequent deposition of carbon structures. This makes it possible to grow a spike of a few hundred nanometers in length on the tip. We further demonstrate that probes obtained by this method can be used for high-resolution scanning. It is important to note that all process operations are carried out in-situ, in air and do not require the use of closed chambers or any additional equipment beyond the atomic force microscope itself.


title: Classic balayage of measures and subharmonic functions to a system of rays

abstract: We construct and apply the classic balayage (sweeping out) of measures and subharmonic functions on closed system of rays in the complex plane with vertex at the origin, including measures and subharmonic functions and infinite order. The need for such a procedure occurs in the study of the behavior of entire and subharmonic functions on systems of rays. The results apply to the complete regularity of the growth of subharmonic and entire functions on a system of rays, the study of the distribution of zeros of entire functions of exponential type of class A, touches upon the incompleteness of exponential systems. The paper contains a survey character sections and new results.


title: 5D Maximally Supersymmetric Yang-Mills on the Lattice

abstract: We provide details of the lattice construction of five-dimensional maximally supersymmetric Yang-Mills theory. The lattice theory is supersymmetric, gauge invariant and free from spectrum doublers. Such a supersymmetric lattice formulation is interesting as it can be used for non-perturbative explorations of the five-dimensional theory, which has a known gravitational dual.


title: Assisted Dictionary Learning for fMRI Data Analysis

abstract: Extracting information from functional magnetic resonance (fMRI) images has been a major area of research for more than two decades. The goal of this work is to present a new method for the analysis of fMRI data sets, that is capable to incorporate a priori available information, via an efficient optimization framework. Tests on synthetic data sets demonstrate significant performance gains over existing methods of this kind.


title: Converse theorems and the local Langlands correspondence in families

abstract: We prove a descent criterion for certain families of smooth representations of GL_n(F) (F a p-adic field) in terms of the gamma factors of pairs constructed in previous work of the second author. We then use this descent criterion, together with a theory of gamma factors for families of representations of the Weil group W_F (developed previously by both authors), to prove a series of conjectures, due to the first author, that give a complete description of the integral Bernstein center in terms of Galois theory and the local Langlands correspondence. An immediate consequence is the conjectural "local Langlands correspondence in families" of Emerton and Helm.


title: On strict convergence of stochastic gradients

abstract: We discuss conditions ensuring the (strict) convergence of stochastic gradient algorithms.


title: $\mathcal H_2$-Quasi-Optimal Model Order Reduction for Quadratic-Bilinear Control Systems

abstract: We investigate the optimal model reduction problem for large-scale quadratic-bilinear (QB) control systems. Our contributions are threefold. First, we discuss the variational analysis and the Volterra series formulation for QB systems. We then define the $\mathcal H_2$-norm for a QB system based on the kernels of the underlying Volterra series and also propose a truncated $\mathcal H_2$-norm. Next, we derive first-order necessary conditions for an optimal approximation, where optimality is measured in term of the truncated $\mathcal H_2$-norm of the error system. We then propose an iterative model reduction algorithm, which upon convergence yields a reduced-order system that approximately satisfies the newly derived optimality conditions. We also discuss an efficient computation of the reduced Hessian, using the special Kronecker structure of the Hessian of the system. We illustrate the efficiency of the proposed method by means of several numerical examples resulting from semi-discretized nonlinear partial differential equations and show its competitiveness with the existing model reduction schemes for QB systems such as moment-matching methods and balanced truncation.


title: Guidelines for designing 2D and 3D plasmonic stub resonators

abstract: In this work, we compare the performance of plasmonic waveguide integrated stub resonators based on 2D metal-dielectric-metal and 3D slot waveguide (SWG) geometries. We show that scattering matrix theory can be extended to 3D devices, and by employing scattering matrix theory, we provide the guidelines for designing plasmonic 2D and 3D single-stub and double-stub resonators with a desired spectral response at the design wavelength. We provide transmission maps of 2D and 3D double-stub resonators versus stub lengths, and we specify the different regions on these maps that result in a minimum, a maximum, or a plasmonically induced transparency shape in the transmission spectrum. Radiation loss from waveguide terminations leads to a degradation of the 3D SWG-based resonators. We illustrate improved waveguide terminations that boost resonator properties. We verify our results with 3D finite-difference time-domain (FDTD) simulations.


title: Noether's theorems and conserved currents in gauge theories in the presence of fixed fields

abstract: We extend the standard construction of conserved currents for matter fields in general relativity to general gauge theories. In the original construction the conserved current associated with a spacetime symmetry generated by a Killing field $h^\mu$ is given by $\sqrt{-g}\,T^{\mu\nu}h_\nu$, where $T^{\mu\nu}$ is the energy-momentum tensor of the matter. We show that if in a Lagrangian field theory that has gauge symmetry in the general Noetherian sense some of the elementary fields are fixed and are invariant under a particular infinitesimal gauge transformation, then there is a current $\mathcal{B}^\mu$ that is analogous to $\sqrt{-g}\,T^{\mu\nu}h_\nu$ and is conserved if the non-fixed fields satisfy their Euler-Lagrange equations. The conservation of $\mathcal{B}^\mu$ can be seen as a consequence of an identity that is a generalization of $\nabla_\mu T^{\mu\nu}=0$ and is a consequence of the gauge symmetry of the Lagrangian. This identity holds in any configuration of the fixed fields if the non-fixed fields satisfy their Euler-Lagrange equations. We also show that $\mathcal{B}^\mu$ differs from the relevant canonical Noether current by the sum of an identically conserved current and a term that vanishes if the non-fixed fields are on-shell. As example we discuss the case of general, possibly fermionic, matter fields propagating in fixed gravitational and Yang-Mills background. We find that in this case the generalization of $\nabla_\mu T^{\mu\nu}=0$ is the Lorentz law $\nabla_\mu T^{\mu\nu} - F^{a\nu\lambda}\mathcal{J}_{a\lambda} = 0$, which holds as a consequence of the diffeomorphism, local Lorentz and Yang-Mills gauge symmetry of the matter Lagrangian. As a second simple example we consider the case of general fields propagating in a background that consists of a gravitational and a real scalar field.


title: Skew derivations on generalized Weyl algebras

abstract: A wide class of skew derivations on degree-one generalized Weyl algebras $R(a,\varphi)$ over a ring $R$ is constructed. All these derivations are twisted by a degree-counting extensions of automorphisms of $R$. It is determined which of the constructed derivations are $Q$-skew derivations. The compatibility of these skew derivations with the natural ${\mathbb Z}$-grading of $R(a,\varphi)$ is studied. Additional classes of skew derivations are constructed for generalized Weyl algebras given by an automorphism $\varphi$ of a finite order. Conditions that the central element $a$ that forms part of the structure of $R(a,\varphi)$ need to satisfy for the orthogonality of pairs of aforementioned skew derivations are derived. General constructions are illustrated by classification of skew derivations of generalized Weyl algebras over the polynomial ring in one variable and with a linear polynomial as the central element.


title: Non-integer characterizing slopes for torus knots

abstract: A slope $p/q$ is a characterizing slope for a knot $K$ in $S^3$ if the oriented homeomorphism type of $p/q$-surgery on $K$ determines $K$ uniquely. We show that for each torus knot its set of characterizing slopes contains all but finitely many non-integer slopes. This generalizes work of Ni and Zhang who established such a result for $T_{5,2}$. Along the way we show that if two knots $K$ and $K'$ in $S^3$ have homeomorphic $p/q$-surgeries, then for $q\geq 3$ and $p$ sufficiently large we can conclude that $K$ and $K'$ have the same genera and Alexander polynomials. This is achieved by consideration of the absolute grading on Heegaard Floer homology.


title: Small quenches and thermalization

abstract: We study the expectation values of observables and correlation functions at long times after a global quantum quench. Our focus is on metallic (`gapless') fermionic many-body models and small quenches. The system is prepared in an eigenstate of an initial Hamiltonian, and the time evolution is performed with a final Hamiltonian which differs from the initial one in the value of one global parameter. We first derive general relations between time-averaged expectation values of observables as well as correlation functions and those obtained in an eigenstate of the final Hamiltonian. Our results are valid to linear and quadratic order in the quench parameter g and generalize prior insights in several essential ways. This allows us to develop a phenomenology for the thermalization of local quantities up to a given order in g. Our phenomenology is put to a test in several case studies of one-dimensional models representative of four distinct classes of Hamiltonians: quadratic ones, effectively quadratic ones, those characterized by an extensive set of (quasi-) local integrals of motion, and those for which no such set is known (and believed to be nonexistent). We show that for each of these models, all observables and correlation functions thermalize to linear order in g. The more local a given quantity, the longer the linear behavior prevails when increasing g. Typical local correlation functions and observables for which the term O(g) vanishes thermalize even to order g^2. Our results show that lowest order thermalization of local observables is an ubiquitous phenomenon even in models with extensive sets of integrals of motion.


title: The Bramson logarithmic delay in the cane toads equations

abstract: We study a nonlocal reaction-diffusion-mutation equation modeling the spreading of a cane toads population structured by a phenotypical trait responsible for the spatial diffusion rate. When the trait space is bounded, the cane toads equation admits traveling wave solutions [7]. Here, we prove a Bramson type spreading result: the lag between the position of solutions with localized initial data and that of the traveling waves grows as (3/(2$\lambda$ *)) log t. This result relies on a present-time Harnack inequality which allows to compare solutions of the cane toads equation to those of a Fisher-KPP type equation that is local in the trait variable.


title: Nitrogen depletion in field red giants: mixing during the He flash?

abstract: We combine simultaneous constraints on stellar evolutionary status from asteroseismology, and on nitrogen abundances derived from large spectroscopic surveys, to follow nitrogen surface abundances all along the evolution of a low-mass star, comparing model expectations with data. After testing and calibrating the observed yields from the APOGEE survey, we first show that nitrogen surface abundances follow the expected trend after the first dredge-up occurred, i.e. that the more massive is the star the more nitrogen is enhanced. Moreover, the behaviour of nitrogen data along the evolution confirms the existence of non-canonical extra-mixing on the RGB for all low-mass stars in the field. But more surprisingly, the data indicate that nitrogen has been depleted between the RGB tip and the red clump. This may suggest that some nitrogen has been burnt near or at the He flash episode.


title: Inference for Empirical Wasserstein Distances on Finite Spaces

abstract: The Wasserstein distance is an attractive tool for data analysis but statistical inference is hindered by the lack of distributional limits. To overcome this obstacle, for probability measures supported on finitely many points, we derive the asymptotic distribution of empirical Wasserstein distances as the optimal value of a linear program with random objective function. This facilitates statistical inference (e.g. confidence intervals for sample based Wasserstein distances) in large generality. Our proof is based on directional Hadamard differentiability. Failure of the classical bootstrap and alternatives are discussed. The utility of the distributional results is illustrated on two data sets.


title: Embeddings and the (virtual) cohomological dimension of the braid and mapping class groups of surfaces

abstract: In this paper, we make use of the relations between the braid and mapping class groups of a compact, connected, non-orientable surface N without boundary and those of its orientable double covering S to study embeddings of these groups and their (virtual) cohomological dimensions. We first generalise results of Birman and Chillingworth and of Gon\c{c}alves and Guaschi to show that the mapping class group MCG(N ; k) of N relative to a k-point subset embeds in the mapping class group MCG(S; 2k) of S relative to a 2k-point subset. We then compute the cohomological dimension of the braid groups of all compact, connected aspherical surfaces without boundary. Finally, if the genus of N is greater than or equal to 2, we give upper bounds for the virtual cohomological dimension of MCG(N ; k).


title: Crunch-in regime - Non-linearly driven hollow-channel plasma

abstract: Plasma wakefields driven inside a hollow-channel plasma are significantly different from those driven in a homogeneous plasma. This work investigates the scaling laws of the accelerating and focusing fields in the "crunch-in" regime. This regime is excited due to the collapse of the electron-rings from the channel walls onto the propagation axis of the energy-source, in its wake. This regime is thus the non-linearly driven hollow channel, since the electron-ring displacement is of the order of the channel radius. We present the properties of the coherent structures in the "crunch-in" regime where the channel radius is matched to the beam properties such that channel-edge to on-axis collapse time has a direct correspondence to the energy source intensity. We also investigate the physical mechanisms that underlie the "crunch-in" wakefields by tuning the channel radius. Using a theoretical framework and results from PIC simulations the possible applications of the "crunch-in" regime for acceleration of positron beams with collider-scale parameters is presented.


title: From Kinetic Theory of Multicellular Systems to Hyperbolic Tissue Equations: Asymptotic Limits and Computing

abstract: This paper deals with the analysis of the asymptotic limit toward the derivation of macroscopic equations for a class of equations modeling complex multicellular systems by methods of the kinetic theory. After having chosen an appropriate scaling of time and space, a Chapman-Enskog expansion is combined with a closed, by minimization, technique to derive hyperbolic models at the macroscopic level. The resulting macroscopic equations show how the macroscopic tissue behavior can be described by hyperbolic systems which seem the most natural in this context. We propose also an asymptotic-preserving well-balanced scheme for the one-dimensional hyperbolic model, in the two dimensional case, we consider a time splitting method between the conservative part and the source term where the conservative equation is approximated by the Lax-Friedrichs scheme.


title: Learning an unknown transformation via a genetic approach

abstract: Recent developments in integrated photonics technology are opening the way to the fabrication of complex linear optical interferometers. The application of this platform is ubiquitous in quantum information science, from quantum simulation to quantum metrology, including the quest for quantum supremacy via the boson sampling problem. Within these contexts, the capability to learn efficiently the unitary operation of the implemented interferometers becomes a crucial requirement. In this letter we develop a reconstruction algorithm based on a genetic approach, which can be adopted as a tool to characterize an unknown linear optical network. We report an experimental test of the described method by performing the reconstruction of a 7-mode interferometer implemented via the femtosecond laser writing technique. Further applications of genetic approaches can be found in other contexts, such as quantum metrology or learning unknown general Hamiltonian evolutions.


title: Oscillating heat kernels on ultrametric spaces

abstract: Let $(X,d)$ be a proper ultrametric space. Given a measure $m$ on $X$ and a function $B \mapsto C(B)$ defined on the collection of all non-singleton balls $B$ of $X$, we consider the associated hierarchical Laplacian $L=L_{C}\,$. The operator $L$ acts in $\mathcal{L}^{2}(X,m),$ is essentially self-adjoint and has a pure point spectrum. It admits a continuous heat kernel $\mathfrak{p}(t,x,y)$ with respect to $m$. We consider the case when $X$ has a transitive group of isometries under which the operator $L$ is invariant and study the asymptotic behaviour of the function $t\mapsto \mathfrak{p}(t,x,x)=\mathfrak{p}(t)$. It is completely monotone, but does not vary regularly. When $X=\mathbb{Q}_{p}\,$, the ring of $p$-adic numbers, and $L=\mathcal{D}^{\alpha} $, the operator of \ fractional derivative of order $\alpha,$ we show that $\mathfrak{p}(t)=t^{-1/\alpha}\mathcal{A}% (\log_{p}t)$, where $\mathcal{A}(\tau)$ is a continuous non-constant $\alpha$-periodic function. We also study asymptotic behaviour of $\min\mathcal{A}$ and $\max\mathcal{A}$ as the space parameter $p$ tends to $\infty$. When $X=S_{\infty}\,$, the infinite symmetric group, and $L$ is a hierarchical Laplacian with metric structure analogous to $\mathcal{D}^{\alpha},$ we show that, contrary to the previous case, the completely monotone function $\mathfrak{p}(t)$ oscillates between two functions $\psi(t)$ and $\Psi(t)$ such that $\psi(t)/\Psi(t)\to 0$ as $t \to \infty\,$.


title: Few-particle systems: An analysis of some strongly correlated states

abstract: The analysis of the quantum Hall response of a small system of ultracold bosonic atoms through the variation of its Hall resistivity against the applied gauge magnetic field, provides a powerful method to unmask its strongly correlated states in a quite eshaustive way. Within a fixed range of values of the magnetic field in the lowest Landau-level regime, where the resistivity displays two successive plateaux, we identify the implied states as the Pfaffian and the state with filling factor \nu=2/3 in the thermodynamic limit. We fix the conditions to have good observability.


title: Observer's observables. Residual diffeomorphisms

abstract: We investigate the fate of diffeomorphisms when the radial gauge is imposed in canonical general relativity. As shown elsewhere, the radial gauge is closely related to the observer's observables. These observables are invariant under a large subgroup of diffeomorphisms which results in their usefulness for canonical general relativity. There are, however, some diffeomorphisms, called residual diffeomorphisms, which might be "observed" by the observer as they do not preserve her observables. The present paper is devoted to the analysis of these diffeomorphisms in the case of the spatial and spacetime radial gauges. Although the residual diffeomorphisms do not form a subgroup of all diffeomorphisms, we show that their induced action in the phase space does form a group. We find the generators of the induced transformations and compute the structure functions of the algebras they form. The obtained algebras are deformations of the algebra of the Euclidean group and the algebra of the Poincar\'e group in the spatial and spacetime case, respectively. In both cases the deformation depends only on the Riemann curvature tensor and in particular vanishes when the space or spacetime is flat.


title: Safe, Multi-Agent, Reinforcement Learning for Autonomous Driving

abstract: Autonomous driving is a multi-agent setting where the host vehicle must apply sophisticated negotiation skills with other road users when overtaking, giving way, merging, taking left and right turns and while pushing ahead in unstructured urban roadways. Since there are many possible scenarios, manually tackling all possible cases will likely yield a too simplistic policy. Moreover, one must balance between unexpected behavior of other drivers/pedestrians and at the same time not to be too defensive so that normal traffic flow is maintained. In this paper we apply deep reinforcement learning to the problem of forming long term driving strategies. We note that there are two major challenges that make autonomous driving different from other robotic tasks. First, is the necessity for ensuring functional safety - something that machine learning has difficulty with given that performance is optimized at the level of an expectation over many instances. Second, the Markov Decision Process model often used in robotics is problematic in our case because of unpredictable behavior of other agents in this multi-agent scenario. We make three contributions in our work. First, we show how policy gradient iterations can be used without Markovian assumptions. Second, we decompose the problem into a composition of a Policy for Desires (which is to be learned) and trajectory planning with hard constraints (which is not learned). The goal of Desires is to enable comfort of driving, while hard constraints guarantees the safety of driving. Third, we introduce a hierarchical temporal abstraction we call an "Option Graph" with a gating mechanism that significantly reduces the effective horizon and thereby reducing the variance of the gradient estimation even further.


title: A panorama of dispersion curves for the weighted isotropic relaxed micromorphic model

abstract: We consider the weighted isotropic relaxed micromorphic model and provide an in depth investigation of the characteristic dispersion curves when the constitutive parameters of the model are varied. The weighted relaxed micromorphic model generalizes the classical relaxed micromorphic model previously introduced by the authors, since it features the Cartan-Lie decomposition of the tensors $P_{,t}$ and Curl $P$ in their dev, dev sym, skew and spheric part. It is shown that the split of the tensor $P_{,t}$ in the micro-inertia provide an independent control of the cut-offs of the optic benches. This is crucial for the future calibration of the relaxed micromorphic model on real band-gap metamaterials. Even if the physical interest of the introduction of the split of the tensor Curl $P$ is less evident than in the previous case, we discuss in detail which is its effect on the dispersion curves. Finally, we also provide a complete parametric study involving all the constitutive parameters of the introduced model, so giving rise to an exhaustive panorama of dispersion curves for the relaxed micromorphic model.


title: Real submanifolds of maximum complex tangent space at a CR singular point, II

abstract: We study a germ of real analytic n-dimensional submanifold of $C^n$ that has a complex tangent space of maximal dimension at a CR singularity. Under the condition that its complexification admits the maximum number of deck transformations, we first classify holomorphically its quadratic CR singularity. We then study its transformation to a normal form under the action of local (possibly formal) biholomorphisms at the singularity. We first conjugate formally its associated reversible map $\sigma$ to suitable normal forms and show that all these normal forms can be divergent. We then construct a unique formal normal form under a non degeneracy condition.


title: Combining Treewidth and Backdoors for CSP

abstract: We show that CSP is fixed-parameter tractable when parameterized by the treewidth of a backdoor into any tractable CSP problem over a finite constraint language. This result combines the two prominent approaches for achieving tractability for CSP: (i) by structural restrictions on the interaction between the variables and the constraints and (ii) by language restrictions on the relations that can be used inside the constraints. Apart from defining the notion of backdoor-treewidth and showing how backdoors of small treewidth can be used to efficiently solve CSP, our main technical contribution is a fixed-parameter algorithm that finds a backdoor of small treewidth.


title: Small slot waveguide rings for on-chip quantum optical circuits

abstract: Nanophotonic interfaces between single emitters and light promise to enable new quantum optical technologies. Here, we use a combination of finite element simulations and analytic quantum theory to investigate the interaction of various quantum emitters with slot-waveguide rings. We predict that for rings with radii as small as 1.44 $\mu$m (Q = 27,900), near-unity emitter-waveguide coupling efficiencies and emission enhancements on the order of 1300 can be achieved. By tuning the ring geometry or introducing losses, we show that realistic emitter-ring systems can be made to be either weakly or strongly coupled, so that we can observe Rabi oscillations in the decay dynamics even for micron-sized rings. Moreover, we demonstrate that slot waveguide rings can be used to directionally couple emission, again with near-unity efficiency. Our results pave the way for integrated solid-state quantum circuits involving various emitters.


title: Stability, convergence to equilibrium and simulation of non-linear Hawkes Processes with memory kernels given by the sum of Erlang kernels

abstract: Non-linear Hawkes processes with memory kernels given by the sum of Erlang kernels are considered. It is shown that their stability properties can be studied in terms of an associated class of piecewise deterministic Markov processes, called Markovian cascades of successive memory terms. Explicit conditions implying the positive Harris recurrence of these processes are presented. The proof is based on integration by parts with respect to the jump times. A crucial property is the non-degeneracy of the transition semigroup which is obtained thanks to the invertibility of an associated Vandermonde matrix. For Lipschitz continuous rate functions we also show that these Markovian cascades converge to equilibrium exponentially fast with respect to the Wasserstein distance. Finally, an extension of the classical thinning algorithm is proposed to simulate such Markovian cascades.


title: An analogue to Dixon's theorem for automaton groups

abstract: Dixon's famous theorem states that the group generated by two random permutations of a finite set is generically either the whole symmetric group or the alternating group. In the context of random generation of finite groups this means that it is hopeless to wish for a uniform distribution -- or even a non-trivial one -- by drawing random permutations and looking at the generated group. Mealy automata are a powerful tool to generate groups, including all finite groups and many interesting infinite ones, whence the idea of generating random finite groups by drawing random Mealy automata. In this paper we show that, for a special class of Mealy automata that generate only finite groups, the distribution is far from being uniform since the obtained groups are generically a semi-direct product between a direct product of alternating groups and a group generated by a tuple of transpositions.


title: The universal C*-algebra of the electromagnetic field II. Topological charges and spacelike linear fields

abstract: Conditions for the appearance of topological charges are studied in the framework of the universal C*-algebra of the electromagnetic field, which is represented in any theory describing electromagnetism. It is shown that non-trivial topological charges, described by pairs of fields localised in certain topologically non-trivial spacelike separated regions, can appear in regular representations of the algebra only if the fields depend non-linearly on the mollifying test functions. On the other hand, examples of regular vacuum representations with non-trivial topological charges are constructed, where the underlying field still satisfies a weakened form of "spacelike linearity". Such representations also appear in the presence of electric currents. The status of topological charges in theories with several types of electromagnetic fields, which appear in the short distance (scaling) limit of asymptotically free non-abelian gauge theories, is also briefly discussed.


title: The Phase Transition in 5 Point Energy Minimization

abstract: Let R_s(r)=sign(s)/r^s be the Riesz s-energy potential. (This is the usual power-law potential.) This monograph proves the existence of a computable number S=15.048... such that the triangular bi-pyramid is the unique minimizer with respect to R_s, amongst all 5-point configurations on the sphere, if and only if s lies in (-2,0) or (0,S). This establishes the existence of the long-conjectured phase transition constant in 5-point energy minimization.


title: Inflation via logarithmic entropy-corrected holographic dark energy model

abstract: We study the inflation via logarithmic entropy-corrected holographic dark energy LECHDE model with future event horizon, particle horizon and Hubble horizon cut-offs, and compare the results with those of obtained in the study of inflation by holographic dark energy HDE model. In comparison, the spectrum of primordial scalar power spectrum in the LECHDE model becomes redder than the spectrum in HDE model. Moreover, the consistency with the observational data in LECHDE model of inflation, constrains the reheating temperature and Hubble parameter by one parameter of holographic dark energy and two new parameters of logarithmic corrections.


title: Model Based In Situ Calibration of Six Axis Force Torque Sensors

abstract: This paper proposes and validates an in situ calibration method to calibrate six axis force torque (F/T) sensors once they are mounted on the system. This procedure takes advantage of the knowledge of the model of the robot to generate the expected wrenches of the sensors during some arbitrary motions. It then uses this information to train and validate new calibration matrices, taking into account the calibration matrix obtained with a classical Workbench calibration. The proposed calibration algorithm is validated on the F/T sensors mounted on the iCub humanoid robot legs.


title: Projective dimension and regularity of path ideals of cycles

abstract: In this paper, we give a formula to compute all the top degree graded Betti numbers of the path ideals of a cycle. As a consequence we can give a formula to compute its projective dimension and regularity.


title: Ball Intersection Properties in Metric Spaces

abstract: We show that in complete metric spaces, $4$-hyperconvexity is equivalent to finite hyperconvexity. Moreover, every complete, almost $n$-hyperconvex metric space is $n$-hyperconvex. This generalizes among others results of Lindenstrauss and answers questions of Aronszajn-Panitchpakdi. Furthermore, we prove local-to-global results for externally and weakly externally hyperconvex subsets of hyperconvex metric spaces and find sufficient conditions in order for those classes of subsets to be convex with respect to a geodesic bicombing.


title: Exact ensemble density functional theory for excited states in a model system: investigating the weight dependence of the correlation energy

abstract: Ensemble density functional theory (eDFT) is an exact time-independent alternative to time-dependent DFT (TD-DFT) for the calculation of excitation energies. Despite its formal simplicity and advantages in contrast to TD-DFT (multiple excitations, for example, can be easily taken into account in an ensemble), eDFT is not standard which is essentially due to the lack of reliable approximate exchange-correlation (xc) functionals for ensembles. Following Smith et al. [Phys. Rev. B 93, 245131 (2016)], we propose in this work to construct an exact eDFT for the nontrivial asymmetric Hubbard dimer, thus providing more insight into the weight dependence of the ensemble xc energy in various correlation regimes. For that purpose, an exact analytical expression for the weight-dependent ensemble exchange energy has been derived. The complementary exact ensemble correlation energy has been computed by means of Legendre--Fenchel transforms. Interesting features like discontinuities in the ensemble xc potential in the strongly correlated limit have been rationalized by means of a generalized adiabatic connection formalism. Finally, functional-driven errors induced by ground-state density-functional approximations have been studied. In the strictly symmetric case or in the weakly correlated regime, combining ensemble exact exchange with ground-state correlation functionals gives relatively accurate ensemble energies. However, when approaching the equiensemble in the strongly correlated regime, this approximation leads to highly curved ensemble energies with negative slope which is unphysical. Using both ground-state exchange and correlation functionals gives much better results in that case. In fact, exact ensemble energies are almost recovered in some density domains. The analysis of density-driven errors is left for future work.


title: A Low-Overhead Hybrid Approach for Universal Fault-Tolerant Quantum Computation

abstract: As there is no quantum error correction code with universal set of transversal gates, several approaches have been proposed which, in combination of transversal gates, make universal fault-tolerant quantum computation possible. Magic state distillation, code switching, code concatenation and pieceable fault-tolerance are well-known examples of such approaches. However, the overhead of these approaches is one of the main bottlenecks for large-scale quantum computation. In this paper, a hybrid approach is proposed which combines the code concatenation technique with the other mentioned approaches. The proposed approach outperforms code concatenation in terms of both number of qubits and error threshold and also significantly reduces the resource overhead of code switching, magic state distillation and pieceable fault-tolerance at the cost of reducing the effective distance of the concatenated code for implementing non-transversal gates.


title: The Relativistic Hamilton-Jacobi Equation for a Massive, Charged and Spinning Particle, its Equivalent Dirac Equation and the de Broglie-Bohm Theory

abstract: Using the Clifford and the Spin-Clifford formalisms we prove that the classical relativistic Hamilton Jacobi equation for a charged massive (and spinning) particle interacting with an external electromagnetic field is equivalent to a Dirac-Hestenes equation satisfied by a class of spinor fields that we call classical spinor fields, characterized for having the Takabayashi angle function constant (equal to 0 or {\pi}). We also investigate a nonlinear Dirac-Hestenes like equation that comes from some class of generalized classical spinor fields. Finally we show that the general Dirac-Hestenes equation (which is a representative in the Clifford bundle of the usual Dirac equation) gives a generalized Hamilton-Jacobi equation where the quantum potential satisfy a severe constraint and the "mass of the particle" becomes a variable. Our results can then explain the experimental discrepancies found between prediction of the de Broglie-Bohm theory and recent experiments. We briefly discuss also the de Broglie's double solution theory in view of our results showing that it can be realized, at least in the case of spinning free particles..The paper contains several Appendices where notation and proofs of some results of the text are presented.


title: Characterizing the gamma-ray long-term variability of PKS 2155-304 with H.E.S.S. and Fermi-LAT

abstract: Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain. Over the course of ~9 yr of H.E.S.S observations the VHE light curve in the quiescent state is consistent with a log-normal behavior. The VHE variability in this state is well described by flicker noise (power-spectral-density index {\ss}_VHE = 1.10 +0.10 -0.13) on time scales larger than one day. An analysis of 5.5 yr of HE Fermi LAT data gives consistent results ({\ss}_HE = 1.20 +0.21 -0.23, on time scales larger than 10 days) compatible with the VHE findings. The HE and VHE power spectral densities show a scale invariance across the probed time ranges. A direct linear correlation between the VHE and HE fluxes could neither be excluded nor firmly established. These long-term-variability properties are discussed and compared to the red noise behavior ({\ss} ~ 2) seen on shorter time scales during VHE-flaring states. The difference in power spectral noise behavior at VHE energies during quiescent and flaring states provides evidence that these states are influenced by different physical processes, while the compatibility of the HE and VHE long-term results is suggestive of a common physical link as it might be introduced by an underlying jet-disk connection.


title: Nuclear magnetic resonance study of the magnetic-field-induced ordered phase in the NiCl2-4SC(NH2)2 compound

abstract: Nuclear magnetic resonance (NMR) study of the high magnetic field (H) part of the Bose-Einstein condensed (BEC) phase of the quasi-onedimensional (quasi-1D) antiferromagnetic quantum spin-chain compound NiCl2-4SC(NH2)2 (DTN) was performed. We precisely determined the phase boundary, Tc(H), down to 40 mK; the critical boson density, n_c(Tc); and the absolute value of the BEC order parameter S_perp at very low temperature (T = 0.12 K). All results are accurately reproduced by numerical quantum Monte Carlo simulations of a realistic three-dimensional (3D) model Hamiltonian. Approximate analytical predictions based on the 1D Tomonaga-Luttinger liquid description are found to be precise for Tc(H), but less so for S_perp(H), which is more sensitive to the strength of 3D couplings, in particular close to the critical field. A mean-field treatment, based on the Hartree-Fock-Popov description, is found to be valid only up to n_c = 4% (T < 0.3 K), while for higher n_c boson interactions appear to modify the density of states.


title: Stability Limits and Surface Chemistry of Ag Nanoparticles in Non-Adsorbing Electrolytes Probed by Bragg Coherent Diffractive Imaging

abstract: Surface chemistry is important across diverse fields such as corrosion and nanostructure synthesis. Unfortunately, many as-synthesized nanomaterials, including partially dealloyed nanoparticle catalysts for fuel cells, with highly active surfaces are not stable in their reactive environments, preventing widespread application. Thus, understanding instability by focusing on the structure-stability and defect-stability relationship at the nanoscale is crucial and will likely play an important role in meeting grand challenges. To this end, recent advances in imaging nanostructure stability have come via both electron, x-ray, and other techniques such as atomic force microscopy, but tend to be limited to specific sample environments and/or two-dimensional images. Here, we report investigations into the defect-stability relationship of silver nanoparticles to voltage-induced electrochemical dissolution imaged in-situ in three-dimensional (3D) detail by Bragg Coherent Diffractive Imaging (BCDI). We first determine the average dissolution kinetics by Stationary Probe Rotating Disk Electrode (SPRDE) in combination with inductively coupled plasma mass spectrometry (ICP-MS), which allows real-time in-situ measurement of Ag+ ions formation and the corresponding electrochemical current. We then observe the dissolution and redeposition processes in 3D with BCDI in single nanocrystals, providing unique insight about the role of surface strain, defects, and their coupling to the dissolution chemistry. The methods developed and the knowledge gained go well beyond a "simple" silver electrochemistry and are applicable to all electrocatalytic reactions where functional links between activity and stability are controlled by structure and defect dynamics.


title: Spontaneous particle-hole symmetry breaking of correlated fermions on the Lieb lattice

abstract: We study spinless fermions with nearest-neighbor repulsive interactions ($t$-$V$ model) on the two-dimensional three-band Lieb lattice. At half-filling, the free electronic band structure consists of a flat band at zero energy and a single cone with linear dispersion. The flat band is expected to be unstable upon inclusion of electronic correlations, and a natural channel is charge order. However, due to the three-orbital unit cell, commensurate charge order implies an imbalance of electron and hole densities and therefore doping away from half-filling. Our numerical results show that below a finite-temperature Ising transition a charge density wave with one electron and two holes per unit cell and its partner under particle-hole transformation are spontaneously generated. Our calculations are based on recent advances in auxiliary-field and continuous-time quantum Monte Carlo simulations that allow sign-free simulations of spinless fermions at half-filling. It is argued that particle-hole symmetry breaking provides a route to access levels of finite doping, without introducing a sign problem.


title: Lorentz-violating neutral-pion decays in isotropic modified Maxwell theory

abstract: We consider an extension of the Standard Model with isotropic nonbirefringent Lorentz violation in the photon sector and restrict the discussion to the case of a "fast" photon with a phase velocity larger than the maximum attainable velocity of the quarks and leptons. With our conventions, this case corresponds to a negative Lorentz-violating parameter $\kappa$ in the action. The decay rate of a neutral pion into two nonstandard photons is calculated as a function of the 3-momentum of the initial pion and the negative Lorentz-violating parameter $\kappa$ of the final photons.


title: Learning from Survey Training Samples: Rate Bounds for Horvitz-Thompson Risk Minimizers

abstract: The generalization ability of minimizers of the empirical risk in the context of binary classification has been investigated under a wide variety of complexity assumptions for the collection of classifiers over which optimization is performed. In contrast, the vast majority of the works dedicated to this issue stipulate that the training dataset used to compute the empirical risk functional is composed of i.i.d. observations. Beyond the cases where training data are drawn uniformly without replacement among a large i.i.d. sample or modelled as a realization of a weakly dependent sequence of r.v.'s, statistical guarantees when the data used to train a classifier are drawn by means of a more general sampling/survey scheme and exhibit a complex dependence structure have not been documented yet. It is the main purpose of this paper to show that the theory of empirical risk minimization can be extended to situations where statistical learning is based on survey samples and knowledge of the related inclusion probabilities. Precisely, we prove that minimizing a weighted version of the empirical risk, refered to as the Horvitz-Thompson risk (HT risk), over a class of controlled complexity lead to a rate for the excess risk of the order $O_{\mathbb{P}}((\kappa_N (\log N)/n)^{1/2})$ with $\kappa_N=(n/N)/\min_{i\leq N}\pi_i$, when data are sampled by means of a rejective scheme of (deterministic) size $n$ within a statistical population of cardinality $N\geq n$, a generalization of basic {\it sampling without replacement} with unequal probability weights $\pi_i>0$. Extension to other sampling schemes are then established by a coupling argument. Beyond theoretical results, numerical experiments are displayed in order to show the relevance of HT risk minimization and that ignoring the sampling scheme used to generate the training dataset may completely jeopardize the learning procedure.


title: A Greedy Approach for Budgeted Maximum Inner Product Search

abstract: Maximum Inner Product Search (MIPS) is an important task in many machine learning applications such as the prediction phase of a low-rank matrix factorization model for a recommender system. There have been some works on how to perform MIPS in sub-linear time recently. However, most of them do not have the flexibility to control the trade-off between search efficient and search quality. In this paper, we study the MIPS problem with a computational budget. By carefully studying the problem structure of MIPS, we develop a novel Greedy-MIPS algorithm, which can handle budgeted MIPS by design. While simple and intuitive, Greedy-MIPS yields surprisingly superior performance compared to state-of-the-art approaches. As a specific example, on a candidate set containing half a million vectors of dimension 200, Greedy-MIPS runs 200x faster than the naive approach while yielding search results with the top-5 precision greater than 75\%.


title: Determination of hole g-factor in InAs/InGaAs/InAlAs quantum wells by magneto-photoluminescence studies

abstract: Circularly-polarized magneto-photoluminescence (magneto-PL) technique has been applied to investigate Zeeman effect in InAs/InGaAs/InAlAs quantum wells (QWs) in Faraday geometry. Structures with different thickness of the QW barriers have been studied in magnetic field parallel and tilted with respect to the sample normal. Effective electron-hole g-factor has been found by measurement of splitting of polarized magneto-PL lines. Lande factors of electrons have been calculated using the 14-band kp method and g-factor of holes was determined by subtracting the calculated contribution of the electrons from the effective electron-hole g-factor. Anisotropy of the hole g-factor has been studied applying tilted magnetic field.


title: A locally p-adaptive approach for Large Eddy Simulation of compressible flows in a DG framework

abstract: We investigate the possibility of reducing the computational burden of LES models by employing local polynomial degree adaptivity in the framework of a high order DG method. A novel degree adaptation technique especially featured to be effective for LES applications is proposed and its effectiveness is compared to that of other criteria already employed in the literature. The resulting locally adaptive approach allows to achieve significant reductions in computational cost of representative LES computations.


title: The Galactic O-Star Catalog (GOSC) and the Galactic O-Star Spectroscopic Survey (GOSSS): current status

abstract: We present the updates of the Galactic O-Star Catalog (GOSC) that we have undertaken in the last two years: new spectral types, more objects, additional information, and coordination with CDS. We also present updates for the Galactic O-Star Spectroscopic Survey (GOSSS). A new paper (GOSSS-III) has been published and ~1000 targets have been observed since 2014. Four new setups have been added to our lineup and for two of them we have already obtained over 100 spectra: with OSIRIS at the 10.4 m GTC we are observing northern dim stars and with FRODOspec at the 2.0 m Liverpool Telescope we are observing northern bright stars. Finally, we also make available new versions of MGB, the spectral classification tool associated with the project, and of the GOSSS grid of spectroscopic standards.


title: Keystroke dynamics as signal for shallow syntactic parsing

abstract: Keystroke dynamics have been extensively used in psycholinguistic and writing research to gain insights into cognitive processing. But do keystroke logs contain actual signal that can be used to learn better natural language processing models? We postulate that keystroke dynamics contain information about syntactic structure that can inform shallow syntactic parsing. To test this hypothesis, we explore labels derived from keystroke logs as auxiliary task in a multi-task bidirectional Long Short-Term Memory (bi-LSTM). Our results show promising results on two shallow syntactic parsing tasks, chunking and CCG supertagging. Our model is simple, has the advantage that data can come from distinct sources, and produces models that are significantly better than models trained on the text annotations alone.


title: Optimal parameters uncoupling vibration modes of oscillators

abstract: A novel optimization concept for an oscillator with two degrees of freedom is proposed. By using specially defined motion ratios, we control the action of springs and dampers to each degree of freedom of the oscillator. If the potential action of the springs in one period of vibration, used as the payoff function for the conservative oscillator, is maximized, then the optimal motion ratios uncouple vibration modes. The same result holds true for the dissipative oscillator. The application to optimal design of vehicle suspension is discussed.


title: How the nonlinear coupled oscillators modelization explains the Blazhko effect, the synchronisation of layers, the mode selection, the limit cycle, and the red limit of the instability strip

abstract: Context. The Blazhko effect, in RR Lyrae type stars, is a century old mystery. Dozens of theory exists, but none have been able to entirely reproduce the observational facts associated to this modulation phenomenon. Existing theory all rely on the usual continuous modelization of the star. Aims. We present a new paradigm which will not only explain the Blazhko effect, but at the same time, will give us alternative explanations to the red limit of the instability strip, the synchronization of layers, the mode selection and the existence of a limit cycle for radially pulsating stars. Methods. We describe the RR Lyrae type pulsating stars as a system of coupled nonlinear oscillators. Considering a spatial discretisation of the star, supposing a spherical symmetry, we develop the equation of motion and energy up to the third order in the radial and adiabatic case. Then, we include the influence of the ionization region as a relaxation oscillator by including elements from synchronisation theory. Results. This discrete approach allows us to exploit existing results in the coupled nonlinear oscillator field. For instance, the study of synchronicity leads to an explanation of the mode selection, the layers synchronisation, the limit cycle and the red limit of the instability strip. But, most of all, the analogy with the Fermi-Pasta-Ulam (FPU) experiment enables us to understand the Blazhko effect. The transfer of energy between different modes, as induced by solitons, not only gives a plausible theory for lightcurve modulation, but also explains the asymmetry of sidelobes.


title: A grid of 1D low-mass star formation collapse models

abstract: The current study was developed to provide a database of relatively simple numerical simulations of protostellar collapse, as a template library for observations of cores and very young protostars, and for researchers who wish to test their chemical modeling under dynamic astrophysical conditions. It was also designed to identify statistical trends that may appear when running many models of the formation of low-mass stars by varying the initial conditions. A large set of 143 calculations of the gravitational collapse of an isolated sphere of gas with uniform temperature and a Bonnor-Ebert like density profile was undertaken using a 1D fully implicit Lagrangian radiation hydrodynamics code. The parameter space covered initial masses from 0.2 to 8 Msun, temperatures of 5-30 K and radii between 3000 and 30,000 AU. A spread in the thermal evolutionary tracks of the runs was found, due to differing initial conditions and optical depths. Within less than an order of magnitude, all first and second Larson cores had masses and radii independent of the initial conditions. The time elapsed between the formation of the first and second cores was found to strongly depend on the first core mass accretion rate, and no first core in our grid of models lived for longer than 2000 years, before the onset of the second collapse. The end product of a protostellar cloud collapse, the second Larson core, is, at birth, a canonical object with a mass and radius of about 3 Mjup and 8 Rjup, independent of its initial conditions. The evolution sequence which brings the gas to stellar densities can however proceed in a variety of scenarios, on different timescales, along different isentropes, but each story line can largely be predicted by the initial conditions. All the data from the simulations are publicly available at this address: http://starformation.hpc.ku.dk/grid-of-protostars.


title: High temperature condensate clouds in super-hot Jupiter atmospheres

abstract: Deciphering the role of clouds is central to our understanding of exoplanet atmospheres, as they have a direct impact on the temperature and pressure structure, and observational properties of the planet. Super-hot Jupiters occupy a temperature regime similar to low mass M-dwarfs, where minimal cloud condensation is expected. However, observations of exoplanets such as WASP-12b (Teq ~ 2500 K) result in a transmission spectrum indicative of a cloudy atmosphere. We re-examine the temperature and pressure space occupied by these super-hot Jupiter atmospheres, to explore the role of the initial Al- and Ti-bearing condensates as the main source of cloud material. Due to the high temperatures a majority of the more common refractory material is not depleted into deeper layers and would remain in the vapor phase. The lack of depletion into deeper layers means that these materials with relatively low cloud masses can become significant absorbers in the upper atmosphere. We provide condensation curves for the initial Al- and Ti-bearing condensates that may be used to provide quantitative estimates of the effect of metallicity on cloud masses, as planets with metal-rich hosts potentially form more opaque clouds because more mass is available for condensation. Increased metallicity also pushes the point of condensation to hotter, deeper layers in the planetary atmosphere further increasing the density of the cloud. We suggest that planets around metal-rich hosts are more likely to have thick refractory clouds, and discuss the implication on the observed spectra of WASP-12b.


title: Stochastic continuity equation with non-smooth velocity

abstract: In this article we study the existence and uniqueness of solutions of stochastic continuity equation with irregular coefficients.