title: Testing the threshold expansion for three-particle energies at fourth order in $\phi^4$ theory

abstract: A relativistic formalism for relating the the energies of the states of three scalar particles in finite volume to infinite volume scattering amplitudes has recently been developed. This formalism has been used to predict the energy of the state closest to threshold in an expansion in powers of $1/L$, with $L$ the box length. This expansion has been tested previously by a perturbative calculation of the threshold energy in $\lambda \phi^4$ theory, working to third order in $\lambda$ and up to $\mathcal O(1/L^6)$ in the volume expansion. However, several aspects of the predicted threshold behavior do not enter until fourth (three-loop) order in perturbation theory. Here I extend the perturbative calculation to fourth order and find agreement with the general prediction. This check also requires a two-loop calculation of the infinite-volume off-shell two-particle scattering amplitude near threshold. As a spin-off, I check the threshold expansion for two particles to the same order, finding agreement with the result that follows from L\"uscher's formalism.


title: GeV detection of HESS J0632+057

abstract: HESS J0632+057 is the only gamma-ray binary that has been detected at TeV energies, but not at GeV energies yet. Based on nearly nine years of Fermi Large Area Telescope (LAT) Pass 8 data, we report here on a deep search for the gamma-ray emission from HESS J0632+057 in the 0.1-300 GeV energy range. We find a previously unknown gamma-ray source, Fermi J0632.6+0548, spatially coincident with HESS J0632+057. The measured flux of Fermi J0632.6+0548 is consistent with the previous flux upper limit on HESS J0632+057 and shows variability that can be related to the HESS J0632+057 orbital phase. We propose that Fermi J0632.6+0548 is the GeV counterpart of HESS J0632+057. Considering the Very High Energy (VHE) spectrum of HESS J0632+057, a possible spectral turnover above 10 GeV may exist in Fermi J0632.6+0548, as appears to be common in other established gamma-ray binaries.


title: Exploring Dimensionality Reductions with Forward and Backward Projections

abstract: Dimensionality reduction is a common method for analyzing and visualizing high-dimensional data across domains. Dimensionality-reduction algorithms involve complex optimizations and the reduced dimensions computed by these algorithms generally lack clear relation to the initial data dimensions. Therefore, interpreting and reasoning about dimensionality reductions can be difficult. In this work, we introduce two interaction techniques, \textit{forward projection} and \textit{backward projection}, for reasoning dynamically about scatter plots of dimensionally reduced data. We also contribute two related visualization techniques, \textit{prolines} and \textit{feasibility map} to facilitate and enrich the effective use of the proposed interactions, which we integrate in a new tool called \textit{Praxis}. To evaluate our techniques, we first analyze their time and accuracy performance across varying sample and dimension sizes. We then conduct a user study in which twelve data scientists use \textit{Praxis} so as to assess the usefulness of the techniques in performing exploratory data analysis tasks. Results suggest that our visual interactions are intuitive and effective for exploring dimensionality reductions and generating hypotheses about the underlying data.


title: Polynomial Counting in Anonymous Dynamic Networks with Applications to Anonymous Dynamic Algebraic Computations

abstract: Starting with Michail, Chatzigiannakis, and Spirakis work, the problem of Counting the number of nodes in Anonymous Dynamic Networks has attracted a lot of attention. The problem is challenging because nodes are indistinguishable (they lack identifiers and execute the same program) and the topology may change arbitrarily from round to round of communication, as long as the network is connected in each round. The problem is central in distributed computing as the number of participants is frequently needed to make important decisions, such as termination, agreement, synchronization, and many others. A variety of algorithms built on top of mass-distribution techniques have been presented, analyzed, and also experimentally evaluated; some of them assumed additional knowledge of network characteristics, such as bounded degree or given upper bound on the network size. However, the question of whether Counting can be solved deterministically in sub-exponential time remained open. In this work, we answer this question positively by presenting Methodical Counting, which runs in polynomial time and requires no knowledge of network characteristics. Moreover, we also show how to extend Methodical Counting to compute the sum of input values and more complex functions without extra cost. Our analysis leverages previous work on random walks in evolving graphs, combined with carefully chosen alarms in the algorithm that control the process and its parameters. To the best of our knowledge, our Counting algorithm and its extensions to other algebraic and Boolean functions are the first that can be implemented in practice with worst-case guarantees.


title: A mathematical model for DNA

abstract: Recently, some authors have shown that a DNA molecule produces electromagnetic signals and communicates with other DNA molecules or other molecules. In fact, a DNA acts like a receiver or transmitter of radio waves. In this paper, we suggest a mathematical model for the DNA molecule and use of its communication to cure some diseases like cancer. In this model, first, by using concepts from string theory and M-theory, we calculate the energy of a DNA in terms of interactions between free electrons and bound electrons. We show that when a DNA is damaged, its energy changes and an extra current is produced. This extra current causes the electromagnetic signals of a damaged DNA molecule to be different when compared to the electromagnetic signals of a normal DNA molecule. The electromagnetic signals of a damaged DNA molecule induces an extra current in a normal DNA molecule and leads to its destruction. By sending crafted electromagnetic signals to normal DNA molecules and inducing an opposite current with respect to this extra current, we can prevent the destruction of normal DNA. Finally, we argue that the type of packing of DNA in chromosomes of men and women are different. This causes radiated waves from DNAs of men and women to have opposite signs and cancel the effect of each other in a pair. Using this property, we suggest another mechanism to cancel the effect of extra waves, which are produced by DNAs in cancer cells of a male or a female, by extra waves which are produced by DNAs in similar cells of a female or a male and prevent the progression of the disease.


title: Who Will You Share a Ride With: Factors that Influence Trust of Potential Rideshare Partners

abstract: Ridesharing can save energy, reduce pollution, and alleviate congestion. Ridesharing systems allow peer-to-peer ride matching through mobile-phone applications. In many systems people have little information about the driver and co-riders. Trust heavily influences the usage of a ridesharing system, however, work in this area is limited. In this paper, we investigate the role of social media in shaping trust and delineate factors that influence the decisions of trust placed on a driver or co-rider. We design a ridesharing system where a rider can see the Instagram profile and 10 recent pictures of a co-rider and ask them to make judgments on six trust related questions. We discovered three important factors that shape trust forming: \textit{social proof}, \textit{social approval}, and \textit{self-disclosure}. These three factors account for 40\%, 33\%, and 13\% of the variance, and subsequent factors extracted accounted for substantially less variance. Using these three factors, we successfully build a trust prediction model with F-measure of 87.4\%, an improvement from 78.6\% on using all factors collected.


title: Zipf's Law for Atlas Models

abstract: A set of data with positive values follows a Pareto distribution if the log-log plot of value versus rank is approximately a straight line. A Pareto distribution satisfies Zipf's law if the log-log plot has a slope of -1. Since many types of ranked data follow Zipf's law, it is considered a form of universality. We propose a mathematical explanation for this phenomenon based on Atlas models and first-order models, systems of positive continuous semimartingales with parameters that depend only on rank. We show that the stable distribution of an Atlas model will follow Zipf's law if and only if two natural conditions, conservation and completeness, are satisfied. Since Atlas models and first-order models can be constructed to approximate systems of time-dependent rank-based data, our results can explain the universality of Zipf's law for such systems. However, ranked data generated by other means may follow non-Zipfian Pareto distributions. Hence, our results explain why Zipf's law holds for word frequency, firm size, household wealth, and city size, while it does not hold for earthquake magnitude, cumulative book sales, the intensity of solar flares, and the intensity of wars, all of which follow non-Zipfian Pareto distributions.


title: Sign-change of the Fourier coefficients of a Hauptmodul for $\Gamma_{0}(2)$

abstract: In this short note, we aim to prove that the Fourier coefficients of the modular function $ \frac{\eta^{24}(\tau)}{\eta^{24}(2\tau)} $ possess a sign-change property.


title: Towards the ab initio based theory of the phase transformations in iron and steel

