title: Coupling and competition between ferroelectricity, magnetism, strain and oxygen vacancies in AMnO3 perovskites

abstract: We use first-principles calculations based on density functional theory to investigate the interplay between oxygen vacancies, A-site cation size / tolerance factor, epitaxial strain, ferroelectricity and magnetism in the perovskite manganite series, AMnO3 (A=Ca2+, Sr2+, Ba2+). We find that, as expected, increasing the volume through either chemical pressure or tensile strain generally lowers the formation energy of neutral oxygen vacancies consistent with their established tendency to expand the lattice. Increased volume also favors polar distortions, both because competing rotations of the oxygen octahedra are suppressed and because coulomb repulsion associated with cation off-centering is reduced. Interestingly, the presence of ferroelectric polarization favors ferromagnetic over antiferromagnetic ordering due to suppressed antiferromagnetic superexchange as the polar distortion bends the Mn-O-Mn bond angles away from the optimal 180deg. Intriguingly, we find that polar distortions compete with the formation of oxygen vacancies, which have a higher formation energy in the polar phases; conversely the presence of oxygen vacancies suppresses the onset of polarization. In contrast, oxygen vacancy formation energies are lower for ferromagnetic than antiferromagnetic orderings of the same structure type. Our findings suggest a rich and complex phase diagram, in which defect chemistry, polarization, structure and magnetism can be modified using chemical potential, strain, and electric or magnetic fields.


title: A Burau-Alexander 2-functor on tangles

abstract: We construct a weak 2-functor from the bicategory of oriented tangles to a bicategory of Lagrangian cospans. This functor simultaneously extends the Burau representation of the braid groups, its generalization to tangles due to Turaev and the first-named author, and the Alexander module of 1 and 2-dimensional links.


title: On goodness-of-fit tests for parametric hypotheses in perturbed dynamical systems using a minimum distance estimator

abstract: We consider the problem of the construction of the Goodness-of-Fit test in the case of continuous time observations of a diffusion process with small noise. The null hypothesis is parametric and we use a minimum distance estimator of the unknown parameter. We propose an asymptotically distribution free test for this model.


title: Energy Conservation and Coupling Error Reduction in Non-Iterative Co-Simulations

abstract: When simulators are energetically coupled in a co-simulation, residual energies alter the total energy of the full coupled system. This distorts the system dynamics, lowers the quality of the results, and can lead to instability. By using power bonds to realize simulator coupling, the Energy-Conservation-based Co-Simulation method (ECCO) [Sadjina et al. 2016] exploits these concepts to define non-iterative global error estimation and adaptive step size control relying on coupling variable data alone. Following similar argumentation, the Nearly Energy Preserving Coupling Element (NEPCE) [Benedikt et al. 2013] uses corrections to the simulator inputs to approximately ensure energy conservation. Here, we discuss a modification to NEPCE for when direct feed-through is present in one of the coupled simulators. We further demonstrate how accuracy and efficiency in non-iterative co-simulations are substantially enhanced when combining NEPCE with ECCO's adaptive step size controller. A quarter car model with linear and nonlinear damping characteristics serves as a co-simulation benchmark, and we observe reductions of the coupling errors of up to 98% utilizing the concepts discussed here.


title: Learning from Non-Stationary Stream Data in Multiobjective Evolutionary Algorithm

abstract: Evolutionary algorithms (EAs) have been well acknowledged as a promising paradigm for solving optimisation problems with multiple conflicting objectives in the sense that they are able to locate a set of diverse approximations of Pareto optimal solutions in a single run. EAs drive the search for approximated solutions through maintaining a diverse population of solutions and by recombining promising solutions selected from the population. Combining machine learning techniques has shown great potentials since the intrinsic structure of the Pareto optimal solutions of an multiobjective optimisation problem can be learned and used to guide for effective recombination. However, existing multiobjective EAs (MOEAs) based on structure learning spend too much computational resources on learning. To address this problem, we propose to use an online learning scheme. Based on the fact that offsprings along evolution are streamy, dependent and non-stationary (which implies that the intrinsic structure, if any, is temporal and scale-variant), an online agglomerative clustering algorithm is applied to adaptively discover the intrinsic structure of the Pareto optimal solution set; and to guide effective offspring recombination. Experimental results have shown significant improvement over five state-of-the-art MOEAs on a set of well-known benchmark problems with complicated Pareto sets and complex Pareto fronts.


title: Non-convex balls in the Teichm\"uller metric

abstract: We prove that the Teichm\"uller space of surfaces of genus $\mathbf{g}$ with $\mathbf{p}$ punctures contains balls which are not convex in the Teichm\"uller metric whenever $3\mathbf{g}-3+\mathbf{p} > 1$.


title: Realistic model for a fifth force explaining anomaly in ${^8Be^*} \to {^8Be} \;{e^+e^-}$ Decay

abstract: A $6.8\,\sigma$ anomaly has been reported in the opening angle and invariant mass distributions of $e^+e^-$ pairs produced in ${^8Be}$ nuclear transitions. It has been shown that a protophobic fifth force mediated by a $17\,\textrm{MeV}$ gauge boson $X$ with pure vector current interactions can explain the data through the decay of an excited state to the ground state, ${^8Be^*} \to {^8Be}\, X$, and then the followed saturating decay $X \to e^+e^-$. In this work we propose a renormalizable model to realize this fifth force. Although axial-vector current interactions also exist in our model, their contributions cancel out in the iso-scalar interaction for ${^8Be^*} \to {^8Be} \,X$. Within the allowed parameter space, this model can alleviate the $(g-2)_\mu$ anomaly problem and can be probed by the LHCb experiment. Several other implications are discussed.


title: Asynchronous simulation of Boolean networks by monotone Boolean networks

abstract: We prove that the fully asynchronous dynamics of a Boolean network $f:\{0,1\}^n\to\{0,1\}^n$ without negative loop can be simulated, in a very specific way, by a monotone Boolean network with $2n$ components. We then use this result to prove that, for every even $n$, there exists a monotone Boolean network $f:\{0,1\}^n\to\{0,1\}^n$, an initial configuration $x$ and a fixed point $y$ of $f$ such that: (i) $y$ can be reached from $x$ with a fully asynchronous updating strategy, and (ii) all such strategies contains at least $2^{\frac{n}{2}}$ updates. This contrasts with the following known property: if $f:\{0,1\}^n\to\{0,1\}^n$ is monotone, then, for every initial configuration $x$, there exists a fixed point $y$ such that $y$ can be reached from $x$ with a fully asynchronous strategy that contains at most $n$ updates.


title: Partial $W^{2,p}$ regularity for optimal transport maps

abstract: We prove that, in the optimal transportation problem with general costs and positive continuous densities, the potential function is always of class $W^{2,p}_{loc}$ for any $p \geq 1$ outside of a closed singular set of measure zero. We also establish global $W^{2,p}$ estimates when the cost is a small perturbation of the quadratic cost. The latter result is new even when the cost is exactly the quadratic cost.


title: Deep Reinforcement Learning Discovers Internal Models

abstract: Deep Reinforcement Learning (DRL) is a trending field of research, showing great promise in challenging problems such as playing Atari, solving Go and controlling robots. While DRL agents perform well in practice we are still lacking the tools to analayze their performance. In this work we present the Semi-Aggregated MDP (SAMDP) model. A model best suited to describe policies exhibiting both spatial and temporal hierarchies. We describe its advantages for analyzing trained policies over other modeling approaches, and show that under the right state representation, like that of DQN agents, SAMDP can help to identify skills. We detail the automatic process of creating it from recorded trajectories, up to presenting it on t-SNE maps. We explain how to evaluate its fitness and show surprising results indicating high compatibility with the policy at hand. We conclude by showing how using the SAMDP model, an extra performance gain can be squeezed from the agent.


title: Haldane linearisation done right: Solving the nonlinear recombination equation the easy way

abstract: The nonlinear recombination equation from population genetics has a long history and is notoriously difficult to solve, both in continuous and in discrete time. This is particularly so if one aims at full generality, thus also including degenerate parameter cases. Due to recent progress for the continuous time case via the identification of an underlying stochastic fragmentation process, it became clear that a direct general solution at the level of the corresponding ODE itself should also be possible. This paper shows how to do it, and how to extend the approach to the discrete-time case as well.


title: Contact matrix in dilute quantum systems

abstract: Contact has been well established as an important quantity to govern dilute quantum systems, in which the pairwise correlation at short distance traces a broad range of thermodynamic properties. So far, studies have been focusing on contact in individual angular momentum channels. Here, we point out that, to have a complete description of the pairwise correlation in a general dilute quantum systems, contact should be defined as a matrix. Whereas the diagonal terms of such matrix include contact of all partial wave scatterings, the off-diagonal terms, which elude previous studies in the literature, characterise the coherence of the asymptotic pairwise wavefunction in the angular momentum space and determine important thermodynamic quantities including the momentum distribution. Contact matrix allows physicists to access unexplored connections between short-range correlations and macroscopic quantum phenomena. As an example, we show the direct connection between contact matrix and order parameters of a superfluid with mixed partial waves.


title: HARQ and AMC: Friends or Foes?

abstract: To ensure reliable communication in randomly varying and error-prone channels, wireless systems use adaptive modulation and coding (AMC) as well as hybrid ARQ (HARQ). In order to elucidate their compatibility and interaction, we compare the throughput provided by AMC, HARQ, and their combination (AMC-HARQ) under two operational conditions: in slow- and fast block-fading channels. Considering both, incremental redundancy HARQ (HARQ-IR) and repetition redundancy HARQ (HARQ-RR) we optimize the rate-decision regions for AMC/HARQ and compare them in terms of attainable throughput. Under a fairly general model of the channel variation and the decoding functions, we conclude that i) adding HARQ on top of AMC may be counterproductive in the high average signal- to-noise ratio regime for fast fading channels, and ii) HARQ is useful for slow fading channels, but it provides moderate throughput gains. We provide explanations for these results which allow us to propose paths to improve AMC-HARQ systems.


title: Coexistence of Weyl Physics and Planar Defects in Semimetals TaP and TaAs