abstract: Despite of the appearance of numerous new materials, the iron based alloys and steels continue to play an essential role in modern technology. The properties of a steel are determined by its structural state (ferrite, cementite, pearlite, bainite, martensite, and their combination) that is formed under thermal treatment as a result of the shear lattice reconstruction "gamma" (fcc) -> "alpha" (bcc) and carbon diffusion redistribution. We present a review on a recent progress in the development of a quantitative theory of the phase transformations and microstructure formation in steel that is based on an ab initio parameterization of the Ginzburg-Landau free energy functional. The results of computer modeling describe the regular change of transformation scenario under cooling from ferritic (nucleation and diffusion-controlled growth of the "alpha" phase to martensitic (the shear lattice instability "gamma" -> "alpha"). It has been shown that the increase in short-range magnetic order with decreasing the temperature plays a key role in the change of transformation scenarios. Phase-field modeling in the framework of a discussed approach demonstrates the typical transformation patterns.


title: Shortcut to adiabaticity in a Stern-Gerlach apparatus

abstract: We show that the performances of a Stern-Gerlach apparatus can be improved by using a magnetic field profile for the atomic spin evolution designed through shortcut to adiabaticity technique. Interestingly, it can be made more compact - for atomic beams propagating at a given velocity - and more resilient to a dispersion in velocity, in comparison with the results obtained with a standard uniform rotation of the magnetic field. Our results are obtained using a reverse engineering approach based on Lewis-Riesenfeld invariants. We discuss quantitatively the advantages offered by our configuration in terms of the resources involved and show that it drastically enhances the fidelity of the quantum state transfer achieved by the Stern-Gerlach device.


title: Band structure of a IV-VI black phosphorus analogue, the thermoelectric SnSe

abstract: The success of black phosphorus in fast electronic and photonic devices is hindered by its rapid degradation in presence of oxygen. Orthorhombic tin selenide is a representative of group IV-VI binary compounds that are robust, isoelectronic, and share the same structure with black phosphorus. We measured the band structure of SnSe and found highly anisotropic valence bands that form several valleys having fast dispersion within the layers and negligible dispersion across. This is exactly the band structure desired for efficient thermoelectric generation where SnSe has shown a great promise.


title: Numerical Approximation of the Integral Fractional Laplacian

abstract: We propose a new nonconforming finite element algorithm to approximate the solution to the elliptic problem involving the fractional Laplacian. We first derive an integral representation of the bilinear form corresponding to the variational problem. The numerical approximation of the action of the corresponding stiffness matrix consists of three steps: (i) apply a sinc quadrature scheme to approximate the integral representation by a finite sum where each term involves the solution of an elliptic partial differential equation defined on the entire space, (ii) truncate each elliptic problem to a bounded domain, (iii) use the finite element method for the space approximation on each truncated domain. The consistency error analysis for the three steps is discussed together with the numerical implementation of the entire algorithm. The results of computations are given illustrating the error behavior in terms of the mesh size of the physical domain, the domain truncation parameter and the quadrature spacing parameter.


title: Predicting Abandonment in Online Coding Tutorials

abstract: Learners regularly abandon online coding tutorials when they get bored or frustrated, but there are few techniques for anticipating this abandonment to intervene. In this paper, we examine the feasibility of predicting abandonment with machine-learned classifiers. Using interaction logs from an online programming game, we extracted a collection of features that are potentially related to learner abandonment and engagement, then developed classifiers for each level. Across the first five levels of the game, our classifiers successfully predicted 61% to 76% of learners who did not complete the next level, achieving an average AUC of 0.68. In these classifiers, features negatively associated with abandonment included account activation and help-seeking behaviors, whereas features positively associated with abandonment included features indicating difficulty and disengagement. These findings highlight the feasibility of providing timely intervention to learners likely to quit.


title: Characterization of the K2-18 multi-planetary system with HARPS: A habitable zone super-Earth and discovery of a second, warm super-Earth on a non-coplanar orbit

abstract: The bright M dwarf K2-18 at 34 pc is known to host a transiting super-Earth-sized planet orbiting within the star's habitable zone; K2-18b. Given the superlative nature of this system for studying an exoplanetary atmosphere receiving similar levels of insolation as the Earth, we aim to characterize the planet's mass which is required to interpret atmospheric properties and infer the planet's bulk composition. We obtain precision radial velocity measurements with the HARPS spectrograph and couple those measurements with the K2 photometry to jointly model the observed radial velocity variation with planetary signals and a radial velocity jitter model based on Gaussian process regression. We measure the mass of K2-18b to be $8.0 \pm 1.9$ M$_{\oplus}$ with a bulk density of $3.7 \pm 0.9$ g/cm$^3$ which may correspond to a predominantly rocky planet with a significant gaseous envelope or an ocean planet with a water mass fraction $\gtrsim 50$%. We also find strong evidence for a second, warm super-Earth K2-18c at $\sim 9$ days with a semi-major axis 2.4 times smaller than the transiting K2-18b. After re-analyzing the available light curves of K2-18 we conclude that K2-18c is not detected in transit and therefore likely has an orbit that is non-coplanar with K2-18b. A suite of dynamical integrations with varying simulated orbital eccentricities of the two planets are used to further constrain each planet's eccentricity posterior from which we measure $e_b < 0.43$ and $e_c < 0.47$ at 99% confidence. The discovery of the inner planet K2-18c further emphasizes the prevalence of multi-planet systems around M dwarfs. The characterization of the density of K2-18b reveals that the planet likely has a thick gaseous envelope which along with its proximity to the Solar system makes the K2-18 planetary system an interesting target for the atmospheric study of an exoplanet receiving Earth-like insolation.


title: A short introduction to quasi-Monte Carlo option pricing

abstract: One of the main practical applications of quasi-Monte Carlo (QMC) methods is the valuation of financial derivatives. We aim to give a short introduction into option pricing and show how it is facilitated using QMC. We give some practical examples for illustration.


title: Ergodic Coverage In Constrained Environments Using Stochastic Trajectory Optimization

abstract: In search and surveillance applications in robotics, it is intuitive to spatially distribute robot trajectories with respect to the probability of locating targets in the domain. Ergodic coverage is one such approach to trajectory planning in which a robot is directed such that the percentage of time spent in a region is in proportion to the probability of locating targets in that region. In this work, we extend the ergodic coverage algorithm to robots operating in constrained environments and present a formulation that can capture sensor footprint and avoid obstacles and restricted areas in the domain. We demonstrate that our formulation easily extends to coordination of multiple robots equipped with different sensing capabilities to perform ergodic coverage of a domain.


title: Approximation Schemes for Clustering with Outliers

abstract: Clustering problems are well-studied in a variety of fields such as data science, operations research, and computer science. Such problems include variants of centre location problems, $k$-median, and $k$-means to name a few. In some cases, not all data points need to be clustered; some may be discarded for various reasons. We study clustering problems with outliers. More specifically, we look at Uncapacitated Facility Location (UFL), $k$-Median, and $k$-Means. In UFL with outliers, we have to open some centres, discard up to $z$ points of $\cal X$ and assign every other point to the nearest open centre, minimizing the total assignment cost plus centre opening costs. In $k$-Median and $k$-Means, we have to open up to $k$ centres but there are no opening costs. In $k$-Means, the cost of assigning $j$ to $i$ is $\delta^2(j,i)$. We present several results. Our main focus is on cases where $\delta$ is a doubling metric or is the shortest path metrics of graphs from a minor-closed family of graphs. For uniform-cost UFL with outliers on such metrics we show that a multiswap simple local search heuristic yields a PTAS. With a bit more work, we extend this to bicriteria approximations for the $k$-Median and $k$-Means problems in the same metrics where, for any constant $\epsilon > 0$, we can find a solution using $(1+\epsilon)k$ centres whose cost is at most a $(1+\epsilon)$-factor of the optimum and uses at most $z$ outliers. We also show that natural local search heuristics that do not violate the number of clusters and outliers for $k$-Median (or $k$-Means) will have unbounded gap even in Euclidean metrics. Furthermore, we show how our analysis can be extended to general metrics for $k$-Means with outliers to obtain a $(25+\epsilon,1+\epsilon)$ bicriteria.


title: Self-Gravity in Magnetized Neutrino-Dominated Accretion Discs

abstract: What is conducted in the present work is the study of self-gravity effects on vertical structure of magnetized neutrino-dominated accretion disc (termed NDAF and considered as a central engine for Gamma-ray bursts (GRBs)). Some of the disc physical time scales such as viscous, cooling and diffusion time scales, that are supposed to play a pivotal role in the late time evolutions of the disc, have been studied. It is of our interest to investigate the possibility of the X-ray flares occurrence, observed in late time GRB's extended emission through the "magnetic barrier" and "fragmentation" processes in our model. The results lead us to interpret self-gravity as an amplifier for Blandford-Payne luminosity (BP power) and the generated magnetic field, but a suppressor for neutrino luminosity and magnetic barrier process via highlighting fragmentation mechanism in the outer disc, especially for the higher mass accretion rates.


title: The size-Ramsey number of powers of paths

abstract: Given graphs $G$ and $H$ and a positive integer $q$ say that $G$ is $q$-Ramsey for $H$, denoted $G\rightarrow (H)_q$, if every $q$-colouring of the edges of $G$ contains a monochromatic copy of $H$. The size-Ramsey number $\hat{r}(H)$ of a graph $H$ is defined to be $\hat{r}(H)=\min\{|E(G)|\colon G\rightarrow (H)_2\}$. Answering a question of Conlon, we prove that, for every fixed $k$, we have $\hat{r}(P_n^k)=O(n)$, where $P_n^k$ is the $k$-th power of the $n$-vertex path $P_n$ (i.e. , the graph with vertex set $V(P_n)$ and all edges $\{u,v\}$ such that the distance between $u$ and $v$ in $P_n$ is at most $k$). Our proof is probabilistic, but can also be made constructive.


title: Strategic Coalitions with Perfect Recall

abstract: The paper proposes a bimodal logic that describes an interplay between distributed knowledge modality and coalition know-how modality. Unlike other similar systems, the one proposed here assumes perfect recall by all agents. Perfect recall is captured in the system by a single axiom. The main technical results are the soundness and the completeness theorems for the proposed logical system.


title: A grid of Synthetic Spectra for Hot DA White Dwarfs and Its Application in Stellar Population Synthesis

abstract: In this work we present a grid of LTE and non-LTE synthetic spectra of hot DA white dwarfs (WDs). In addition to its usefulness for the determination of fundamental stellar parameters of isolated WDs and in binaries, this grid will be of interest for the construction of theoretical libraries for stellar studies from integrated light. The spectral grid covers both a wide temperature and gravity range, with 17,000 K <= T_eff <= 100,000 K and 7.0 <= log(g) <= 9.5. The stellar models are built for pure hydrogen and the spectra cover a wavelength range from 900 A to 2.5 microns. Additionally, we derive synthetic HST/ACS, HST/WFC3, Bessel UBVRI and SDSS magnitudes. The grid was also used to model integrated spectral energy distributions of simple stellar populations and our modeling suggests that DAs might be detectable in ultraviolet bands for populations older than ~8 Gyr.


title: Coalescent-based species tree estimation: a stochastic Farris transform

abstract: The reconstruction of a species phylogeny from genomic data faces two significant hurdles: 1) the trees describing the evolution of each individual gene--i.e., the gene trees--may differ from the species phylogeny and 2) the molecular sequences corresponding to each gene often provide limited information about the gene trees themselves. In this paper we consider an approach to species tree reconstruction that addresses both these hurdles. Specifically, we propose an algorithm for phylogeny reconstruction under the multispecies coalescent model with a standard model of site substitution. The multispecies coalescent is commonly used to model gene tree discordance due to incomplete lineage sorting, a well-studied population-genetic effect. In previous work, an information-theoretic trade-off was derived in this context between the number of loci, $m$, needed for an accurate reconstruction and the length of the locus sequences, $k$. It was shown that to reconstruct an internal branch of length $f$, one needs $m$ to be of the order of $1/[f^{2} \sqrt{k}]$. That previous result was obtained under the molecular clock assumption, i.e., under the assumption that mutation rates (as well as population sizes) are constant across the species phylogeny. Here we generalize this result beyond the restrictive molecular clock assumption, and obtain a new reconstruction algorithm that has the same data requirement (up to log factors). Our main contribution is a novel reduction to the molecular clock case under the multispecies coalescent. As a corollary, we also obtain a new identifiability result of independent interest: for any species tree with $n \geq 3$ species, the rooted species tree can be identified from the distribution of its unrooted weighted gene trees even in the absence of a molecular clock.


title: Mellin-Meijer-kernel density estimation on $\mathbb{R}^+$

abstract: Nonparametric kernel density estimation is a very natural procedure which simply makes use of the smoothing power of the convolution operation. Yet, it performs poorly when the density of a positive variable is to be estimated (boundary issues, spurious bumps in the tail). So various extensions of the basic kernel estimator allegedly suitable for $\mathbb{R}^+$-supported densities, such as those using Gamma or other asymmetric kernels, abound in the literature. Those, however, are not based on any valid smoothing operation analogous to the convolution, which typically leads to inconsistencies. By contrast, in this paper a kernel estimator for $\mathbb{R}^+$-supported densities is defined by making use of the Mellin convolution, the natural analogue of the usual convolution on $\mathbb{R}^+$. From there, a very transparent theory flows and leads to new type of asymmetric kernels strongly related to Meijer's $G$-functions. The numerous pleasant properties of this `Mellin-Meijer-kernel density estimator' are demonstrated in the paper. Its pointwise and $L_2$-consistency (with optimal rate of convergence) is established for a large class of densities, including densities unbounded at 0 and showing power-law decay in their right tail. Its practical behaviour is investigated further through simulations and some real data analyses.


title: The role of grain boundary scattering in reducing the thermal conductivity of polycrystalline XNiSn (X = Hf, Zr, Ti) half-Heusler alloys

abstract: Thermoelectric application of half-Heusler compounds suffers from their fairly high thermal conductivities. Insight into how effective various scattering mechanisms are in reducing the thermal conductivity of fabricated XNiSn compounds (X = Hf, Zr, Ti, and mixtures thereof) is therefore crucial. Here, we show that such insight can be obtained through a concerted theory-experiment comparison of how the lattice thermal conductivity kLat(T) depends on temperature and crystallite size. Comparing theory and experiment for a range of Hf0.5Zr0.5NiSn and ZrNiSn samples reported in the literature and in the present paper revealed that grain boundary scattering plays the most important role in bringing down kLat, in particular so for unmixed compounds. Our concerted analysis approach was corroborated by a good qualitative agreement between the measured and calculated kLat of polycrystalline samples, where the experimental average crystallite size was used as an input parameter for the calculations. The calculations were based on the Boltzmann transport equation and ab initio density functional theory. Our analysis explains the significant variation of reported kLat of nominally identical XNiSn samples and is expected to provide valuable insights into the dominant scattering mechanisms even for other materials.


title: A numerical study of Gibbs $u$-measures for partially hyperbolic diffeomorphisms on $\mathbb T^3$

abstract: We consider a hyperbolic automorphism $A\colon\mathbb T^3\to\mathbb T^3$ of the 3-torus whose 2-dimensional unstable distribution splits into weak and strong unstable subbundles. We unfold $A$ into two one-parameter families of Anosov diffeomorphisms --- a conservative family and a dissipative one. For diffeomorphisms in these families we numerically calculate the strong unstable manifold of the fixed point. Our calculations strongly suggest that the strong unstable manifold is dense in $\mathbb T^3$. Further, we calculate push-forwards of the Lebesgue measure on a local strong unstable manifold. These numeric data indicate that the sequence of push-forwards converges to the SRB measure.


title: Multiscale Surrogate Modeling and Uncertainty Quantification for Periodic Composite Structures

abstract: Computational modeling of the structural behavior of continuous fiber composite materials often takes into account the periodicity of the underlying micro-structure. A well established method dealing with the structural behavior of periodic micro-structures is the so- called Asymptotic Expansion Homogenization (AEH). By considering a periodic perturbation of the material displacement, scale bridging functions, also referred to as elastic correctors, can be derived in order to connect the strains at the level of the macro-structure with micro- structural strains. For complicated inhomogeneous micro-structures, the derivation of such functions is usually performed by the numerical solution of a PDE problem - typically with the Finite Element Method. Moreover, when dealing with uncertain micro-structural geometry and material parameters, there is considerable uncertainty introduced in the actual stresses experienced by the materials. Due to the high computational cost of computing the elastic correctors, the choice of a pure Monte-Carlo approach for dealing with the inevitable material and geometric uncertainties is clearly computationally intractable. This problem is even more pronounced when the effect of damage in the micro-scale is considered, where re-evaluation of the micro-structural representative volume element is necessary for every occurring damage. The novelty in this paper is that a non-intrusive surrogate modeling approach is employed with the purpose of directly bridging the macro-scale behavior of the structure with the material behavior in the micro-scale, therefore reducing the number of costly evaluations of corrector functions, allowing for future developments on the incorporation of fatigue or static damage in the analysis of composite structural components.


title: Typically bounding torsion

abstract: We formulate the notion of \emph{typical boundedness} of torsion on a family of abelian varieties defined over number fields. This means that the torsion subgroups of elements in the family can be made uniformly bounded by removing from the family all abelian varieties defined over number fields of degree lying in a set of arbitrarily small density. We show that for each fixed $g$, torsion is typically bounded on the family of all $g$-dimensional CM abelian varieties. We show that torsion is \emph{not} typically bounded on the family of all elliptic curves, and we establish results -- some unconditional and some conditional -- on typical boundedness of torsion of elliptic curves for which the degree of the $j$-invariant is fixed.


title: A Scalable Algorithm for Gaussian Graphical Models with Change-Points

abstract: Graphical models with change-points are computationally challenging to fit, particularly in cases where the number of observation points and the number of nodes in the graph are large. Focusing on Gaussian graphical models, we introduce an approximate majorize-minimize (MM) algorithm that can be useful for computing change-points in large graphical models. The proposed algorithm is an order of magnitude faster than a brute force search. Under some regularity conditions on the data generating process, we show that with high probability, the algorithm converges to a value that is within statistical error of the true change-point. A fast implementation of the algorithm using Markov Chain Monte Carlo is also introduced. The performances of the proposed algorithms are evaluated on synthetic data sets and the algorithm is also used to analyze structural changes in the S&P 500 over the period 2000-2016.


title: Clock spectroscopy of interacting bosons in deep optical lattices

abstract: We report on high-resolution optical spectroscopy of interacting bosonic $^{174}$Yb atoms in deep optical lattices with negligible tunneling. We prepare Mott insulator phases with singly- and doubly-occupied isolated sites and probe the atoms using an ultra-narrow "clock" transition. Atoms in singly-occupied sites undergo long-lived Rabi oscillations. Atoms in doubly-occupied sites are strongly affected by interatomic interactions, and we measure their inelastic decay rates and energy shifts. We deduce from these measurements all relevant collisional parameters involving both clock states, in particular the intra- and inter-state scattering lengths.


title: Optimizing EDELWEISS detectors for low-mass WIMP searches

abstract: The physics potential of EDELWEISS detectors for the search of low-mass Weakly Interacting Massive Particles (WIMPs) is studied. Using a data-driven background model, projected exclusion limits are computed using frequentist and multivariate analysis approaches, namely profile likelihood and boosted decision tree. Both current and achievable experimental performance are considered. The optimal strategy for detector optimization depends critically on whether the emphasis is put on WIMP masses below or above $\sim$ 5 GeV/c$^2$. The projected sensitivity for the next phase of the EDELWEISS-III experiment at the Modane Underground Laboratory (LSM) for low-mass WIMP search is presented. By 2018 an upper limit on the spin-independent WIMP-nucleon cross-section of $\sigma_{SI} = 7 \times 10^{-42}$ cm$^2$ is expected for a WIMP mass in the range 2$-$5 GeV/c$^2$. The requirements for a future hundred-kilogram scale experiment designed to reach the bounds imposed by the coherent scattering of solar neutrinos are also described. By improving the ionization resolution down to 50 eV$_{ee}$, we show that such an experiment installed in an even lower background environment (e.g. at SNOLAB) should allow to observe about 80 $^8$B neutrino events after discrimination.


title: A Dolbeault-Grothendieck Resolution for Singular Spaces

abstract: We construct a generalization of the Dolbeault-Grothendieck resolution on a singular complex space. The same construction yields, for each morphism of analytic spaces, a pullback mapping between the respective Dolbeault-Grothendieck resolutions. As in the smooth case, the terms of the resolution are soft sheaves with stalks which are flat with respect to the sheaf of holomorphic sections. If, moreover, the complex space $(X,\mathcal{O}_{X})$ is countable at infinity then the global section spaces of the terms of the resolution are endowed with natural Fr\'{e}chet-Schwarz topologies which induce the natural topology on the cohomology groups $H^{\bullet }(X,\mathcal{O}_{X})$. The construction is an exercise in globalization using semi-simplicial techniques. Using the above construction one can produce, for instance, a soft resolution with $\mathcal{O}_{X}$-flat stalks for the de Rham complex on the analytic space $X$.


title: Topological Sorting under Regular Constraints

abstract: We introduce the constrained topological sorting problem (CTS): given a regular language K and a directed acyclic graph G with labeled vertices, determine if G has a topological sort that forms a word in K. This natural problem applies to several settings, e.g., scheduling with costs or verifying concurrent programs. We consider the problem CTS[K] where the target language K is fixed, and study its complexity depending on K. We show that CTS[K] is tractable when K falls in several language families, e.g., unions of monomials, which can be used for pattern matching. However, we show that CTS[K] is NP-hard for K = (ab)^* and introduce a shuffle reduction technique to show hardness for more languages. We also study the special case of the constrained shuffle problem (CSh), where the input graph is a disjoint union of strings, and show that CSh[K] is additionally tractable when K is a group language or a union of district group monomials. We conjecture that a dichotomy should hold on the complexity of CTS[K] or CSh[K] depending on K, and substantiate this by proving a coarser dichotomy under a different problem phrasing which ensures that tractable languages are closed under common operators.