abstract: We report a structural study of the Weyl semimetals TaAs and TaP, utilizing diffraction and imaging techniques, where we show that they contain a high density of defects, leading to non-stoichiometric single crystals of both semimetals. Despite the observed defects and non-stoichiometry on samples grown using techniques already reported in the literature, de Haas-van Alphen measurements on TaP reveal quantum oscillations and a high carrier mobility, an indication that the crystals are of quality comparable to those reported elsewhere. Electronic structure calculations on TaAs reveal that the position of the Weyl points relative to the Fermi level shift with the introduction of vacancies and stacking faults. In the case of vacancies the Fermi surface becomes considerably altered, while the effect of stacking faults on the electronic structure is to allow the Weyl pockets to remain close to the Fermi surface. The observation of quantum oscillations in a non-stoichiometric crystal and the persistence of Weyl fermion pockets near the Fermi surface in a crystal with stacking faults point to the robustness of these quantum phenomena in these materials.


title: Incentivized Campaigning in Social Networks

abstract: Campaigners, advertisers and activists are increasingly turning to social recommendation mechanisms, provided by social media, for promoting their products, services, brands and even ideas. However, many times, such social network based campaigns perform poorly in practice because the intensity of the recommendations drastically reduces beyond a few hops from the source. A natural strategy for maintaining the intensity is to provide incentives. In this paper, we address the problem of minimizing the cost incurred by the campaigner for incentivizing a fraction of individuals in the social network, while ensuring that the campaign message reaches a given expected fraction of individuals. We also address the dual problem of maximizing the campaign penetration for a resource constrained campaigner. To help us understand and solve the above mentioned problems, we use percolation theory to formally state them as optimization problems. These problems are not amenable to traditional approaches because of a fixed point equation that needs to be solved numerically. However, we use results from reliability theory to establish some key properties of the fixed point, which in turn enables us to solve these problems using algorithms that are linearithmic in maximum node degree. Furthermore, we evaluate the efficacy of the analytical solution by performing simulations on real world networks.


title: Partial redistribution in 3D non-LTE radiative transfer in solar atmosphere models

abstract: Resonance spectral lines such as H I Ly {\alpha}, Mg II h&k, and Ca II H&K that form in the solar chromosphere are influenced by the effects of 3D radiative transfer as well as partial redistribution (PRD). So far no one has modeled these lines including both effects simultaneously owing to the high computing demands of existing algorithms. Such modeling is however indispensable for accurate diagnostics of the chromosphere. We present a computationally tractable method to treat PRD scattering in 3D model atmospheres using a 3D non-LTE radiative transfer code. To make the method memory-friendly, we use the hybrid approximation of Leenaarts et al. (2012) for the redistribution integral. To make it fast, we use linear interpolation on equidistant frequency grids. We verify our algorithm against computations with the RH code and analyze it for stability, convergence, and usefulness of acceleration using model atoms of Mg II with the h&k lines and H I with the Ly {\alpha} line treated in PRD. A typical 3D PRD solution can be obtained in a model atmosphere with $252 \times 252 \times 496$ coordinate points in 50 000--200 000 CPU hours, which is a factor ten slower than computations assuming complete redistribution. We illustrate the importance of the joint action of PRD and 3D effects for the Mg II h&k lines for disk-center intensities as well as the center-to-limb variation. The proposed method allows simulating PRD lines in time series of radiation-MHD models in order to interpret observations of chromospheric lines at high spatial resolution.


title: Hysteresis of nanocylinders with Dzyaloshinskii-Moriya interaction

abstract: The potential for application of magnetic skyrmions in high density storage devices provides a strong drive to investigate and exploit their stability and manipulability. Through a three-dimensional micromagnetic hysteresis study, we investigate the question of existence of skyrmions in cylindrical nanostructures of variable thickness. We quantify the applied field and thickness dependence of skyrmion states, and show that these states can be accessed through relevant practical hysteresis loop measurement protocols. As skyrmionic states have yet to be observed experimentally in confined helimagnetic geometries, our work opens prospects for developing viable hysteresis process-based methodologies to access and observe skyrmionic states.


title: Adaptive Cross-Packet HARQ

abstract: In this work, we investigate a coding strategy devised to increase the throughput in hybrid ARQ (HARQ) transmission over block fading channel. In our approach, the transmitter jointly encodes a variable number of bits for each round of HARQ. The parameters (rates) of this joint coding can vary and may be based on the negative acknowledgment (NACK) provided by the receiver or, on the past (outdated) information about the channel states. These new degrees of freedom allow us to improve the match between the codebook and the channel states experienced by the receiver. The results indicate that significant gains can be obtained using the proposed coding strategy, particularly notable when the conventional HARQ fails to offer throughput improvement even if the number of transmission rounds is increased. The new cross-packet HARQ is also implemented using turbo codes where we show that the theoretically predicted throughput gains materialize in practice, and we discuss the implementation challenges.


title: Reconstructing a Bounded-Degree Directed Tree Using Path Queries

abstract: We present a randomized algorithm for reconstructing directed rooted trees of $n$ nodes and node degree at most $d$, by asking at most $O(dn\log^2 n)$ path queries. Each path query takes as input an origin node and a target node, and answers whether there is a directed path from the origin to the target. Regarding lower bounds, we show that any randomized algorithm requires at least $\Omega(n\log n)$ queries, while any deterministic algorithm requires at least $\Omega(dn)$ queries. Additionally, we present a $O(dn\log^3 n)$ randomized algorithm for noisy queries, and a $O(dn\log^2 n)$ randomized algorithm for additive queries on weighted trees.


title: A Complete Characterization of Determinantal Quadratic Polynomials

abstract: The problem of expressing a multivariate polynomial as the determinant of a monic (definite) symmetric or Hermitian linear matrix polynomial (LMP) has drawn a huge amount of attention due to its connection with optimization problems. In this paper we provide a necessary and sufficient condition for the existence of \textit{monic Hermitian determinantal representation} as well as \textit{monic symmetric determinantal representation} of size $2$ for a given quadratic polynomial. Further we propose a method to construct such a monic determinantal representtaion (MDR) of size $2$ if it exists. It is known that a quadratic polynomial $f(\x)=\x^{T}A\x+b^{T}\x+1$ has a symmetric MDR of size $n+1$ if $A$ is \textit{negative semidefinite}. We prove that if a quadratic polynomial $f(\x)$ with $A$ which is not negative semidefinite has an MDR of size greater than $2$, then it has an MDR of size $2$ too. We also characterize quadratic polynomials which exhibit diagonal MDRs.


title: Regularity of the level set flow

abstract: We showed earlier that the level set function of a monotonic advancing front is twice differentiable everywhere with bounded second derivative. We show here that the second derivative is continuous if and only if the flow has a single singular time where it becomes extinct and the singular set consists of a closed $C^1$ manifold with cylindrical singularities.


title: Characterization of Three High Efficiency and Blue Sensitive Silicon Photomultipliers

abstract: We report about the optical and electrical characterization of three high efficiency and blue sensitive Silicon photomultipliers from FBK, Hamamatsu, and SensL. Key features of the tested devices when operated at 90% breakdown probability are peak photon detection efficiencies between 40% and 55%, temperature dependencies of gain and PDE that are less than 1%/$^{\circ}$C, dark rates of $\sim$50\,kHz/mm$^{2}$ at room temperature, afterpulsing of about 2%, and direct optical crosstalk between 6% and 20%. The characteristics of all three devices impressively demonstrate how the Silicon-photomultiplier technology has improved over the past ten years. It is further demonstrated how the voltage and temperature characteristics of a number of quantities can be parameterized on the basis of physical models. The models provide a deeper understanding of the device characteristics over a wide bias and temperature range. They also serve as examples how producers could provide the characteristics of their SiPMs to users. A standardized parameterization of SiPMs would enable users to find the optimal SiPM for their application and the operating point of SiPMs without having to perform measurements thus significantly reducing design and development cycles.


title: Origin and evolution of two-component debris discs and an application to the q$^1$ Eridani system

abstract: Many debris discs reveal a two-component structure, with an outer Kuiper-belt analogue and a warm inner component whose origin is still a matter of debate. One possibility is that warm emission stems from an "asteroid belt" closer in to the star. We consider a scenario in which a set of giant planets is formed in an initially extended planetesimal disc. These planets carve a broad gap around their orbits, splitting up the disc into the outer and the inner belts. After the gas dispersal, both belts undergo collisional evolution in a steady-state regime. This scenario is explored with detailed collisional simulations involving realistic physics to describe a long-term collisional depletion of the two-component disc. We find that the inner disc may be able to retain larger amounts of material at older ages than thought before on the basis of simplified analytic models. We show that the proposed scenario is consistent with a suite of thermal emission and scattered light observational data for a bright two-temperature debris disc around a nearby solar-type star q$^1$ Eridani. This implies a Solar System-like architecture of the system, with an outer massive "Kuiper belt", an inner "asteroid belt", and a few Neptune- to Jupiter-mass planets in between.


title: A Framework for Optimal Matching for Causal Inference

abstract: We propose a novel framework for matching estimators for causal effect from observational data that is based on minimizing the dual norm of estimation error when expressed as an operator. We show that many popular matching estimators can be expressed as optimal in this framework, including nearest-neighbor matching, coarsened exact matching, and mean-matched sampling. This reveals their motivation and aptness as structural priors formulated by embedding the effect in a particular functional space. This also gives rise to a range of new, kernel-based matching estimators that arise when one embeds the effect in a reproducing kernel Hilbert space. Depending on the case, these estimators can be found using either quadratic optimization or integer optimization. We show that estimators based on universal kernels are universally consistent without model specification. In empirical results using both synthetic and real data, the new, kernel-based estimators outperform all standard causal estimators in estimation error.


title: On the Astrid asteroid family