abstract: We prove existence of equilibrium states with special properties for a class of distance expanding local homeomorphisms on compact metric spaces and continuous potentials. Moreover, we formulate a C$^1$ generalization of Pesin's Entropy Formula: all ergodic weak-SRB-like measures satisfy Pesin's Entropy Formula for $C^1$ non-uniformly expanding maps. We show that for weak-expanding maps such that $\Leb$-a.e $x$ has positive frequency of hyperbolic times, then all the necessarily existing ergodic weak-SRB-like measures satisfy Pesin's Entropy Formula and are equilibrium states for the potential $\psi=-\log|\det Df|$. In particular, this holds for any $C^1$-expanding map and in this case the set of invariant probability measures that satisfy Pesin's Entropy Formula is the weak$^*$-closed convex hull of the ergodic weak-SRB-like measures.


title: Lempel-Ziv: a "one-bit catastrophe" but not a tragedy

abstract: The so-called "one-bit catastrophe" for the compression algorithm LZ'78 asks whether the compression ratio of an infinite word can change when a single bit is added in front of it. We answer positively this open question raised by Lutz and others: we show that there exists an infinite word $w$ such that $\rho_{sup}(w)=0$ but $\rho_{inf}(0w)>0$, where $\rho_{sup}$ and $\rho_{inf}$ are respectively the $\limsup$ and the $\liminf$ of the compression ratios $\rho$ of the prefixes. To that purpose we explore the behaviour of LZ'78 on finite words and show the following results: - There is a constant $C>0$ such that, for any finite word $w$ and any letter $a$, $\rho(aw)\leq C\sqrt{\rho(w)\log|w|}$. Thus, sufficiently compressible words ($\rho(w)=o(1/\log|w|)$) remain compressible with a letter in front; - The previous result is tight up to a multiplicative constant for any compression ratio $\rho(w)=O(1/\log|w|)$. In particular, there are infinitely many words $w$ satisfying $\rho(w)=O(1/\log|w|)$ but $\rho(0w)=\Omega(1)$.


title: Density matrix renormalization group study of a three-orbital Hubbard model with spin-orbit coupling in one dimension

abstract: Using the Density Matrix Renormalization Group technique we study the effect of spin-orbit coupling on a three-orbital Hubbard model in the $(t_{2g})^{4}$ sector and in one dimension. Fixing the Hund coupling to a robust value compatible with some multiorbital materials, we present the phase diagram varying the Hubbard $U$ and spin-orbit coupling $\lambda$, at zero temperature. Our results are shown to be qualitatively similar to those recently reported using the Dynamical Mean Field Theory in higher dimensions, providing a robust basis to approximate many-body techniques. Among many results, we observe an interesting transition from an orbital-selective Mott phase to an excitonic insulator with increasing $\lambda$ at intermediate $U$. In the strong $U$ coupling limit, we find a non-magnetic insulator with an effective angular momentum $\langle(\textbf{J}^{eff})^{2}\rangle \ne 0$ near the excitonic phase, smoothly connected to the $\langle(\textbf{J}^{eff})^{2}\rangle = 0$ regime. We also provide a list of quasi-one dimensional materials where the physics discussed in this publication could be realized.


title: Bayesian Optimization for Probabilistic Programs

abstract: We present the first general purpose framework for marginal maximum a posteriori estimation of probabilistic program variables. By using a series of code transformations, the evidence of any probabilistic program, and therefore of any graphical model, can be optimized with respect to an arbitrary subset of its sampled variables. To carry out this optimization, we develop the first Bayesian optimization package to directly exploit the source code of its target, leading to innovations in problem-independent hyperpriors, unbounded optimization, and implicit constraint satisfaction; delivering significant performance improvements over prominent existing packages. We present applications of our method to a number of tasks including engineering design and parameter optimization.


title: High-energy resummation in semi-hard processes at the LHC

abstract: Semi-hard processes in the large COM energy limit offer us an exclusive chance to test the dynamics behind strong interactions in kinematical sectors so far unexplored. In the Regge limit, fixed-order calculations in pQCD based on collinear factorisation miss the effect of large energy logs, which are so large to compensate the small QCD coupling $\alpha_s$ and must therefore be accounted for to all perturbative orders.. The BFKL approach represents the most powerful tool to resum to all orders these large logs both in the LLA and in the NLA. The inclusive hadroproduction of forward jets with high transverse momenta separated by a large rapidity interval at the LHC (Mueller--Navelet jets) has been one of the most studied reactions so far. Interesting observables associated to this process are the azimuthal correlations, showing a very good agreement with experimental data at the LHC. However, new BFKL-sensitive observables should be considered in the context of the LHC physics program. With the aim the to further and deeply probe the dynamics of QCD in the Regge limit, we give phenomenological predictions for four distinct semi-hard process. On one hand, we continue the analysis of reactions with two objects identified in the final state $(i)$ by addressing open problems in the Mueller--Navelet sector and $(ii)$ by studying the inclusive dihadron production in the full NLA BKFL accuracy. On the other hand, we show how inclusive multi-jet production processes allow us to define new, generalised and suitable BFKL observables, where transverse momenta and rapidities of the tagged jets, well separated in rapidity from each other, appear in new combinations. We give the first predictions for the inclusive three-jet production, encoding the effects of higher-order BFKL corrections. Then, making use of the same formalism, we present the first complete BFKL analysis for the four-jet production.


title: How hard is it to satisfy (almost) all roommates?

abstract: The classic Stable Roommates problem (which is the non-bipartite generalization of the well-known Stable Marriage problem) asks whether there is a stable matching for a given set of agents, i.e. a partitioning of the agents into disjoint pairs such that no two agents induce a blocking pair. Herein, each agent has a preference list denoting who it prefers to have as a partner, and two agents are blocking if they prefer to be with each other rather than with their assigned partners. Since stable matchings may not be unique, we study an NP-hard optimization variant of Stable Roommates, called Egal Stable Roommates, which seeks to find a stable matching with a minimum egalitarian cost {\gamma}, i.e. the sum of the dissatisfaction of the agents is minimum. The dissatisfaction of an agent is the number of agents that this agent prefers over its partner if it is matched; otherwise it is the length of its preference list. We also study almost stable matchings, called Min-Block-Pair Stable Roommates, which seeks to find a matching with a minimum number {\beta} of blocking pairs. Our main result is that Egal Stable Roommates parameterized by {\gamma} is fixed-parameter tractable, while Min-Block-Pair Stable Roommates parameterized by {\beta} is W[1]-hard, even if the length of each preference list is at most five.


title: Infinite rate symbiotic branching on the real line: The tired frogs model

abstract: Consider a population of infinitesimally small frogs on the real line. Initially the frogs on the positive half-line are dormant while those on the negative half-line are awake and move according to the heat flow. At the interface, the incoming wake frogs try to wake up the dormant frogs and succeed with a probability proportional to their amount among the total amount of involved frogs at the specific site. Otherwise, the incoming frogs also fall asleep. This frog model is a special case of the infinite rate symbiotic branching process on the real line with different motion speeds for the two types. We construct this frog model as the limit of approximating processes and compute the structure of jumps. We show that our frog model can be described by a stochastic partial differential equation on the real line with a jump type noise.


title: Discriminative Optimization: Theory and Applications to Computer Vision Problems

abstract: Many computer vision problems are formulated as the optimization of a cost function. This approach faces two main challenges: (i) designing a cost function with a local optimum at an acceptable solution, and (ii) developing an efficient numerical method to search for one (or multiple) of these local optima. While designing such functions is feasible in the noiseless case, the stability and location of local optima are mostly unknown under noise, occlusion, or missing data. In practice, this can result in undesirable local optima or not having a local optimum in the expected place. On the other hand, numerical optimization algorithms in high-dimensional spaces are typically local and often rely on expensive first or second order information to guide the search. To overcome these limitations, this paper proposes Discriminative Optimization (DO), a method that learns search directions from data without the need of a cost function. Specifically, DO explicitly learns a sequence of updates in the search space that leads to stationary points that correspond to desired solutions. We provide a formal analysis of DO and illustrate its benefits in the problem of 3D point cloud registration, camera pose estimation, and image denoising. We show that DO performed comparably or outperformed state-of-the-art algorithms in terms of accuracy, robustness to perturbations, and computational efficiency.


title: Model compression as constrained optimization, with application to neural nets. Part II: quantization

abstract: We consider the problem of deep neural net compression by quantization: given a large, reference net, we want to quantize its real-valued weights using a codebook with $K$ entries so that the training loss of the quantized net is minimal. The codebook can be optimally learned jointly with the net, or fixed, as for binarization or ternarization approaches. Previous work has quantized the weights of the reference net, or incorporated rounding operations in the backpropagation algorithm, but this has no guarantee of converging to a loss-optimal, quantized net. We describe a new approach based on the recently proposed framework of model compression as constrained optimization \citep{Carreir17a}. This results in a simple iterative "learning-compression" algorithm, which alternates a step that learns a net of continuous weights with a step that quantizes (or binarizes/ternarizes) the weights, and is guaranteed to converge to local optimum of the loss for quantized nets. We develop algorithms for an adaptive codebook or a (partially) fixed codebook. The latter includes binarization, ternarization, powers-of-two and other important particular cases. We show experimentally that we can achieve much higher compression rates than previous quantization work (even using just 1 bit per weight) with negligible loss degradation.


title: Residual Phase Noise Measurement of Optical Second Harmonic Generation in PPLN Waveguides

abstract: We report on the characterization, including residual phase noise and fractional frequency instability, of fiber-coupled PPLN non-linear crystals. These components are devoted to frequency doubling 871 nm light from an extended-cavity diode laser to produce a 435.5 nm beam, corresponding to the ytterbium ion electric quadrupole clock transition. We measure doubling efficiencies of up to 117.5 %/W. Using a Mach-Zehnder interferometer and an original noise rejection technique, the residual phase noise of the doublers is estimated to be lower than ${\rm -35\, dBrad^2/Hz}$ at 1 Hz, making these modules compatible with up-to-date optical clocks and ultra-stable cavities. The influence of external parameters such as pump laser frequency and intensity is investigated, showing that they do not limit the stability of the frequency-doubled signal. Our results demonstrate that such compact, fiber-coupled modules are suitable for use in ultra-low phase noise metrological experiments, including transportable optical atomic clocks.


title: Feasibility and applications of the spin-echo modulation option for a small angle neutron scattering instrument at the European Spallation Source

abstract: We describe theoretical and practical aspects of spin-echo modulated small-angle neutron scattering (SEMSANS) as well as the potential combination with SANS. Based on the preliminary technical designs of SKADI (a SANS instrument proposed for the European Spallation Source) and a SEMSANS add-on, we assess the practicability, feasibility and scientific merit of a combined SANS and SEMSANS setup by calculating tentative SANS and SEMSANS results for soft matter, geology and advanced material samples that have been previously studied by scattering methods. We conclude that lengths from 1 nm up to 0.01 mm can be observed simultaneously in a single measurement. Thus, the combination of SANS and SEMSANS instrument is suited for the simultaneous observation of a wide range of length scales, e.g. for time-resolved studies of kinetic processes in complex multiscale systems.


title: Microresonator-based high-resolution gas spectroscopy

abstract: In recent years, microresonator-based optical frequency combs have created up opportunities for developing a spectroscopy laboratory on a chip due to its broadband emission and high comb power. However, with mode spacings typically in the range of 10 - 1000 GHz, the realization of a chip-based high-resolution spectrometer suitable for gas-phase spectroscopy has proven to be difficult. Here, we show mode-hop-free tuning of a microresonator-based frequency comb over 16 GHz by simultaneously tuning both the pump laser and the cavity resonance. We illustrate the power of this scanning technique by demonstrating gas-phase molecular fingerprinting of acetylene with a high-spectral-resolution of < 80 MHz over a 45-THz optical bandwidth in the mid-IR. Our technique represents a significant step towards on-chip gas sensing with an ultimate spectral resolution given by the comb linewidth.


title: Ab-Initio Pulsar Magnetosphere: Particle acceleration in Oblique Rotators and High-energy Emission Modeling

abstract: We perform global particle-in-cell simulations of pulsar magnetospheres including pair production, ion extraction from the surface, frame dragging corrections, and high energy photon emission and propagation. In the case of oblique rotators, effects of general relativity increase the fraction of open field lines which support active pair discharge. We find that the plasma density and particle energy flux in the pulsar wind are highly non-uniform with latitude. Significant fraction of the outgoing particle energy flux is carried by energetic ions, which are extracted from the stellar surface. Their energies may extend up to a large fraction of the open field line voltage, making them interesting candidates for ultra-high-energy cosmic rays. We show that pulsar gamma-ray radiation is dominated by synchrotron emission, produced by particles that are energized by relativistic magnetic reconnection close to the Y-point and in the equatorial current sheet. In most cases, calculated light curves contain two strong peaks, in general agreement with Fermi observations. The radiative efficiency decreases with increasing pulsar inclination and increasing efficiency of pair production in the current sheet, explaining the observed scatter in $L_{\gamma}$ vs $\dot{E}$. We find that the high-frequency cutoff in the spectra is regulated by the pair loading of the current sheet. Our findings lay the foundation for quantitative interpretation of Fermi observations of gamma-ray pulsars.


title: Tensor-Based Backpropagation in Neural Networks with Non-Sequential Input

abstract: Neural networks have been able to achieve groundbreaking accuracy at tasks conventionally considered only doable by humans. Using stochastic gradient descent, optimization in many dimensions is made possible, albeit at a relatively high computational cost. By splitting training data into batches, networks can be distributed and trained vastly more efficiently and with minimal accuracy loss. We have explored the mathematics behind efficiently implementing tensor-based batch backpropagation algorithms. A common approach to batch training is iterating over batch items individually. Explicitly using tensor operations to backpropagate allows training to be performed non-linearly, increasing computational efficiency.


title: Tracking capacitance of liquid crystal devices to improve polarization rotation accuracy

abstract: We report a capacitance tracking method for achieving arbitrary polarization rotation from nematic liquid crystals. By locking to the unique capacitance associated with the molecular orientation, any polarization rotation can be achieved with improved accuracy over a wide temperature range. A modified relaxation oscillator circuit that can simultaneously determine the capacitance and drive the rotator is presented.


title: Quantitative isoperimetry \`a la Levy-Gromov

abstract: On a Riemannian manifold with a positive lower bound on the Ricci tensor, the distance of isoperimetric sets from geodesic balls is quantitatively controlled in terms of the gap between the isoperimetric profile of the manifold and that of a round sphere of suitable radius. The deficit between the diameters of the manifold and of the corresponding sphere is bounded likewise. These results are actually obtained in the more general context of (possibly non-smooth) metric measure spaces with curvature-dimension conditions through a quantitative analysis of the transport-rays decompositions obtained by the localization method.


title: Human-Level Intelligence or Animal-Like Abilities?

abstract: The vision systems of the eagle and the snake outperform everything that we can make in the laboratory, but snakes and eagles cannot build an eyeglass or a telescope or a microscope. (Judea Pearl)


title: Generalized stealthy hyperuniform processes : maximal rigidity and the bounded holes conjecture

abstract: We study translation invariant stochastic processes on $\mathbb{R}^d$ or $\mathbb{Z}^d$ whose diffraction spectrum or structure function $S(k)$, i.e. the Fourier transform of the truncated total pair correlation function, vanishes on an open set $U$ in the wave space. A key family of such processes are stealthy hyperuniform point processes, for which the origin $k=0$ is in $U$; these are of much current physical interest. We show that all such processes exhibit the following remarkable maximal rigidity : namely, the configuration outside a bounded region determines, with probability 1, the exact value (or the exact locations of the points) of the process inside the region. In particular, such processes are completely determined by their tail. In the 1D discrete setting (i.e. $\mathbb{Z}$-valued processes on $\mathbb{Z}$), this can also be seen as a consequence of a recent theorem of Borichev, Sodin and Weiss; in higher dimensions or in the continuum, such a phenomenon seems novel. For stealthy hyperuniform point processes, we prove the Zhang-Stillinger-Torquato conjecture that such processes have bounded holes (empty regions), with a universal bound that depends inversely on the size of $U$.

categories: math.PR cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.stat-mech math-ph math.MP


title: Monitor variability of millimeter lines in IRC +10216

abstract: A single dish monitoring of millimeter maser lines SiS J=14-13 and HCN nu_2 = 1^f J=3-2 and several other rotational lines is reported for the archetypal carbon star IRC+10216. Relative line strength variations of 5%~30% are found for eight molecular line features with respect to selected reference lines. Definite line-shape variation is found in limited velocity intervals of the SiS and HCN line profiles. The asymmetrical line profiles of the two lines are mainly due to the varying components. Their dominant varying components of the line profiles have similar periods and phases as the IR light variation, although both quantities show some degree of velocity dependence; there is also variability asymmetry between the blue and red line wings of both lines. Combining the velocities and amplitudes with a wind velocity model, we suggest that the line profile variations are due to SiS and HCN masing lines emanating from the wind acceleration zone. The possible link of the variabilities to thermal, dynamical and/or chemical processes within or under this region is also discussed.


title: Open Chemistry: RESTful Web APIs, JSON, NWChem and the Modern Web Application

abstract: An end-to-end platform for chemical science research has been developed that integrates data from computational and experimental approaches through a modern web-based interface. The platform offers a highly interactive visualization and analytics environment that functions well on mobile, laptop and desktop devices. It offers pragmatic solutions to ensure that large and complex data sets are more accessible. Existing desktop applications/frameworks were extended to integrate with high-performance computing (HPC) resources, and offer command-line tools to automate interaction---connecting distributed teams to this software platform on their own terms. The platform was developed openly, and all source code hosted on the GitHub platform with automated deployment possible using Ansible coupled with standard Ubuntu-based machine images deployed to cloud machines. The platform is designed to enable teams to reap the benefits of the connected web---going beyond what conventional search and analytics platforms offer in this area. It also has the goal of offering federated instances, that can be customized to the sites/research performed. Data gets stored using JSON, extending upon previous approaches using XML, building structures that support computational chemistry calculations. These structures were developed to make it easy to process data across different languages, and send data to a JavaScript web client.


title: A Tight Approximation for Co-flow Scheduling for Minimizing Total Weighted Completion Time

abstract: Co-flows model a modern scheduling setting that is commonly found in a variety of applications in distributed and cloud computing. In co-flow scheduling, there are $m$ input ports and $m$ output ports. Each co-flow $j \in J$ can be represented by a bipartite graph between the input and output ports, where each edge $(i,o)$ with demand $d_{i,o}^j$ means that $d_{i,o}^j$ units of packets must be delivered from port $i$ to port $o$. To complete co-flow $j$, we must satisfy all of its demands. Due to capacity constraints, a port can only transmit (or receive) one unit of data in unit time. A feasible schedule at each time $t$ must therefore be a bipartite matching. We consider co-flow scheduling and seek to optimize the popular objective of total weighted completion time. Our main result is a $(2+\epsilon)$-approximation for this problem, which is essentially tight, as the problem is hard to approximate within a factor of $(2 - \epsilon)$. This improves upon the previous best known 4-approximation. Further, our result holds even when jobs have release times without any loss in the approximation guarantee. The key idea of our approach is to construct a continuous-time schedule using a configuration linear program and interpret each job's completion time therein as the job's deadline. The continuous-time schedule serves as a witness schedule meeting the discovered deadlines, which allows us to reduce the problem to a deadline-constrained scheduling problem. * This result is flawed; see the first page for the details.


title: Brittle to Quasi-Brittle Transition and Crack Initiation Precursors in Disordered Crystals

abstract: Crack initiation emerges due to a combination of elasticity, plasticity, and disorder, and it is heavily dependent on the material's microstructural details. In this paper, we investigate brittle metals with coarse-grained, microstructural disorder that could originate in a material's manufacturing process, such as alloying. As an investigational tool, we consider crack initiation from a surface, ellipsoidal notch: As the radius of curvature at the notch increases, there is a dynamic transition from notch-induced crack initiation to bulk-disorder crack nucleation. We perform extensive and realistic simulations using a phase-field approach coupled to crystal plasticity. Furthermore, the microstructural disorder and notch width are varied in order to study the transition. We identify this transition for various disorder strengths in terms of the damage evolution. Above the transition, we identify detectable precursors to crack initiation that we quantify in terms of the expected stress drops during mode I fracture loading. We discuss ways to observe and analyze this brittle to quasi-brittle transition in experiments.

categories: cond-mat.mtrl-sci cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech


title: Relativistic Extension of a Charge-Conservative Finite Element Solver for Time-Dependent Maxwell-Vlasov Equations

abstract: In many problems involving particle accelerators and relativistic plasmas, the accurate modeling of relativistic particle motion is essential for accurate physical predictions. Here, we extend a charge-conserving finite element time-domain (FETD) particle-in-cell (PIC) algorithm for the time-dependent Maxwell-Vlasov equations on irregular (unstructured) meshes to the relativistic regime by implementing and comparing three particle pushers: (relativistic) Boris, Vay, and Higuera-Cary. We illustrate the application of the proposed relativistic FETD-PIC algorithm for the analysis of particle cyclotron motion at relativistic speeds, harmonic particle oscillation in the Lorentz-boosted frame, and relativistic Bernstein modes in magnetized charge-neutral (pair) plasmas.


title: Implications of a wavelength dependent PSF for weak lensing measurements

abstract: The convolution of galaxy images by the point-spread function (PSF) is the dominant source of bias for weak gravitational lensing studies, and an accurate estimate of the PSF is required to obtain unbiased shape measurements. The PSF estimate for a galaxy depends on its spectral energy distribution (SED), because the instrumental PSF is generally a function of the wavelength. In this paper we explore various approaches to determine the resulting `effective' PSF using broad-band data. Considering the Euclid mission as a reference, we find that standard SED template fitting methods result in biases that depend on source redshift, although this may be remedied if the algorithms can be optimised for this purpose. Using a machine-learning algorithm we show that, at least in principle, the required accuracy can be achieved with the current survey parameters. It is also possible to account for the correlations between photometric redshift and PSF estimates that arise from the use of the same photometry. We explore the impact of errors in photometric calibration, errors in the assumed wavelength dependence of the PSF model and limitations of the adopted template libraries. Our results indicate that the required accuracy for Euclid can be achieved using the data that are planned to determine photometric redshifts.