abstract: Among asteroid families, the Astrid family is peculiar because of its unusual inclination distribution. Objects at $a\simeq$~2.764 au are quite dispersed in this orbital element, giving the family a "crab-like" appearance. Recent works showed that this feature is caused by the interaction of the family with the $s-s_C$ nodal secular resonance with Ceres, that spreads the inclination of asteroids near its separatrix. As a consequence, the currently observed distribution of the $v_W$ component of terminal ejection velocities obtained from inverting Gauss equation is quite leptokurtic, since this parameter mostly depends on the asteroids inclination. The peculiar orbital configuration of the Astrid family can be used to set constraints on key parameters describing the strength of the Yarkovsky force, such as the bulk and surface density and the thermal conductivity of surface material. By simulating various fictitious families with different values of these parameters, and by demanding that the current value of the kurtosis of the distribution in $v_W$ be reached over the estimated lifetime of the family, we obtained that the thermal conductivity of Astrid family members should be $\simeq$ 0.001 W/m/K, and that the surface and bulk density should be higher than 1000 kg/m$^{3}$. Monte Carlo methods simulating Yarkovsky and stochastic YORP evolution of the Astrid family show its age to be $T$ = 140$\pm$30 Myr old, in good agreement with estimates from other groups. Its terminal ejection velocity parameter is in the range $V_{EJ}= 5^{+17}_{-5}$~m/s.] Values of $V_{EJ}$ larger than 25 m/s are excluded from constraints from the current inclination distribution.


title: Formation of the G-ring arc

abstract: Since 2004, the images obtained by Cassini spacecraft's on-board cameras have revealed the existence of several small satellites in the Saturn system. Some of these small satellites are embedded in arcs of particles. While these satellites and their arcs are known to be in corotation resonances with Mimas, their origin remains unknown. This work investigates one possible process for capturing bodies into a corotation resonance, which involves raising the eccentricity of a perturbing body. Therefore, through numerical simulations and analytical studies, we show a scenario that the excitation of Mimas' eccentricity could capture particles in a corotation resonance and given a possible explanation for the origin for the arcs.


title: The Picard group of motivic A(1)

abstract: We show that the Picard group $Pic(A(1))$ of the stable category of modules over $\mathbb{C}$-motivic $A(1)$ is isomorphic to $\mathbb{Z}^4$. By comparison, the Picard group of classical $A(1)$ is $\mathbb{Z}^2 \oplus \mathbb{Z}/2$. One extra copy of $\mathbb{Z}$ arises from the motivic bigrading. The joker is a well-known exotic element of order $2$ in the Picard group of classical $A(1)$. The $\mathbb{C}$-motivic joker has infinite order.


title: Polarimetry of comets 67P/Churyumov-Gerasimenko, 74P/Smirnova-Chernykh, and 152P/Helin-Lawrence

abstract: Aims. Polarimetric characteristics of comets at large heliocentric distances is a relatively unexplored area; we extend the idea by adding and analysing the data for three Jupiter family comets (JFCs). Methods. With the FORS2 instrument of the ESO VLT, we performed quasi-simultaneous photometry and polarimetry of three active JFCs 67P/Churyumov-Gerasimenko, 74P/Smirnova-Chernykh, and 152P/Helin-Lawrence. Results. We obtained in total 23 polarimetric measurements at different epochs, covering a phase-angle range ~1 -16 degrees and heliocentric distances from 3 to 4.6 au. From our observations we obtained both colour and polarimetric maps to look for structures in the comae and tails of the comets. Conclusions. 74P/Smirnova-Chernykh and 152P/Helin-Lawrence exhibit enough activity at large heliocentric distances to be detectable in polarimetric measurements. Polarimetric and colour maps indicate no evidence of dust particle


title: Detailed study of the microwave emission of the supernova remnant 3C 396

abstract: We have observed the supernova remnant 3C~396 in the microwave region using the Parkes 64-m telescope. Observations have been made at 8.4 GHz, 13.5 GHz, and 18.6 GHz and in polarisation at 21.5 GHz. We have used data from several other observatories, including previously unpublished observations performed by the Green Bank Telescope at 31.2 GHz, to investigate the nature of the microwave emission of 3C 396. Results show a spectral energy distribution dominated by a single component power law emission with $\alpha=(-0.364 \pm 0.017)$. Data do not favour the presence of anomalous microwave emission coming from the source. Polarised emission at 21.5 GHz is consistent with synchrotron-dominated emission. We present microwave maps and correlate them with infrared (IR) maps in order to characterise the interplay between thermal dust and microwave emission. IR vs. microwave TT plots reveal poor correlation between mid-infrared and microwave emission from the core of the source. On the other hand, a correlation is detected in the tail emission of the outer shell of 3C 396, which could be ascribed to Galactic contamination.


title: Asymptotic Theory of Channeling in the Field of an Atomic Chain and an Atomic Plane

abstract: A rigorous theory of diffraction scattering from extended objects is proposed. The present theory is based on a multiple asymptotic expansion of an integral equation for the exact wave function in terms of the large parameters of the problem, which are the range of the potential and the momentum components of the incident particle. For small angles of incidence the density of positively charged particles on the axis of a chain is always lower than unit y and the density of negatively charged particles has a maximum for certain strength of the potential which can exceed considerably unity. The conditions of validity of the proposed approach are obtained.


title: Effective Chabauty for symmetric powers of curves

abstract: Faltings' theorem states that curves of genus $g \geq 2$ have finitely many rational points. Using the ideas of Faltings, Mumford, Parshin and Raynaud, one obtains an upper bound on the number of rational points, but this bound is too large to be used in any reasonable sense. In 1985, Coleman showed that Chabauty's method, which works when the Mordell-Weil rank of the Jacobian of the curve is smaller than $g$, can be used to give a good effective bound on the number of rational points of curves of genus $g > 1$. We draw ideas from nonarchimedean geometry to show that we can also give an effective bound on the number of rational points outside of the special set of $d$-th symmetric power of $X$, where $X$ is a curve of genus $g > d$, when the Mordell-Weil rank of the Jacobian of the curve is at most $g-d$ and the curve further satisfies certain rigid analytic conditions.


title: A Goldilocks principle for modeling radial velocity noise

abstract: The doppler measurements of stars are diluted and distorted by stellar activity noise. Different choices of noise models and statistical methods have led to much controversy in the confirmation of exoplanet candidates obtained through analysing radial velocity data. To quantify the limitation of various models and methods, we compare different noise models and signal detection criteria for various simulated and real data sets in the Bayesian framework. According to our analyses, the white noise model tend to interpret noise as signal, leading to false positives. On the other hand, the red noise models are likely to interprete signal as noise, resulting in false negatives. We find that the Bayesian information criterion combined with a Bayes factor threshold of 150 can efficiently rule out false positives and confirm true detections. We further propose a Goldilocks principle aimed at modeling radial velocity noise to avoid too many false positives and too many false negatives. We propose that the noise model with RHK-dependent jitter is used in combination with the moving average model to detect planetary signals for M dwarfs. Our work may also shed light on the noise modeling for hotter stars, and provide a valid approach for finding similar principles in other disciplines.


title: On Properties of Differential Inclusions with Prox-regular Sets