title: Reconciling solar and stellar magnetic cycles with nonlinear dynamo simulations

abstract: The magnetic fields of solar-type stars are observed to cycle over decadal periods -11 years in the case of the Sun. The fields originate in the turbulent convective layers of stars and have a complex dependency upon stellar rotation rate. We have performed a set of turbulent global simulations that exhibit magnetic cycles varying systematically with stellar rotation and luminosity. We find that the magnetic cycle period is inversely proportional to the Rossby number, which quantifies the influence of rotation on turbulent convection. The trend relies on a fundamentally non-linear dynamo process and is compatible with the Sun's cycle and those of other solar-type stars.


title: Charge order in SrCoO$_{3-y}$ investigated through resonant soft x-ray scattering

abstract: Resonant soft X-ray scattering was used to determine the presence of more subtle orderings not detected in standard structural analyses. By tuning to specific Co absorption edges, arrangements particular to the electronic states of those elements are enhanced. We have discovered an ordering commensurate to the lattice at the ($\frac{1}{4}~\frac{1}{4}~\frac{1}{4}$) position. Incommensurate peaks near this position were also observed. The intensity of these peaks depends on the oxygen concentration of the sample, and can be suppressed at temperatures above 320 K. Regular orderings of charge density which closely match the underlying lattice may help to explain the observed propensity for SrCoO$_{3-y}$ (0 $\leq$ y $\leq\frac{1}{2}$) to stabilize at particular phases.


title: Films of Bacteria at Interfaces (FBI): Remodeling of Fluid Interfaces by Pseudomonas aeruginosa

abstract: Bacteria at fluid interfaces endure physical and chemical stresses unique to these highly asymmetric environments. The responses of Pseudomonas aeruginosa PAO1 and PA14 to a hexadecane-water interface are compared. PAO1 cells form elastic films of bacteria, excreted polysaccharides and proteins, whereas PA14 cells move actively without forming an elastic film. Studies of PAO1 mutants show that, unlike solid-supported biofilms, elastic interfacial film formation occurs in the absence of flagella, pili, or certain polysaccharides. Highly induced genes identified in transcriptional profiling include those for putative enzymes and a carbohydrate metabolism enzyme, alkB2; this latter gene is not upregulated in PA14 cells. Notably, PAO1 mutants lacking the alkB2 gene fail to form an elastic layer. Rather, they form an active film like that formed by PA14. These findings demonstrate that genetic expression is altered by interfacial confinement, and suggest that the ability to metabolize alkanes may play a role in elastic film formation at oil-water interfaces.


title: ADMM Based Privacy-preserving Decentralized Optimization

abstract: Privacy preservation is addressed for decentralized optimization, where $N$ agents cooperatively minimize the sum of $N$ convex functions private to these individual agents. In most existing decentralized optimization approaches, participating agents exchange and disclose states explicitly, which may not be desirable when the states contain sensitive information of individual agents. The problem is more acute when adversaries exist which try to steal information from other participating agents. To address this issue, we propose a privacy-preserving decentralized optimization approach based on ADMM and partially homomorphic cryptography. To our knowledge, this is the first time that cryptographic techniques are incorporated in a fully decentralized setting to enable privacy preservation in decentralized optimization in the absence of any third party or aggregator. To facilitate the incorporation of encryption in a fully decentralized manner, we introduce a new ADMM which allows time-varying penalty matrices and rigorously prove that it has a convergence rate of $O(1/t)$. Numerical and experimental results confirm the effectiveness and low computational complexity of the proposed approach.


title: Some cases of a conjecture on L-functions of twisted Carlitz modules

abstract: We prove two polynomial identities which are particular cases of a conjecture arising in the theory of L-functions of twisted Carlitz modules. This conjecture is stated in earlier papers of the second author.


title: Discord, quantum knowledge and private communications

abstract: In this brief review, we discuss the role that quantum correlations, as quantified by quantum discord, play in two interesting settings. The first one is discerning which unitaries have been applied on a quantum system, by taking advantage of knowledge regarding its initial configuration. Here discord captures the "quantum" component of this knowledge, useful only when we have access to a quantum memory. In particular, discord can be used to detect whether an untrusted party has certain quantum capabilities. The second setting is quantum cryptography. Here discord represents an important resource for trusted-noise quantum key distribution and also provides a general upper bound for the optimal secret key rates that are achievable by ideal protocols. In particular, the (two-way assisted) secret key capacity of a lossy bosonic channel exactly coincides with the maximum discord that can be distributed between the remote parties at the two ends of the channel.


title: Associative idempotent nondecreasing functions are reducible

abstract: An $n$-ary associative function is called reducible if it can be written as a composition of a binary associative function. We summarize known results when the function is defined on a chain and is nondecreasing. Our main result shows that associative idempotent and nondecreasing functions are uniquely reducible.


title: Constructions of cyclic constant dimension codes

abstract: Subspace codes and particularly constant dimension codes have attracted much attention in recent years due to their applications in random network coding. As a particular subclass of subspace codes, cyclic subspace codes have additional properties that can be applied efficiently in encoding and decoding algorithms. It is desirable to find cyclic constant dimension codes such that both the code sizes and the minimum distances are as large as possible. In this paper, we explore the ideas of constructing cyclic constant dimension codes proposed in \big([2], IEEE Trans. Inf. Theory, 2016\big) and \big([17], Des. Codes Cryptogr., 2016\big) to obtain further results. Consequently, new code constructions are provided and several previously known results in [2] and [17] are extended.


title: Stable processes, self-similarity and the unit ball

abstract: Around the 1960s a celebrated collection of papers emerged offering a number of explicit identities for the class of isotropic stable processes in one and higher dimensions; these include, for example, the lauded works of Blumenthal, Getoor, Ray, Port and Rogozin. Amongst other things, these results nicely exemplify the use of standard Riesz potential theory on the unit open ball $\mathbb{B}_d :=\{ x\in \mathbb{R}^d :|x|<1\}$, $\mathbb{R}^d\backslash\mathbb{B}_d$ and $\mathbb{S}_d :=\partial \mathbb{B}_d$ with the, then, modern theory of potential analysis for Markov processes. Following initial observations of Lamperti in 1972, with the occasional sporadic work of Kiu, Vuolle-Apiala and Graversen in the 1980s, an alternative understanding of stable processes through the theory of self-similar Markov processes has prevailed in the last decade or more. This point of view offers deeper probabilistic insights into some of the aforementioned potential analytical relations. In this review article, we will rediscover many of the aforementioned classical identities in relation to the unit ball and recent extensions thereof by combining elements of these two theories, which have otherwise been largely separated by decades in the literature. We present a dialogue that appeals as much as possible to path decompositions. Most notable in this respect is the Lamperti-Kiu decomposition of self-similar Markov processes and the Riesz-Bogdan-\.Zak transformation. Some of the proofs we give are known, some are known proofs mixed with new methods and some proofs are completely new. We assume that the reader has a degree of familiarity with the bare basics of L\'evy processes but nonetheless, we often include reminders of standard material.


title: Probing many-body dynamics on a 51-atom quantum simulator

abstract: Controllable, coherent many-body systems can provide insights into the fundamental properties of quantum matter, enable the realization of new quantum phases and could ultimately lead to computational systems that outperform existing computers based on classical approaches. Here we demonstrate a method for creating controlled many-body quantum matter that combines deterministically prepared, reconfigurable arrays of individually trapped cold atoms with strong, coherent interactions enabled by excitation to Rydberg states. We realize a programmable Ising-type quantum spin model with tunable interactions and system sizes of up to 51 qubits. Within this model, we observe phase transitions into spatially ordered states that break various discrete symmetries, verify the high-fidelity preparation of these states and investigate the dynamics across the phase transition in large arrays of atoms. In particular, we observe robust manybody dynamics corresponding to persistent oscillations of the order after a rapid quantum quench that results from a sudden transition across the phase boundary. Our method provides a way of exploring many-body phenomena on a programmable quantum simulator and could enable realizations of new quantum algorithms.


title: Gaussian Graphical Models: An Algebraic and Geometric Perspective