abstract: In this paper, some regularity properties of solutions of the following differential inclusion \begin{equation}\nonumber \left\{ \begin{array}{l} \dot{x}(t) \in f\big(x(t)\big) -N_{C}\big(x(t)\big)\; {\rm a.e.} \; t \in [0,+\infty), x(0) = x_0\in C, \end{array}\right. \end{equation} are analyzed where $f: H\to H$ is Lipschitz continuous and $C$ is closed, uniformly prox-regular subset of a Hilbet space $H$. Here $N_{C}(\cdot)$ denotes the proximal normal cone of $C$. This work can be considered as an improvement of [Hantoute-Mazade 2013] since these properties are established without the additional tangential condition at each point in $C$.


title: LP-Based Robust Algorithms for Noisy Minor-Free and Bounded Treewidth Graphs

abstract: We give a general approach for solving optimization problems on noisy minor free graphs, where a \delta-fraction of edges and vertices are adversarially corrupted. The noisy setting was first considered by Magen and Moharrami and they gave a (1 + \delta)-estimation algorithm for the independent set problem. Later, Chan and Har-Peled designed a local search algorithm that finds a (1 + O(\delta))-approximate independent set. However, nothing was known regarding other problems in the noisy setting. Our main contribution is a general LP-based framework that yields a (1 + O(\delta log m log log m))-approximation algorithm for noisy MAX-k-CSPs on m clauses.


title: Tuning of Near- and Far-Field Properties of All-dielectric Dimer Nanoantennas via Ultrafast Electron-Hole Plasma Photoexcitation

abstract: Achievement of all-optical ultrafast signal modulation and routing by a low-loss nanodevice is a crucial step towards an ultracompact optical chip with high performance. Here, we propose a specifically designed silicon dimer nanoantenna, which is tunable via photoexcitation of dense electron-hole plasma with ultrafast relaxation rate. Basing on this concept, we demonstrate the effect of beam steering up to 20 degrees via simple variation of incident intensity, being suitable for ultrafast light routing in an optical chip. The effect is demonstrated both in the visible and near-IR spectral regions for silicon and germanium based nanoantennas. We also reveal the effect of electron-hole plasma photoexcitation on local density of states (LDOS) in the dimer gap and find that the orientation averaged LDOS can be altered by 50\%, whereas modification of the projected LDOS can be even more dramatic: almost 500\% for transverse dipole orientation. Moreover, our analytical model sheds light on transient dynamics of the studied nonlinear nanoantennas, yielding all temporal characteristics of the proposed ultrafast nanodevice. The proposed concept paves the ways to creation of low-loss, ultrafast, and compact devices for optical signal modulation and routing.


title: High-speed real-time single-pixel microscopy based on Fourier sampling

abstract: Single-pixel cameras based on the concepts of compressed sensing (CS) leverage the inherent structure of images to retrieve them with far fewer measurements and operate efficiently over a significantly broader spectral range than conventional silicon-based cameras. Recently, photonic time-stretch (PTS) technique facilitates the emergence of high-speed single-pixel cameras. A significant breakthrough in imaging speed of single-pixel cameras enables observation of fast dynamic phenomena. However, according to CS theory, image reconstruction is an iterative process that consumes enormous amounts of computational time and cannot be performed in real time. To address this challenge, we propose a novel single-pixel imaging technique that can produce high-quality images through rapid acquisition of their effective spatial Fourier spectrum. We employ phase-shifting sinusoidal structured illumination instead of random illumination for spectrum acquisition and apply inverse Fourier transform to the obtained spectrum for image restoration. We evaluate the performance of our prototype system by recognizing quick response (QR) codes and flow cytometric screening of cells. A frame rate of 625 kHz and a compression ratio of 10% are experimentally demonstrated in accordance with the recognition rate of the QR code. An imaging flow cytometer enabling high-content screening with an unprecedented throughput of 100,000 cells/s is also demonstrated. For real-time imaging applications, the proposed single-pixel microscope can significantly reduce the time required for image reconstruction by two orders of magnitude, which can be widely applied in industrial quality control and label-free biomedical imaging.


title: Assessing and tuning brain decoders: cross-validation, caveats, and guidelines

abstract: Decoding, ie prediction from brain images or signals, calls for empirical evaluation of its predictive power. Such evaluation is achieved via cross-validation, a method also used to tune decoders' hyper-parameters. This paper is a review on cross-validation procedures for decoding in neuroimaging. It includes a didactic overview of the relevant theoretical considerations. Practical aspects are highlighted with an extensive empirical study of the common decoders in within-and across-subject predictions, on multiple datasets --anatomical and functional MRI and MEG-- and simulations. Theory and experiments outline that the popular " leave-one-out " strategy leads to unstable and biased estimates, and a repeated random splits method should be preferred. Experiments outline the large error bars of cross-validation in neuroimaging settings: typical confidence intervals of 10%. Nested cross-validation can tune decoders' parameters while avoiding circularity bias. However we find that it can be more favorable to use sane defaults, in particular for non-sparse decoders.


title: Verification of the helioseismic Fourier-Legendre analysis for meridional flow measurements

abstract: Measuring the Sun's internal meridional flow is one of the key issues of helioseismology. Using the Fourier-Legendre analysis is a technique for addressing this problem. We validate this technique with the help of artificial helioseismic data. The analysed data set was obtained by numerically simulating the effect of the meridional flow on the seismic wave field in the full volume of the Sun. In this way, a 51.2-hour long time series was generated. The resulting surface velocity field is then analyzed in various settings: Two $360^\circ \times 90^\circ$ halfspheres, two $120^\circ \times 60^\circ$ patches on the front and farside of the Sun (North and South, respectively) and two $120^\circ \times 60^\circ$ patches on the northern and southern frontside only. We compare two possible measurement setups: observations from Earth and from an additional spacecraft on the solar farside, and observations from Earth only, in which case the full information of the global solar oscillation wave field was available. We find that, with decreasing observing area, the accessible depth range decreases: the $360^\circ \times 90^\circ$ view allows us to probe the meridional flow almost to the bottom of the convection zone, while the $120^\circ \times 60^\circ$ view means only the outer layers can be probed. These results confirm the validity of the Fourier-Legendre analysis technique for helioseismology of the meridional flow. Furthermore these flows are of special interest for missions like Solar Orbiter that promises to complement standard helioseismic measurements from the solar nearside with farside observations.


title: A new asymmetric generalisation of the t-distribution

abstract: A 6-parameter fat-tailed distribution is proposed that generalises the t-distribution and allows asymmetry of scale and also of tail power, whilst avoiding the discontinuity of the second derivative of the split-t (AST) distribution. With the sixth parameter set to unity and no asymmetry, the distribution reduces to a t-distribution, but with the sixth parameter reduced, fatter tails than those of the t-distribution are allowed (the tails start earlier) and the distribution generalises Johnson's $S_U$ distribution. Data fitting is illustrated with examples.


title: Multiple imputation of missing covariate values in multilevel models with random slopes: A cautionary note

abstract: Multiple imputation (MI) has become one of the main procedures used to treat missing data, but the guidelines from the methodological literature are not easily transferred to multilevel research. For models including random slopes, proper MI can be difficult, especially when the covariate values are partially missing. In the present article, we discuss applications of MI in multilevel random-coefficient models, theoretical challenges posed by slope variation, and the current limitations of standard MI software. Our findings from three simulation studies suggest that (a) MI is able to recover most parameters, but is currently not well suited to capture slope variation entirely when covariate values are missing; (b) MI offers reasonable estimates for most parameters, even in smaller samples or when its assumptions are not met; and (c) listwise deletion can be an alternative worth considering when preserving the slope variance is particularly important.


title: Spectral theory for structured perturbations of linear operators

abstract: We characterize the spectrum (and its parts) of operators which can be represented as G=A+BC for a simpler operator A and a structured perturbation BC. The interest in this kind of perturbations is motivated, e.g., by perturbations of the domain of an operator A but also arises in the theory of closed-loop systems in control theory. In many cases our results yield the spectral values of G as zeros of a "characteristic equation".


title: On a geometric inequality related to fractional integration

abstract: In this paper we consider a new kind of inequality related to fractional integration, motivated by Gressman's paper. Based on it we investigate its multilinear analogue inequalities. Combining with the Gressman's work on multilinear integral, we establish this new kind of geometric inequalities with bilinear form and multilinear form in more general settings. Moreover, in some cases we also find the best constants and optimisers for these geometric inequalities on Euclidean spaces with Lebesgue measure settings with $L^p$ bounds.


title: Stationary Schr\"odinger equation in the semi-classical limit: numerical coupling of oscillatory and evanescent regions

abstract: This paper is concerned with a 1D Schr\"odinger scattering problem involving both oscillatory and evanescent regimes, separated by jump discontinuities in the potential function, to avoid "turning points". We derive a non-overlapping domain decomposition method to split the original problem into sub-problems on these regions, both for the continuous and afterwards for the discrete problem. Further, a hybrid WKB-based numerical method is designed for its efficient and accurate solution in the semi-classical limit: a WKB-marching method for the oscillatory regions and a FEM with WKB-basis functions in the evanescent regions. We provide a complete error analysis of this hybrid method and illustrate our convergence results by numerical tests.


title: On some determinant and matrix inequalities with a geometrical flavour

abstract: In this paper we study some determinant inequalities and matrix inequalities which have a geometrical flavour. We first examine some inequalities which place work of Macbeath [13] in a more general setting and also relate to recent work of Gressman [8]. In particular, we establish optimisers for these determinant inequalities. We then use these inequalities to establish our main theorem which gives a geometric inequality of matrix type which improves and extends some inequalities of Christ in [5].


title: Cost Effective Campaigning in Social Networks

abstract: Campaigners are increasingly using online social networking platforms for promoting products, ideas and information. A popular method of promoting a product or even an idea is incentivizing individuals to evangelize the idea vigorously by providing them with referral rewards in the form of discounts, cash backs, or social recognition. Due to budget constraints on scarce resources such as money and manpower, it may not be possible to provide incentives for the entire population, and hence incentives need to be allocated judiciously to appropriate individuals for ensuring the highest possible outreach size. We aim to do the same by formulating and solving an optimization problem using percolation theory. In particular, we compute the set of individuals that are provided incentives for minimizing the expected cost while ensuring a given outreach size. We also solve the problem of computing the set of individuals to be incentivized for maximizing the outreach size for given cost budget. The optimization problem turns out to be non trivial; it involves quantities that need to be computed by numerically solving a fixed point equation. Our primary contribution is, that for a fairly general cost structure, we show that the optimization problems can be solved by solving a simple linear program. We believe that our approach of using percolation theory to formulate an optimization problem is the first of its kind.


title: Weighted Online Problems with Advice

abstract: Recently, the first online complexity class, AOC, was introduced. The class consists of many online problems where each request must be either accepted or rejected, and the aim is to either minimize or maximize the number of accepted requests, while maintaining a feasible solution. All AOC-complete problems (including Independent Set, Vertex Cover, Dominating Set, and Set Cover) have essentially the same advice complexity. In this paper, we study weighted versions of problems in AOC, i.e., each request comes with a weight and the aim is to either minimize or maximize the total weight of the accepted requests. In contrast to the unweighted versions, we show that there is a significant difference in the advice complexity of complete minimization and maximization problems. We also show that our algorithmic techniques for dealing with weighted requests can be extended to work for non-complete AOC problems such as maximum matching (giving better results than what follow from the general AOC results) and even non-AOC problems such as scheduling.


title: Singular inflation from Born-Infeld-f(R) gravity

abstract: Accelerating dynamics from Born-Infeld-$f(R)$ gravity are studied in a simplified conformal approach without matter. Explicit unification of inflation with late-time acceleration is realized within this singular inflation approach, which is similar to Odintsov-Oikonomou singular $f(R)$ inflation. Our model turns out to be consistent with the latest release of Planck data.


title: Inertial displacement of a domain wall excited by ultra-short circularly polarized laser pulses

abstract: Domain wall motion driven by ultra-short laser pulses is a prerequisite for envisaged low-power spintronics combining storage of information in magneto electronic devices with high speed and long distance transmission of information encoded in circularly polarized light. Here we demonstrate the conversion of the circular polarization of incident femtosecond laser pulses into inertial displacement of a domain wall in a ferromagnetic semiconductor. In our study we combine electrical measurements and magneto-optical imaging of the domain wall displacement with micromagnetic simulations. The optical spin transfer torque acts over a picosecond recombination time of the spin polarized photo-carriers which only leads to a deformation of the internal domain wall structure. We show that subsequent depinning and micro-meter distance displacement without an applied magnetic field or any other external stimuli can only occur due to the inertia of the domain wall.


title: Blue straggler formation at core collapse

abstract: Among the most striking feature of blue straggler stars (BSS) is the presence of multiple sequences of BSSs in the colour-magnitude diagrams (CMDs) of several globular clusters. It is often envisaged that such a multiple BSS sequence would arise due a recent core collapse of the host cluster, triggering a number of stellar collisions and binary mass transfers simultaneously over a brief episode of time. Here we examine this scenario using direct N-body computations of moderately massive star clusters (of order 10^4 Msun ). As a preliminary attempt, these models are initiated with approx. 8-10 Gyr old stellar population and King profiles of high concentrations, being "tuned" to undergo core collapse quickly. BSSs are indeed found to form in a "burst" at the onset of the core collapse and several of such BS-bursts occur during the post-core-collapse phase. In those models that include a few percent primordial binaries, both collisional and binary BSSs form after the onset of the (near) core-collapse. However, there is as such no clear discrimination between the two types of BSSs in the corresponding computed CMDs. We note that this may be due to the less number of BSSs formed in these less massive models than that in actual globular clusters.


title: The maximum of the minimal multiplicity of eigenvalues of symmetric matrices whose pattern is constrained by a graph

abstract: In this paper we introduce a parameter $Mm(G)$, defined as the maximum over the minimal multiplicities of eigenvalues among all symmetric matrices corresponding to a graph $G$. We compute $Mm(G)$ for several families of graphs.


title: Stellar populations in the outskirts of M31: the mid-infrared view

abstract: The mid-infrared provides a unique view of galaxy stellar populations, sensitive to both the integrated light of old, low-mass stars and to individual dusty mass-losing stars. We present results from an extended Spitzer/IRAC survey of M31 with total lengths of 6.6 and 4.4 degrees along the major and minor axes, respectively. The integrated surface brightness profile proves to be surprisingly diffcult to trace in the outskirts of the galaxy, but we can also investigate the disk/halo transition via a star count profile, with careful correction for foreground and background contamination. Our point-source catalog allows us to report on mid-infrared properties of individual objects in the outskirts of M31, via cross-correlation with PAndAS, WISE, and other catalogs.


title: On the initial boundary value problem of a Navier-Stokes/$Q$-tensor model for liquid crystals

abstract: This work is concerned with the solvability of a Navier-Stokes/$Q$-tensor coupled system modeling the nematic liquid crystal flow on a bounded domain in three dimensional Euclidian space with strong anchoring boundary condition for the order parameter. We prove the existence of local in time strong solution to the system with the anisotropic elastic energy. The proof is based on mainly two ingredients: first, we show that the Euler-Lagrange operator corresponding to the Landau-de Gennes free energy with general elastic coefficients fulfills the strong Legendre condition. This result together with a higher order energy estimate leads to the well-posedness of the linearized system, and then a local in time solution of the original system which is regular in temporal variable follows via a fixed point argument. Secondly, the hydrodynamic part of the coupled system can be reformulated into a quasi-stationary Stokes type equation to which the regularity theory of the generalized Stokes system, and then a bootstrap argument can be applied to enhance the spatial regularity of the local in time solution.


title: LOLAS-2 : redesign of an optical turbulence profiler

abstract: We present the development, tests and first results of the second generation Low Layer Scidar (LOLAS-2). This instrument constitutes a strongly improved version of the prototype Low Layer Scidar, which is aimed at the measurement of optical turbulence profiles close to the ground, with high altitude-resolution. The method is based on the Generalised Scidar principle which consists in taking double-star scintillation images on a defocused pupil plane and calculating in real time the autocovariance of the scintillation. The main components are an open-truss 40-cm Ritchey-Chr\'etien telescope, a german-type equatorial mount, an Electron Multiplying CCD camera and a dedicated acquisition and real-time data processing software. The new optical design of LOLAS-2 is significantly simplified compared with the prototype. The experiments carried out to test the permanence of the image within the useful zone of the detector and the stability of the telescope focus show that LOLAS-2 can function without the use of the autoguiding and autofocus algorithms that were developed for the prototype version. Optical turbulence profiles obtained with the new Low Layer Scidar have the best altitude-resolution ever achieved with Scidar-like techniques (6.3 m). The simplification of the optical layout and the improved mechanical properties of the telescope and mount make of LOLAS-2 a more robust instrument.


title: Linear radial Regge trajectories for mesons with any quark flavor

abstract: In the Regge phenomenology, the radial spectrum of light mesons is given by a linear relation $M_n^2=a(n+b)$, where $a$ is a universal slope, the dimensionless intercept $b$ depends on quantum numbers, and $n$ enumerates the excited states in radial recurrences. The usual extensions of this relation to heavy quarkonia in the framework of hadron string models typically lead to strong nonlinearities which seem to be at variance with the available experimental data. Introducing a radially static string picture of mesons, we put forward a linear generalization $(M_n-m_1-m_2)^2=a(n+b)$, where $m_{1,2}$ are quark masses. The vector channel contains enough experimental states to check this new relation and a good agreement is observed. It is shown that this generalization leads to a simple estimate of current quark masses from the radial spectra.


title: Generalized Weyl modules for twisted current algebras

abstract: We introduce the notion of generalized Weyl modules for twisted current algebras. We study their representation-theoretic and combinatorial properties and connection to the theory of nonsymmetric Macdonald polynomials. As an application we compute the dimension of the classical Weyl modules in the remaining unknown case.


title: Magnetic Properties from the Viewpoints of Electronic Hamiltonian: Spin Exchange Parameters, Spin Orientation and Spin-Half Misconception

abstract: In this chapter we review the quantitative and qualitative aspects of describing the properties of magnetic solids on the basis of electronic Hamiltonian. We show that a spin Hamiltonian approach becomes consistent with an electronic Hamiltonian approach if the spin lattice and its associated spin exchange parameters, to be used for the spin Hamiltonian, are determined by the energy-mapping analysis based on DFT calculations. The preferred spin orientation (i.e., the magnetic anisotropy) of a magnetic ion is not predicted by a spin Hamiltonian because it does not include the orbital degree of freedom explicitly. In contrast, the magnetic anisotropy is readily predicted by electronic structure theories employing both orbital and spin degrees of freedom, if one takes into consideration the spin-orbit coupling (SOC). It was shown that the preferred spin orientation of a magnetic ion can be predicted and understood in terms of the HOMO-LUMO interactions of the magnetic ion by taking SOC as perturbation. A spin Hamiltonian gives rise to the spin-half misconception, namely, the blind belief that spin-half magnetic ions do not possess magnetic anisotropy that arise from SOC. This misconception is a direct consequence from the limitedness of a spin Hamiltonian that it lacks the orbital degree of freedom. We show that the magnetic properties of 5d ion oxides are better explained by the LS-coupling than by the jj-coupling scheme of SOC, that the spin-orbital entanglement of 5d ions is not as strong as has been assumed.


title: On the Effective Field Theory of Heterotic Vacua

abstract: The effective field theory of heterotic vacua that realise $\mathbb{R}^{3,1}$ preserving $\mathcal{N} =1$ supersymmetry are studied. The vacua in question admit large radius limits taking the form $\mathbb{R}^{3,1}\times {X}$ , with ${X}$ a smooth three-fold with vanishing first Chern class and a stable holomorphic gauge bundle $\mathcal{E}$. In a previous paper we calculated the kinetic terms for moduli, deducing the moduli metric and Kahler potential. In this paper, we compute the remaining couplings in the effective field theory, correct to first order in alpha prime. In particular, we compute the contribution of the matter sector to the Kahler potential, derive the Yukawa couplings and other quadratic fermionic couplings. From this we write down a Kahler potential $\mathcal{K}$ and superpotential $\mathcal{W}$ .


title: On point interactions realised as Ter-Martirosyan-Skornyakov Hamiltonians

abstract: For quantum systems of zero-range interaction we discuss the mathematical scheme within which modelling the two-body interaction by means of the physically relevant ultra-violet asymptotics known as the "Ter-Martirosyan-Skornyakov condition" gives rise to a self-adjoint realisation of the corresponding Hamiltonian. This is done within the self-adjoint extension scheme of Krein, Visik, and Birman. We show that the Ter-Martirosyan-Skornyakov asymptotics is a condition of self-adjointness only when is imposed in suitable functional spaces, and not just as a pointwise asymptotics, and we discuss the consequences of this fact on a model of two identical fermions and a third particle of different nature.


title: Guarded Cubical Type Theory: Path Equality for Guarded Recursion

abstract: This paper improves the treatment of equality in guarded dependent type theory (GDTT), by combining it with cubical type theory (CTT). GDTT is an extensional type theory with guarded recursive types, which are useful for building models of program logics, and for programming and reasoning with coinductive types. We wish to implement GDTT with decidable type-checking, while still supporting non-trivial equality proofs that reason about the extensions of guarded recursive constructions. CTT is a variation of Martin-L\"of type theory in which the identity type is replaced by abstract paths between terms. CTT provides a computational interpretation of functional extensionality, is conjectured to have decidable type checking, and has an implemented type-checker. Our new type theory, called guarded cubical type theory, provides a computational interpretation of extensionality for guarded recursive types. This further expands the foundations of CTT as a basis for formalisation in mathematics and computer science. We present examples to demonstrate the expressivity of our type theory, all of which have been checked using a prototype type-checker implementation, and present semantics in a presheaf category.


title: Waterfront on the Martian Planitia: Algorithmic emergent catchments on disordered terrain

abstract: Under a terraforming scenario, a reactivated hydrological cycle on Mars will result in upwards movement of water due to evaporation and precipitation. If Mars' embryonic fossilized catchments provide inadequate drainage, Mars' limited supplies of water may be absorbed entirely by crater lakes and glaciers, with negative consequences for the terraforming effort. We demonstrate a stable, convergent algorithm for the efficient modeling of water flow over disordered terrain. This model is applied to Mars Orbital Laser Altimeter data and successfully predicts the formation of fossilized waterways and canyons visible only at much higher resolution. This exploratory study suggests that despite its impossibly rugged appearance, ancient water flows have carved channels that provide effective drainage over the majority of Mars' surface. We also provide one possible reconstruction of a terraformed surface water distribution.


title: Geometric Median in Nearly Linear Time

abstract: In this paper we provide faster algorithms for solving the geometric median problem: given $n$ points in $\mathbb{R}^{d}$ compute a point that minimizes the sum of Euclidean distances to the points. This is one of the oldest non-trivial problems in computational geometry yet despite an abundance of research the previous fastest algorithms for computing a $(1+\epsilon)$-approximate geometric median were $O(d\cdot n^{4/3}\epsilon^{-8/3})$ by Chin et. al, $\tilde{O}(d\exp{\epsilon^{-4}\log\epsilon^{-1}})$ by Badoiu et. al, $O(nd+\mathrm{poly}(d,\epsilon^{-1})$ by Feldman and Langberg, and $O((nd)^{O(1)}\log\frac{1}{\epsilon})$ by Parrilo and Sturmfels and Xue and Ye. In this paper we show how to compute a $(1+\epsilon)$-approximate geometric median in time $O(nd\log^{3}\frac{1}{\epsilon})$ and $O(d\epsilon^{-2})$. While our $O(d\epsilon^{-2})$ is a fairly straightforward application of stochastic subgradient descent, our $O(nd\log^{3}\frac{1}{\epsilon})$ time algorithm is a novel long step interior point method. To achieve this running time we start with a simple $O((nd)^{O(1)}\log\frac{1}{\epsilon})$ time interior point method and show how to improve it, ultimately building an algorithm that is quite non-standard from the perspective of interior point literature. Our result is one of very few cases we are aware of outperforming traditional interior point theory and the only we are aware of using interior point methods to obtain a nearly linear time algorithm for a canonical optimization problem that traditionally requires superlinear time. We hope our work leads to further improvements in this line of research.


title: Calculation of thermal conductivity coefficients of electrons in magnetized dense matter

abstract: The solution of Boltzmann equation for plasma in magnetic field, with arbitrarily degenerate electrons and non-degenerate nuclei, is obtained by Chapman-Enskog method. Functions, generalizing Sonin polynomials are used for obtaining an approximate solution. Fully ionized plasma is considered. The tensor of the heat conductivity coefficients in non-quantized magnetic field is calculated. For non-degenerate and strongly degenerate plasma the asymptotic analytic formulas are obtained, which are compared with results of previous authors. The Lorentz approximation, with neglecting of electron-electron encounters, is asymptotically exact for strongly degenerate plasma. We obtain, for the first time, in three polynomial approximation, with account of electron-electron collisions, analytical expressions for the heat conductivity tensor for non-degenerate electrons, in presence of a magnetic field. Account of the third polynomial improved substantially the precision of results. In two polynomial approximation our solution coincides with the published results. For strongly degenerate electrons we obtain, for the first time, an asymptotically exact analytical solution for the heat conductivity tensor in presence of a magnetic field. This solution has considerably more complicated dependence on the magnetic field than those in previous publications, and gives several times smaller relative value of a thermal conductivity across the magnetic field at $\omega\tau\gtrsim 0.8$.


title: Critical slowing down and hyperuniformity on approach to jamming

abstract: Hyperuniformity characterizes a state of matter that is poised at a critical point at which density or volume-fraction fluctuations are anomalously suppressed at infinite wavelengths. Recently, much attention has been given to the link between strict jamming and hyperuniformity in frictionless hard-particle packings. Doing so requires one to study very large packings, which can be difficult to jam properly. We modify the rigorous linear programming method of Donev et al. [J. Comp. Phys. 197, 139 (2004)] in order to test for jamming in putatively jammed packings of hard-disks in two dimensions. We find that various standard packing protocols struggle to reliably create packings that are jammed for even modest system sizes; importantly, these packings appear to be jammed by conventional tests. We present evidence that suggests that deviations from hyperuniformity in putative maximally random jammed (MRJ) packings can in part be explained by a shortcoming in generating exactly-jammed configurations due to a type of "critical slowing down" as the necessary rearrangements become difficult to realize by numerical protocols. Additionally, various protocols are able to produce packings exhibiting hyperuniformity to different extents, but this is because certain protocols are better able to approach exactly-jammed configurations. Nonetheless, while one should not generally expect exact hyperuniformity for disordered packings with rattlers, we find that when jamming is ensured, our packings are very nearly hyperuniform, and deviations from hyperuniformity correlate with an inability to ensure jamming, suggesting that strict jamming and hyperuniformity are indeed linked. This raises the possibility that the ideal MRJ packings have no rattlers. Our work provides the impetus for the development of packing algorithms that produce large disordered strictly jammed packings that are rattler-free.


title: How many faces can be recognized? Performance extrapolation for multi-class classification

abstract: The difficulty of multi-class classification generally increases with the number of classes. Using data from a subset of the classes, can we predict how well a classifier will scale with an increased number of classes? Under the assumption that the classes are sampled exchangeably, and under the assumption that the classifier is generative (e.g. QDA or Naive Bayes), we show that the expected accuracy when the classifier is trained on $k$ classes is the $k-1$st moment of a \emph{conditional accuracy distribution}, which can be estimated from data. This provides the theoretical foundation for performance extrapolation based on pseudolikelihood, unbiased estimation, and high-dimensional asymptotics. We investigate the robustness of our methods to non-generative classifiers in simulations and one optical character recognition example.


title: Estimating mutual information in high dimensions via classification error

abstract: Multivariate pattern analyses approaches in neuroimaging are fundamentally concerned with investigating the quantity and type of information processed by various regions of the human brain; typically, estimates of classification accuracy are used to quantify information. While a extensive and powerful library of methods can be applied to train and assess classifiers, it is not always clear how to use the resulting measures of classification performance to draw scientific conclusions: e.g. for the purpose of evaluating redundancy between brain regions. An additional confound for interpreting classification performance is the dependence of the error rate on the number and choice of distinct classes obtained for the classification task. In contrast, mutual information is a quantity defined independently of the experimental design, and has ideal properties for comparative analyses. Unfortunately, estimating the mutual information based on observations becomes statistically infeasible in high dimensions without some kind of assumption or prior. In this paper, we construct a novel classification-based estimator of mutual information based on high-dimensional asymptotics. We show that in a particular limiting regime, the mutual information is an invertible function of the expected $k$-class Bayes error. While the theory is based on a large-sample, high-dimensional limit, we demonstrate through simulations that our proposed estimator has superior performance to the alternatives in problems of moderate dimensionality.


title: Superconformal SU(1,1|n) mechanics

abstract: Recent years have seen an upsurge of interest in dynamical realizations of the superconformal group SU(1,1|2) in mechanics. Remarking that SU(1,1|2) is a particular member of a chain of supergroups SU(1,1|n) parametrized by an integer n, here we begin a systematic study of SU(1,1|n) multi-particle mechanics. A representation of the superconformal algebra su(1,1|n) is constructed on the phase space spanned by m copies of the (1,2n,2n-1) supermultiplet. We show that the dynamics is governed by two prepotentials V and F, and the Witten-Dijkgraaf-Verlinde-Verlinde equation for F shows up as a consequence of a more general fourth-order equation. All solutions to the latter in terms of root systems reveal decoupled models only. An extension of the dynamical content of the (1,2n,2n-1) supermultiplet by angular variables in a way similar to the SU(1,1|2) case is problematic.


title: Theory of the Lamb shift and Fine Structure in muonic $\mathrm{^4He}$ ions and the muonic $\mathrm{^3He-^4He}$ Isotope Shift

abstract: We provide an up to date summary of the theory contributions to the 2S-2P Lamb shift and the fine structure of the 2P state in the muonic helium ion $(\mathrm{\mu^4He})^+$. This summary serves as the basis for the extraction of the alpha particle charge radius from the muonic helium Lamb shift measurements at the Paul Scherrer Institute, Switzerland. Individual theory contributions needed for a charge radius extraction are compared and compiled into a consistent summary. The influence of the alpha particle charge distribution on the elastic two-photon exchange is studied to take into account possible model-dependencies of the energy levels on the electric form factor of the nucleus. We also discuss the theory uncertainty which enters the extraction of the $\mathrm{^3He-^4He}$ isotope shift from the muonic measurements. The theory uncertainty of the extraction is much smaller than a present discrepancy between previous isotope shift measurements. This work completes our series of $n=2$ theory compilations in light muonic atoms which we have performed already for muonic hydrogen, deuterium, and helium-3 ions.


title: Matters of Gravity, The Newsletter of the Division of Gravitational Physics of the American Physical Society, Volume 47, June 2016

abstract: DGRAV News: DGRAV we hear that.... Research Briefs: GW150914 Obituaries: Remembering Felix Pirani Remembering David Finkelstein Steve, the physicist Remembering Sergio Dain Editorial: Gravitational physics in the modern university


title: Learning feed-forward one-shot learners

abstract: One-shot learning is usually tackled by using generative models or discriminative embeddings. Discriminative methods based on deep learning, which are very effective in other learning scenarios, are ill-suited for one-shot learning as they need large amounts of training data. In this paper, we propose a method to learn the parameters of a deep model in one shot. We construct the learner as a second deep network, called a learnet, which predicts the parameters of a pupil network from a single exemplar. In this manner we obtain an efficient feed-forward one-shot learner, trained end-to-end by minimizing a one-shot classification objective in a learning to learn formulation. In order to make the construction feasible, we propose a number of factorizations of the parameters of the pupil network. We demonstrate encouraging results by learning characters from single exemplars in Omniglot, and by tracking visual objects from a single initial exemplar in the Visual Object Tracking benchmark.


title: High-power ultrafast Yb:fiber laser frequency combs using commercially available components and basic fiber tools

abstract: We present a detailed description of the design, construction, and performance of high-power ultrafast Yb:fiber laser frequency combs in operation in our laboratory. We discuss two such laser systems: an 87 MHz, 9 W, 85 fs laser operating at 1060 nm and an 87 MHz, 80 W, 155 fs laser operating at 1035 nm. Both are constructed using low-cost, commercially available components, and can be assembled using only basic tools for cleaving and splicing single-mode fibers. We describe practical methods for achieving and characterizing low-noise single-pulse operation and long-term stability from Yb:fiber oscillators based on nonlinear polarization evolution. Stabilization of the combs using a variety of transducers, including a new method for tuning the carrier-envelope offset frequency, is discussed. High average power is achieved through chirped-pulse amplification in simple fiber amplifiers based on double-clad photonic crystal fibers. We describe the use of these combs in several applications, including ultrasensitive femtosecond time-resolved spectroscopy and cavity-enhanced high-order harmonic generation.


title: Local geodesics between toric plurisubharmonic functions with infinite energy

abstract: Our interest is the behavior of weak geodesics between two plurisubharmonic functions on pseudoconvex domains. We characterize the convergence condition along the geodesic between toric psh functions with a pole at origin on a unit ball in $\mathbb{C}^n$ by means of Lelong numbers.


title: The Schur-Horn theorem for unbounded operators with discrete spectrum

abstract: We characterize diagonals of unbounded self-adjoint operators on a Hilbert space H that have only discrete spectrum, i.e., with empty essential spectrum. Our result extends the Schur-Horn theorem from a finite dimensional setting to an infinite dimensional Hilbert space, analogous to Kadison's theorem for orthogonal projections, Kaftal and Weiss' results for positive compact operators, and Bownik and Jasper's characterization for operators with finite spectrum. Furthermore, we show that if a symmetric unbounded operator E on H has a nondecreasing unbounded diagonal, then any sequence that weakly majorizes this diagonal is also a diagonal of E.


title: On well-posedness of vector-valued fractional differential-difference equations

abstract: We develop an operator-theoretical method for the analysis on well posedness of partial differential equations that can be modeled in the form \begin{equation*} \left\{ \begin{array}{rll} \Delta^{\alpha} u(n) &= Au(n+2) + f(n,u(n)), \quad n \in \mathbb{N}_0, \,\, 1< \alpha \leq 2; u(0) &= u_0; u(1) &= u_1, \end{array} \right. \end{equation*} where $A$ is an closed linear operator defined on a Banach space $X$. Our ideas are inspired on the Poisson distribution as a tool to sampling fractional differential operators into fractional differences. Using our abstract approach, we are able to show existence and uniqueness of solutions for the problem (*) on a distinguished class of weighted Lebesgue spaces of sequences, under mild conditions on strongly continuous sequences of bounded operators generated by $A,$ and natural restrictions on the nonlinearity $f$. Finally we present some original examples to illustrate our results.


title: Estimation and Uncertainty Quantification of Yield Via Strain Recovery Simulations

abstract: In computational materials science, predicting the yield strain of crosslinked polymers remains a challenging task. A common approach is to identify yield via the first critical point of stress-strain curves produced by molecular dynamics simulations. However, the simulated data can be excessively noisy, making it difficult to extract meaningful results. In this work, we discuss an alternate method for identifying yield on the basis of residual strain computations. Notably, the associated raw data produce a sharper signal for yield through a transition in their global behavior. As we show, this transition can be analyzed in terms of simple functions (e.g. hyperbolas) that admit straightforward uncertainty quantification techniques.


title: Dynamical Contents of Unconventional Supersymmetry

abstract: The Dirac Hamiltonian formalism is applied to a system in $(2+1)$-dimensions consisting of a Dirac field $\psi$ minimally coupled to Chern-Simons $U(1)$ and $SO(2,1)$ connections, $A$ and $\omega$, respectively. This theory is connected to a supersymmetric Chern-Simons form in which the gravitino has been projected out (unconventional supersymmetry) and, in the case of a flat background, corresponds to the low energy limit of graphene. The separation between first-class and second-class constraints is performed explicitly, and both the field equations and gauge symmetries of the Lagrangian formalism are fully recovered. The degrees of freedom of the theory in generic sectors shows that the propagating states correspond to fermionic modes in the background determined by the geometry of the graphene sheet and the nondynamical electromagnetic field. This is shown for the following canonical sectors: i) a conformally invariant generic description where the spinor field and the dreibein are locally rescaled; ii) a specific configuration for the Dirac fermion consistent with its spin, where Weyl symmetry is exchanged by time reparametrizations; iii) the vacuum sector $\psi=0$, which is of interest for perturbation theory. For the latter the analysis is adapted to the case of manifolds with boundary, and the corresponding Dirac brackets together with the centrally extended charge algebra are found. Finally, the $SU(2)$ generalization of the gauge group is briefly treated, yielding analogous conclusions for the degrees of freedom.