abstract: Gaussian graphical models are used throughout the natural sciences, social sciences, and economics to model the statistical relationships between variables of interest in the form of a graph. We here provide a pedagogic introduction to Gaussian graphical models and review recent results on maximum likelihood estimation for such models. Throughout, we highlight the rich algebraic and geometric properties of Gaussian graphical models and explain how these properties relate to convex optimization and ultimately result in insights on the existence of the maximum likelihood estimator (MLE) and algorithms for computing the MLE.


title: Pushing Configuration-Interaction to the Limit: Towards Massively Parallel MCSCF Calculations

abstract: A new large-scale parallel multiconfigurational self-consistent field (MCSCF) implementation in the open-source NWChem computational chemistry code is presented. The generalized active space (GAS) approach is used to partition large configuration interaction (CI) vectors and generate a sufficient number of batches that can be distributed to the available nodes. Massively parallel CI calculations with large active spaces can be treated. The performance of the new parallel MCSCF implementation is presented for the chromium trimer and for an active space of 20 electrons in 20 orbitals. Unprecedented CI calculations with an active space of 22 electrons in 22 orbitals for the pentacene systems were performed and a single CI iteration calculation with an active space of 24 electrons in 24 orbitals for the chromium tetramer was possible. The chromium tetramer corresponds to a CI expansion of one trillion SDs (914 058 513 424) and is largest conventional CI calculation attempted up to date.


title: Weakly Submodular Maximization Beyond Cardinality Constraints: Does Randomization Help Greedy?

abstract: Submodular functions are a broad class of set functions, which naturally arise in diverse areas. Many algorithms have been suggested for the maximization of these functions. Unfortunately, once the function deviates from submodularity, the known algorithms may perform arbitrarily poorly. Amending this issue, by obtaining approximation results for set functions generalizing submodular functions, has been the focus of recent works. One such class, known as weakly submodular functions, has received a lot of attention. A key result proved by Das and Kempe (2011) showed that the approximation ratio of the greedy algorithm for weakly submodular maximization subject to a cardinality constraint degrades smoothly with the distance from submodularity. However, no results have been obtained for maximization subject to constraints beyond cardinality. In particular, it is not known whether the greedy algorithm achieves any non-trivial approximation ratio for such constraints. In this paper, we prove that a randomized version of the greedy algorithm (previously used by Buchbinder et al. (2014) for a different problem) achieves an approximation ratio of $(1 + 1/\gamma)^{-2}$ for the maximization of a weakly submodular function subject to a general matroid constraint, where $\gamma$ is a parameter measuring the distance of the function from submodularity. Moreover, we also experimentally compare the performance of this version of the greedy algorithm on real world problems against natural benchmarks, and show that the algorithm we study performs well also in practice. To the best of our knowledge, this is the first algorithm with a non-trivial approximation guarantee for maximizing a weakly submodular function subject to a constraint other than the simple cardinality constraint. In particular, it is the first algorithm with such a guarantee for the important and broad class of matroid constraints.


title: Natural Boundary Conditions for Smoothing in Geometry Processing

abstract: In geometry processing, smoothness energies are commonly used to model scattered data interpolation, dense data denoising, and regularization during shape optimization. The squared Laplacian energy is a popular choice of energy and has a corresponding standard implementation: squaring the discrete Laplacian matrix. For compact domains, when values along the boundary are not known in advance, this construction bakes in low-order boundary conditions. This causes the geometric shape of the boundary to strongly bias the solution. For many applications, this is undesirable. Instead, we propose using the squared Frobenious norm of the Hessian as a smoothness energy. Unlike the squared Laplacian energy, this energy's natural boundary conditions (those that best minimize the energy) correspond to meaningful high-order boundary conditions. These boundary conditions model free boundaries where the shape of the boundary should not bias the solution locally. Our analysis begins in the smooth setting and concludes with discretizations using finite-differences on 2D grids or mixed finite elements for triangle meshes. We demonstrate the core behavior of the squared Hessian as a smoothness energy for various tasks.


title: Categorical diagonalization

abstract: This paper lays the groundwork for the theory of categorical diagonalization. Given a diagonalizable operator, tools in linear algebra (such as Lagrange interpolation) allow one to construct a collection of idempotents which project to each eigenspace. These idempotents are mutually orthogonal and sum to the identity. We categorify these tools. At the categorical level, one has not only eigenobjects and eigenvalues but also eigenmaps, which relate an endofunctor to its eigenvalues. Given an invertible endofunctor of a triangulated category with a sufficiently nice collection of eigenmaps, we construct idempotent functors which project to eigencategories. These idempotent functors are mutually orthogonal, and a convolution thereof is isomorphic to the identity functor. In several sequels to this paper, we will use this technology to study the categorical representation theory of Hecke algebras. In particular, for Hecke algebras of type A, we will construct categorified Young symmetrizers by simultaneously diagonalizing certain functors associated to the full twist braids.


title: Indication for Double Parton Scatterings in W + Prompt J/Psi Production at the LHC

abstract: We re-analyse the associated production of a prompt J/psi and a W boson in pp collisions at the LHC following the results of the ATLAS Collaboration. We perform the first study of the Single-Parton-Scattering (SPS) contributions at the Next-to-Leading Order (NLO) in alpha_s in the Colour-Evaporation Model (CEM), an approach based on the quark-hadron-duality. Our study provides clear indications for Double-Parton-Scattering (DPS) contributions, in particular at low transverse momenta, since our SPS CEM evaluation, which can be viewed as a conservative upper limit of the SPS yields, falls short compared to the ATLAS experimental data by 3.1 standard deviations. We also determine a finite allowed region for sigma_eff, inversely proportional to the size of the DPS yields, corresponding to the otherwise opposed hypotheses, namely our NLO CEM evaluation and the LO direct Colour-Singlet (CS) Model contribution. In both cases, the resulting DPS yields are significantly larger than that initially assumed by ATLAS based on jet-related analyses but is consistent with their observed raw-yield azimuthal distribution and with their prompt J/psi+J/psi and Z + prompt J/psi data.


title: A Generating Function for the Distribution of Runs in Binary Words

abstract: Let $N(n,r,k)$ denote the number of binary words of length $n$ that begin with $0$ and contain exactly $k$ runs (i.e., maximal subwords of identical consecutive symbols) of length $r$. We show that the generating function for the sequence $N(n,r,0)$, $n=0,1,\ldots$, is $(1-x)(1-2x + x^r-x^{r+1})^{-1}$ and that the generating function for $\{N(n,r,k)\}$ is $x^{kr}$ time the $k+1$ power of this. We extend to counts of words containing exactly $k$ runs of $1$s by using symmetries on the set of binary words.


title: Advances in Artificial Intelligence Require Progress Across all of Computer Science

abstract: Advances in Artificial Intelligence require progress across all of computer science.


title: Effective field theory for triaxially deformed nuclei

abstract: Effective field theory (EFT) is generalized to investigate the rotational motion of triaxially deformed even-even nuclei. A Hamiltonian, called the triaxial rotor model (TRM), is obtained up to next-to-leading order (NLO) within the EFT formalism. Its applicability is examined by comparing with a five-dimensional collective Hamiltonian (5DCH) for the description of the energy spectra of the ground state and $\gamma$ band in Ru isotopes. It is found that by taking into account the NLO corrections, the ground state band in the whole spin region and the $\gamma$ band in the low spin region are well described. The results presented here indicate that it should be possible to further generalize the EFT to triaxial nuclei with odd mass number.


title: Teleportation Through the Wormhole

abstract: ER=EPR allows us to think of quantum teleportation as communication of quantum information through space-time wormholes connecting entangled systems. The conditions for teleportation render the wormhole traversable so that a quantum system entering one end of the ERB will, after a suitable time, appear at the other end. Teleportation requires the transfer of classical information outside the horizon, but the classical bit-string carries no information about the teleported system; the teleported system passes through the ERB leaving no trace outside the horizon. In general the teleported system will retain a memory of what it encountered in the wormhole. This phenomenon could be observable in a laboratory equipped with quantum computers.


title: On the arithmetic of simple singularities of type E

abstract: An ADE Dynkin diagram gives rise to a family of algebraic curves. In this paper, we use arithmetic invariant theory to study the integral points of the curves associated to the exceptional diagrams $E_6, E_7$, $E_8$. These curves are non-hyperelliptic of genus 3 or 4. We prove that a positive proportion of each family consists of curves with integral points everywhere locally but no integral points globally.


title: Gain stabilization in Micro Pattern Gaseous Detectors: methodology and results

abstract: The phenomenon of avalanche-gain variations over time, particularly in Micro Pattern Gaseous Detectors (MPGD) incorporating insulator materials, have been generally attributed to electric-field modifications resulting from "charging-up" effects of the insulator. A robust methodology for characterization of gain-transients in such detectors is presented. It comprises three guidelines: detector initialization, long-gain stabilization monitoring and imposing transients by applying abrupt changes in operation conditions. Using THWELL and RPWELL detectors, we validated the proposed methodology by assessing a charging-up/charging-down model describing the governing processes of gain stabilization. The results provide a deeper insight into these processes, reflected by different transients upon abrupt variations of detector gain or the irradiation rate. This methodology provides a handle for future investigations of the involved physics phenomena in MPGD detectors comprising insulating components.


title: Probing low noise at the MOS interface with a spin-orbit qubit

abstract: The silicon metal-oxide-semiconductor (MOS) material system is technologically important for the implementation of electron spin-based quantum information technologies. Researchers predict the need for an integrated platform in order to implement useful computation, and decades of advancements in silicon microelectronics fabrication lends itself to this challenge. However, fundamental concerns have been raised about the MOS interface (e.g. trap noise, variations in electron g-factor and practical implementation of multi-QDs). Furthermore, two-axis control of silicon qubits has, to date, required the integration of non-ideal components (e.g. microwave strip-lines, micro-magnets, triple quantum dots, or introduction of donor atoms). In this paper, we introduce a spin-orbit (SO) driven singlet-triplet (ST) qubit in silicon, demonstrating all-electrical two-axis control that requires no additional integrated elements and exhibits charge noise properties equivalent to other more model, but less commercially mature, semiconductor systems. We demonstrate the ability to tune an intrinsic spin-orbit interface effect, which is consistent with Rashba and Dresselhaus contributions that are remarkably strong for a low spin-orbit material such as silicon. The qubit maintains the advantages of using isotopically enriched silicon for producing a quiet magnetic environment, measuring spin dephasing times of 1.6 $\mu$s using 99.95% $^{28}$Si epitaxy for the qubit, comparable to results from other isotopically enhanced silicon ST qubit systems. This work, therefore, demonstrates that the interface inherently provides properties for two-axis control, and the technologically important MOS interface does not add additional detrimental qubit noise.