title: A linear time algorithm for a variant of the max cut problem in series parallel graphs

abstract: Given a graph $G=(V, E)$, a connected sides cut $(U, V\backslash U)$ or $\delta (U)$ is the set of edges of E linking all vertices of U to all vertices of $V\backslash U$ such that the induced subgraphs $G[U]$ and $G[V\backslash U]$ are connected. Given a positive weight function $w$ defined on $E$, the maximum connected sides cut problem (MAX CS CUT) is to find a connected sides cut $\Omega$ such that $w(\Omega)$ is maximum. MAX CS CUT is NP-hard. In this paper, we give a linear time algorithm to solve MAX CS CUT for series parallel graphs. We deduce a linear time algorithm for the minimum cut problem in the same class of graphs without computing the maximum flow.


title: The Mondrian Kernel

abstract: We introduce the Mondrian kernel, a fast random feature approximation to the Laplace kernel. It is suitable for both batch and online learning, and admits a fast kernel-width-selection procedure as the random features can be re-used efficiently for all kernel widths. The features are constructed by sampling trees via a Mondrian process [Roy and Teh, 2009], and we highlight the connection to Mondrian forests [Lakshminarayanan et al., 2014], where trees are also sampled via a Mondrian process, but fit independently. This link provides a new insight into the relationship between kernel methods and random forests.


title: Sentiment Aggregate Functions for Political Opinion Polling using Microblog Streams

abstract: The automatic content analysis of mass media in the social sciences has become necessary and possible with the raise of social media and computational power. One particularly promising avenue of research concerns the use of sentiment analysis in microblog streams. However, one of the main challenges consists in aggregating sentiment polarity in a timely fashion that can be fed to the prediction method. We investigated a large set of sentiment aggregate functions and performed a regression analysis using political opinion polls as gold standard. Our dataset contains nearly 233 000 tweets, classified according to their polarity (positive, negative or neutral), regarding the five main Portuguese political leaders during the Portuguese bailout (2011-2014). Results show that different sentiment aggregate functions exhibit different feature importance over time while the error keeps almost unchanged.


title: Orthogonality of divisorial Zariski decompositions for classes with volume zero

abstract: We show that the orthogonality conjecture for divisorial Zariski decompositions on compact Kahler manifolds holds for pseudoeffective (1,1) classes with volume zero.


title: Notes on heat engines and negative temperatures

abstract: We show that a Carnot cycle operating between a positive canonical-temperature bath and a negative canonical-temperature bath has efficiency equal to unity. It follows that a negative canonical-temperature cannot be identified with an absolute temperature. We illustrate this with a spin in a varying magnetic field.