title: The quark propagator in the complex domain: A status report

abstract: In these proceedings I review the status of an ongoing project that aims at solving the quark propagator Dyson-Schwinger equation in the complex domain. The novel aspect of the approach is that the non-analyticities arising throughout the iteration of the equation are to be taken into account in an appropriate way through dynamic contour deformation. Because of the complexity of the approach, these studies are undertaken in a heavily truncated scenario that serves as a toy model for the development of the numerical techniques that are required to treat the system in a mathematically sound way.


title: Spin 1 low lying meson spectra and the subtle link to the spin 0 mesons

abstract: An NJL-type three-flavor quark model with a complete set of explicit chiral symmetry breaking terms is extended to include vector and axial vector effective interaction terms. The bosonized Lagrangian is written up to quadratic order in the bosonic fields and the role of the new interactions is analysed in detail. The model's parameters are fitted to yield reasonable values to the four low-lying spin 0 and spin 1 meson nonets' masses.


title: Derivative Principal Component Analysis for Representing the Time Dynamics of Longitudinal and Functional Data

abstract: We propose a nonparametric method to explicitly model and represent the derivatives of smooth underlying trajectories for longitudinal data. This representation is based on a direct Karhunen--Lo\`eve expansion of the unobserved derivatives and leads to the notion of derivative principal component analysis, which complements functional principal component analysis, one of the most popular tools of functional data analysis. The proposed derivative principal component scores can be obtained for irregularly spaced and sparsely observed longitudinal data, as typically encountered in biomedical studies, as well as for functional data which are densely measured. Novel consistency results and asymptotic convergence rates for the proposed estimates of the derivative principal component scores and other components of the model are derived under a unified scheme for sparse or dense observations and mild conditions. We compare the proposed representations for derivatives with alternative approaches in simulation settings and also in a wallaby growth curve application. It emerges that representations using the proposed derivative principal component analysis recover the underlying derivatives more accurately compared to principal component analysis-based approaches especially in settings where the functional data are represented with only a very small number of components or are densely sampled. In a second wheat spectra classification example, derivative principal component scores were found to be more predictive for the protein content of wheat than the conventional functional principal component scores.


title: Integrability of scalar curvature and normal metric on conformally flat manifolds

abstract: On a manifold $(\mathbb{R}^n, e^{2u} |dx|^2)$, we say $u$ is normal if the $Q$-curvature equation that $u$ satisfies $(-\Delta)^{\frac{n}{2}} u = Q_g e^{nu}$ can be written as the integral form $u(x)=\frac{1}{c_n}\int_{\mathbb R^n}\log\frac{|y|}{|x-y|}Q_g(y)e^{nu(y)}dy+C$. In this paper, we show that the integrability assumption on the negative part of the scalar curvature implies the metric is normal. As an application, we prove a bi-Lipschitz equivalence theorem for conformally flat metrics.


title: Hyppo-X: A Scalable Exploratory Framework for Analyzing Complex Phenomics Data

abstract: Phenomics is an emerging branch of modern biology that uses high throughput phenotyping tools to capture multiple environmental and phenotypic traits, often at massive spatial and temporal scales. The resulting high dimensional data represent a treasure trove of information for providing an in-depth understanding of how multiple factors interact and contribute to the overall growth and behavior of different genotypes. However, computational tools that can parse through such complex data and aid in extracting plausible hypotheses are currently lacking. In this paper, we present Hyppo-X, a new algorithmic approach to visually explore complex phenomics data and in the process characterize the role of environment on phenotypic traits. We model the problem as one of unsupervised structure discovery, and use emerging principles from algebraic topology and graph theory for discovering higher-order structures of complex phenomics data. We present an open source software which has interactive visualization capabilities to facilitate data navigation and hypothesis formulation. We test and evaluate Hyppo-X on two real-world plant (maize) data sets. Our results demonstrate the ability of our approach to delineate divergent subpopulation-level behavior. Notably, our approach shows how environmental factors could influence phenotypic behavior, and how that effect varies across different genotypes and different time scales. To the best of our knowledge, this effort provides one of the first approaches to systematically formalize the problem of hypothesis extraction for phenomics data. Considering the infancy of the phenomics field, tools that help users explore complex data and extract plausible hypotheses in a data-guided manner will be critical to future advancements in the use of such data.


title: ALMA Maps of Dust and Warm Dense Gas Emission in the Starburst Galaxy IC 5179$^\star$

abstract: We present our high-resolution ($0^{\prime\prime}.15\times0^{\prime\prime}.13$, $\sim$34 pc) observations of the CO(6-5) line emission, which probes the warm and dense molecular gas, and the 434 $\mu$m dust continuum emission in the nuclear region of the starburst galaxy IC 5179, conducted with the Atacama Large Millimeter Array (ALMA). The CO(6-5) emission is spatially distributed in filamentary structures with many dense cores and shows a velocity field that is characteristic of a circum-nuclear rotating gas disk, with 90% of the rotation speed arising within a radius of $\lesssim150$ pc. At the scale of our spatial resolution, the CO(6-5) and dust emission peaks do not always coincide, with their surface brightness ratio varying by a factor of $\sim$10. This result suggests that their excitation mechanisms are likely different, as further evidenced by the Southwest to Northeast spatial gradient of both CO-to-dust continuum ratio and Pa-$\alpha$ equivalent width. Within the nuclear region (radius$\sim$300 pc) and with a resolution of $\sim$34 pc, the CO line flux (dust flux density) detected in our ALMA observations is $180\pm18$ Jy km/s ($71\pm7$ mJy), which account for 22% (2.4%) of the total value measured by Herschel.


title: Complex Event Processing of Health Data in Real-time to Predict Heart Failure Risk and Stress

abstract: In this paper, we develop a scalable system which can do real-time analytics for different health applications. The occurrence of different health conditions can be regarded as the complex events and thus this concept can be extended to other use cases easily. Large number of users should be able to send the data in real-time, and should be able to receive the feedback and result. Keeping the requirements in mind we used Kafka and Spark to develop our system. In this setting, multiple users are running Kafka producer clients, which are sending data in real-time. Spark streaming is used to process data from Kafka of different window sizes to analyze the health conditions. We have developed and tested the heart attack risk and stress prediction as our sample complex events detection use cases. We have simulated and tested our system with multiple health datasets.


title: Biological implications of $\mathfrak{so}(2,1)$ symmetry in exact solutions for a self-repressing gene

abstract: We chemically characterize the symmetries underlying the exact solutions of a stochastic negatively self-regulating gene. The breaking of symmetry at low molecular number causes three effects. Two branches of the solution exist, having high and low switching rates, such that the low switching rate branch approaches deterministic behavior and the high switching rate branch exhibits sub-Fano behavior. Average protein number differs from the deterministically expected value. Bimodal probability distributions appear as the protein number becomes a readout of the ON/OFF state of the gene.


title: Continuity of Hilbert-Kunz multiplicity and F-signature

abstract: We establish the continuity of Hilbert-Kunz multiplicity and F-signature as functions from a Cohen-Macaulay local ring $(R,\m,k)$ of prime characteristic to the real numbers at reduced parameter elements with respect to the $\m$-adic topology.


title: On Nematicity and Charge Order in Superoxygenated La$_{2-x}$Sr$_x$CuO$_{4+y}$

abstract: A nematic order in stripe-ordered cuprates was recently identified with (001) reflection at resonant energies associated with the in-plane states. However, whether this resonant reflection is ubiquitous among all 214 cuprates is still unknown. Here we report a Resonant soft X-ray Scattering (RXS) measurement on two La$_{2-x}$Sr$_x$CuO$_{4+y}$ crystals. Charge order was found in La$_2$CuO$_{4+y}$ with a total hole concentration near 0.125/Cu but no measurable (001) peak at any resonance, while in a La$_{1.94}$Sr$_{0.06}$CuO$_{4+y}$ sample near 0.16/Cu a (001) peak resonant was identified to be consistent with the presence of LTT tilting. The lack of such a (001) peak in a compound with stripe-like charge order raises questions about nematicity and the origin of the scattering feature.


title: Kernel Method for Detecting Higher Order Interactions in multi-view Data: An Application to Imaging, Genetics, and Epigenetics

abstract: In this study, we tested the interaction effect of multimodal datasets using a novel method called the kernel method for detecting higher order interactions among biologically relevant mulit-view data. Using a semiparametric method on a reproducing kernel Hilbert space (RKHS), we used a standard mixed-effects linear model and derived a score-based variance component statistic that tests for higher order interactions between multi-view data. The proposed method offers an intangible framework for the identification of higher order interaction effects (e.g., three way interaction) between genetics, brain imaging, and epigenetic data. Extensive numerical simulation studies were first conducted to evaluate the performance of this method. Finally, this method was evaluated using data from the Mind Clinical Imaging Consortium (MCIC) including single nucleotide polymorphism (SNP) data, functional magnetic resonance imaging (fMRI) scans, and deoxyribonucleic acid (DNA) methylation data, respectfully, in schizophrenia patients and healthy controls. We treated each gene-derived SNPs, region of interest (ROI) and gene-derived DNA methylation as a single testing unit, which are combined into triplets for evaluation. In addition, cardiovascular disease risk factors such as age, gender, and body mass index were assessed as covariates on hippocampal volume and compared between triplets. Our method identified $13$-triplets ($p$-values $\leq 0.001$) that included $6$ gene-derived SNPs, $10$ ROIs, and $6$ gene-derived DNA methylations that correlated with changes in hippocampal volume, suggesting that these triplets may be important in explaining schizophrenia-related neurodegeneration. With strong evidence ($p$-values $\leq 0.000001$), the triplet ({\bf MAGI2, CRBLCrus1.L, FBXO28}) has the potential to distinguish schizophrenia patients from the healthy control variations.