title: Networked Control under Random and Malicious Packet Losses

abstract: We study cyber security issues in networked control of a linear dynamical system. Specifically, the dynamical system and the controller are assumed to be connected through a communication channel that face malicious attacks as well as random packet losses due to unreliability of transmissions. We provide a probabilistic characterization for the link failures which allows us to study combined effects of malicious and random packet losses. We first investigate almost sure stabilization under an event-triggered control law, where we utilize Lyapunov-like functions to characterize the triggering times at which the plant and the controller attempt to exchange state and control data over the network. We then provide a look at the networked control problem from the attacker's perspective and explore malicious attacks that cause instability. Finally, we demonstrate the efficacy of our results with numerical examples.


title: Development of n-in-p pixel modules for the ATLAS Upgrade at HL-LHC

abstract: Thin planar pixel modules are promising candidates to instrument the inner layers of the new ATLAS pixel detector for HL-LHC, thanks to the reduced contribution to the material budget and their high charge collection efficiency after irradiation. 100-200 $\mu$m thick sensors, interconnected to FE-I4 read-out chips, have been characterized with radioactive sources and beam tests at the CERN-SPS and DESY. The results of these measurements are reported for devices before and after irradiation up to a fluence of $14\times10^{15}$ n$_{eq}$/cm$^2$. The charge collection and tracking efficiency of the different sensor thicknesses are compared. The outlook for future planar pixel sensor production is discussed, with a focus on sensor design with the pixel pitches (50x50 and 25x100 $\mu$m$^2$) foreseen for the RD53 Collaboration read-out chip in 65 nm CMOS technology. An optimization of the biasing structures in the pixel cells is required to avoid the hit efficiency loss presently observed in the punch-through region after irradiation. For this purpose the performance of different layouts have been compared in FE-I4 compatible sensors at various fluence levels by using beam test data. Highly segmented sensors will represent a challenge for the tracking in the forward region of the pixel system at HL-LHC. In order to reproduce the performance of 50x50 $\mu$m$^2$ pixels at high pseudo-rapidity values, FE-I4 compatible planar pixel sensors have been studied before and after irradiation in beam tests at high incidence angle (80$^\circ$) with respect to the short pixel direction. Results on cluster shapes, charge collection and hit efficiency will be shown.


title: Search for giant planets in M67 III: excess of hot Jupiters in dense open clusters

abstract: Since 2008 we used high-precision radial velocity (RV) measurements obtained with different telescopes to detect signatures of massive planets around main-sequence and evolved stars of the open cluster (OC) M67. We aimed to perform a long-term study on giant planet formation in open clusters and determine how this formation depends on stellar mass and chemical composition. A new hot Jupiter (HJ) around the main-sequence star YBP401 is reported in this work. An update of the RV measurements for the two HJ host-stars YBP1194 and YBP1514 is also discussed. Our sample of 66 main-sequence and turnoff stars includes 3 HJs, which indicates a high rate of HJs in this cluster (~5.6% for single stars and ~4.5% for the full sample ). This rate is much higher than what has been discovered in the field, either with RV surveys or by transits. High metallicity is not a cause for the excess of HJs in M67, nor can the excess be attributed to high stellar masses. When combining this rate with the non-zero eccentricity of the orbits, our results are qualitatively consistent with a HJ formation scenario dominated by strong encounters with other stars or binary companions and subsequent planet-planet scattering, as predicted by N-body simulations.


title: Leadership Network and Team Performance in Interactive Contests

abstract: Over the years, the concept of leadership has experienced a paradigm shift - from solitary leader (centralized leadership) to de-centralized leadership or distributed leadership. This paper explores the idea that centralized leadership, as earlier suggested, negatively impacts team performance. I applied the hypothesis to cricket, a sport in which leaders play an important role in team's success. I generated batting partnership network and evaluated the central-most player in the team, applying tools of social network analysis. Analyzing 3420 matches in one day international cricket and 1979 Test matches involving 10 teams, I examined the impact of centralized leadership in outcome of a contest. I observed that the odds for winning a one day international match under centralized leadership is 30% higher than the odds for winning under de-centralized leadership. In both forms of cricket (Test and one day international ), I failed to find evidence that distributed leadership is associated with higher team performance. These results suggest important implications for cricket administrators in development and management of working teams.


title: The Nature and Frequency of the Gas Outbursts in Comet 67P/Churyumov-Gerasimenko observed by the Alice Far-ultraviolet Spectrograph on Rosetta

abstract: Alice is a far-ultraviolet imaging spectrograph onboard Rosetta that, amongst multiple objectives, is designed to observe emissions from various atomic and molecular species from within the coma of comet 67P/Churyumov-Gerasimenko. The initial observations, made following orbit insertion in August 2014, showed emissions of atomic hydrogen and oxygen spatially localized close to the nucleus and attributed to photoelectron impact dissociation of H2O vapor. Weaker emissions from atomic carbon were subsequently detected and also attributed to electron impact dissociation, of CO2, the relative H I and C I line intensities reflecting the variation of CO2 to H2O column abundance along the line-of-sight through the coma. Beginning in mid-April 2015, Alice sporadically observed a number of outbursts above the sunward limb characterized by sudden increases in the atomic emissions, particularly the semi-forbidden O I 1356 multiplet, over a period of 10-30 minutes, without a corresponding enhancement in long wavelength solar reflected light characteristic of dust production. A large increase in the brightness ratio O I 1356/O I 1304 suggests O2 as the principal source of the additional gas. These outbursts do not correlate with any of the visible images of outbursts taken with either OSIRIS or the navigation camera. Beginning in June 2015 the nature of the Alice spectrum changed considerably with CO Fourth Positive band emission observed continuously, varying with pointing but otherwise fairly constant in time. However, CO does not appear to be a major driver of any of the observed outbursts.


title: SQuAD: 100,000+ Questions for Machine Comprehension of Text

abstract: We present the Stanford Question Answering Dataset (SQuAD), a new reading comprehension dataset consisting of 100,000+ questions posed by crowdworkers on a set of Wikipedia articles, where the answer to each question is a segment of text from the corresponding reading passage. We analyze the dataset to understand the types of reasoning required to answer the questions, leaning heavily on dependency and constituency trees. We build a strong logistic regression model, which achieves an F1 score of 51.0%, a significant improvement over a simple baseline (20%). However, human performance (86.8%) is much higher, indicating that the dataset presents a good challenge problem for future research. The dataset is freely available at https://stanford-qa.com


title: On The Fourier And Wavelet Analysis Of Coronal Time Series

abstract: Using Fourier and wavelet analysis, we critically re-assess the significance of our detection of periodic pulsations in coronal loops. We show that the proper identification of the frequency dependence and statistical properties of the different components of the power spectra provies a strong argument against the common practice of data detrending, which tends to produce spurious detections around the cut-off frequency of the filter. In addition, the white and red noise models built into the widely used wavelet code of Torrence & Compo cannot, in most cases, adequately represent the power spectra of coronal time series, thus also possibly causing false positives. Both effects suggest that several reports of periodic phenomena should be re-examined. The Torrence & Compo code nonetheless effectively computes rigorous confidence levels if provided with pertinent models of mean power spectra, and we describe the appropriate manner in which to call its core routines. We recall the meaning of the default confidence levels output from the code, and we propose new Monte-Carlo-derived levels that take into account the total number of degrees of freedom in the wavelet spectra. These improvements allow us to confirm that the power peaks that we detected have a very low probability of being caused by noise.


title: Division by 2 on hyperelliptic curves and jacobians

abstract: Let $K$ be an algebraically closed field of characteristic different from 2, $g$ a positive integer, $f(x)$ a degree $(2g+1)$ polynomial with coefficients in $K$ and without multiple roots, $C: y^2=f(x)$ the corresponding genus $g$ hyperelliptic curve over $K$ and $J$ the jacobian of $C$. We identify $C$ with the image of its canonical embedding into $J$ (the infinite point of $C$ goes to the zero point of $J$). For each point $P=(a,b)\in C(K)$ there are $2^{2g}$ points $\frac{1}{2}P \in J(K)$. We describe explicitly the Mumford represesentations of all $\frac{1}{2}P$. The rationality questions for $\frac{1}{2}P$ are also discussed.


title: Universal criterion for designability of heteropolymers

abstract: Proteins are an example of heteropolymers able to self-assemble in specific target structures. The self-assembly of designed artificial heteropolymers is still, to the best of our knowledge, not possible with control over the single chain self-assembling properties comparable to what natural proteins can achieve. What artificial heteropolymers lacks compared to bio-heteropolymers that grants the latter such a versatility? Is the geometry of the protein skeleton the only a particular choice to be designable? Here we introduce a general criteria to discriminate which polymer backbones can be designed to adopt a predetermined structure. With our approach we can explore different polymer backbones and different amino acids alphabets. By comparing the radial distribution functions of designable and not-designable scenarios we identify as designability criteria the presence of a particular peak in the radial distribution function that dominates over the random packing of the heteropolymer. We show that the peak is a universal feature of all designable heteropolymers, as it is dominating also the radial distribution function of natural proteins. Our finding can help in understanding the key features that make proteins a highly designable system. The criteria that we present can be applied to engineer new types of self-assembling modular polymers that will open new applications for polymer based material science.


title: Recent advancement in Next Generation Sequencing techniques and its computational analysis

abstract: Next Generation Sequencing (NGS), a recently evolved technology, have served a lot in the research and development sector of our society. This novel approach is a newbie and has critical advantages over the traditional Capillary Electrophoresis (CE) based Sanger Sequencing. The advancement of NGS has led to numerous important discoveries, which could have been costlier and time taking in case of traditional CE based Sanger sequencing. NGS methods are highly parallelized enabling to sequence thousands to millions of molecules simultaneously. This technology results into huge amount of data, which need to be analysed to conclude valuable information. Specific data analysis algorithms are written for specific task to be performed. The algorithms in group, act as a tool in analysing the NGS data. Analysis of NGS data unravels important clues in quest for the treatment of various life-threatening diseases; improved crop varieties and other related scientific problems related to human welfare. In this review, an effort was made to address basic background of NGS technologies, possible applications, computational approaches and tools involved in NGS data analysis, future opportunities and challenges in the area.


title: 3D zigzag for multislicing, multiband and video processing

abstract: We present a 3D zigzag rafter (first in literature) which allows us to obtain the exact sequence of spectral components after application of Discrete Cosine Transform 3D (DCT-2D) over a cube. Such cube represents part of a video or eventually a group of images such as multislicing (e.g., Magnetic Resonance or Computed Tomography imaging) and multi or hyperspectral imagery (optical satellites). Besides, we present a new version of the traditional 2D zigzag, including the case of rectangular blocks. Finally, all the attached code is done in MATLAB, and that code serves both blocks of pixels or blocks of blocks.


title: Effects of Spatial Heterogeneity in Rainfall and Vegetation Type on Soil Moisture and Evapotranspiration

abstract: Nonlinear plant-scale interactions controlling the soil-water balance are generally not valid at larger spatial scales due to spatial heterogeneity in rainfall and vegetation type. The relationships between spatially averaged variables are hysteretic even when unique relationships are imposed at the plant scale. The characteristics of these hysteretic relationships depend on the size of the averaging area and the spatial properties of the soil, vegetation, and rainfall. We upscale the plant-scale relationships to the scale of a regional land-surface model based on simulation data obtained through explicit representation of spatial heterogeneity in rainfall and vegetation type. The proposed upscaled function improves predictions of spatially averaged soil moisture and evapotranspiration relative to the effective-parameter approach for a water-limited Texas shrubland. The degree of improvement is a function of the scales of heterogeneity and the size of the averaging area. We also find that single-valued functions fail to predict spatially averaged leakage accurately. Furthermore, the spatial heterogeneity results in scale-dependent hysteretic relationships for the statistical-dynamic and Montaldo & Albertson approaches.


title: DNA brick self-assembly with an off-lattice potential

abstract: We report Monte Carlo simulations of a simple off-lattice patchy-particle model for DNA `bricks'. We relate the parameters that characterise this model with the binding free energy of pairs of single-stranded DNA molecules. We verify that an off-lattice potential parameterised in this way reproduces much of the behaviour seen with a simpler lattice model we introduced previously, although the relaxation of the geometric constraints leads to a more error-prone self-assembly pathway. We investigate the self-assembly process as a function of the strength of the non-specific interactions. We show that our off-lattice model for DNA bricks results in robust self-assembly into a variety of target structures.


title: Magnitudes and Timescales of Total Solar Irradiance Variability

abstract: The Sun's net radiative output varies on timescales of minutes to gigayears. Direct measurements of the total solar irradiance (TSI) show changes in the spatially- and spectrally-integrated radiant energy on timescales as short as minutes to as long as a solar cycle. Variations of ~0.01 % over a few minutes are caused by the ever-present superposition of convection and oscillations with very large solar flares on rare occasion causing slightly-larger measureable signals. On timescales of days to weeks, changing photospheric magnetic activity affects solar brightness at the ~0.1 % level. The 11-year solar cycle shows variations of comparable magnitude with irradiances peaking near solar maximum. Secular variations are more difficult to discern, being limited by instrument stability and the relatively short duration of the space-borne record. Historical reconstructions of the Sun's irradiance based on indicators of solar-surface magnetic activity, such as sunspots, faculae, and cosmogenic isotope records, suggest solar brightness changes over decades to millennia, although the magnitudes of these variations have high uncertainties due to the indirect historical records on which they rely. Stellar evolution affects yet longer timescales and is responsible for the greatest solar variabilities. In this manuscript I summarize the Sun's variability magnitudes over different temporal regimes and discuss the irradiance record's relevance for solar and climate studies as well as for detections of exo-solar planets transiting Sun-like stars.