title: Translating Event-B machines to Eiffel programs

abstract: Formal modelling languages play a key role in the development of software since they enable users to prove correctness of system properties. However, there is still not a clear understanding on how to map a formal model to a specific programming language. In order to propose a solution, this paper presents a source-to-source mapping between Event- B models and Eiffel programs, therefore enabling the proof of correctness of certain system properties via Design-by-Contract (natively supported by Eiffel), while still making use of all features of O-O programming.


title: Limits of Predictability of Cascading Overload Failures in Spatially-Embedded Networks with Distributed Flows

abstract: Cascading failures are a critical vulnerability of complex information or infrastructure networks. Here we investigate the properties of load-based cascading failures in real and synthetic spatially-embedded network structures, and propose mitigation strategies to reduce the severity of damages caused by such failures. We introduce a stochastic method for optimal heterogeneous distribution of resources (node capacities) subject to a fixed total cost. Additionally, we design and compare the performance of networks with N-stable and (N-1)-stable network-capacity allocations by triggering cascades using various real-world node-attack and node-failure scenarios. We show that failure mitigation through increased node protection can be effectively achieved against single node failures. However, mitigating against multiple node failures is much more difficult due to the combinatorial increase in possible failures. We analyze the robustness of the system with increasing protection, and find that a critical tolerance exists at which the system undergoes a phase transition, and above which the network almost completely survives an attack. Moreover, we show that cascade-size distributions measured in this region exhibit a power-law decay. Finally, we find a strong correlation between cascade sizes induced by individual nodes and sets of nodes. We also show that network topology alone is a weak factor in determining the progression of cascading failures.

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


title: Context-Aware System Synthesis, Task Assignment, and Routing

abstract: The design and organization of complex robotic systems traditionally requires laborious trial-and-error processes to ensure both hardware and software components are correctly connected with the resources necessary for computation. This paper presents a novel generalization of the quadratic assignment and routing problem, introducing formalisms for selecting components and interconnections to synthesize a complete system capable of providing some user-defined functionality. By introducing mission context, functional requirements, and modularity directly into the assignment problem, we derive a solution where components are automatically selected and then organized into an optimal hardware and software interconnection structure, all while respecting restrictions on component viability and required functionality. The ability to generate \emph{complete} functional systems directly from individual components reduces manual design effort by allowing for a guided exploration of the design space. Additionally, our formulation increases resiliency by quantifying resource margins and enabling adaptation of system structure in response to changing environments, hardware or software failure. The proposed formulation is cast as an integer linear program which is provably $\mathcal{NP}$-hard. Two case studies are developed and analyzed to highlight the expressiveness and complexity of problems that can be addressed by this approach: the first explores the iterative development of a ground-based search-and-rescue robot in a variety of mission contexts, while the second explores the large-scale, complex design of a humanoid disaster robot for the DARPA Robotics Challenge. Numerical simulations quantify real world performance and demonstrate tractable time complexity for the scale of problems encountered in many modern robotic systems.


title: A new statistical method for characterizing the atmospheres of extrasolar planets

abstract: By detecting light from extrasolar planets,we can measure their compositions and bulk physical properties. The technologies used to make these measurements are still in their infancy, and a lack of self-consistency suggests that previous observations have underestimated their systemic errors.We demonstrate a statistical method, newly applied to exoplanet characterization, which uses a Bayesian formalism to account for underestimated errorbars. We use this method to compare photometry of a substellar companion, GJ 758b, with custom atmospheric models. Our method produces a probability distribution of atmospheric model parameters including temperature, gravity, cloud model (fsed), and chemical abundance for GJ 758b. This distribution is less sensitive to highly variant data, and appropriately reflects a greater uncertainty on parameter fits.


title: Existence versus Exploitation: The Opacity of Backbones and Backdoors Under a Weak Assumption

abstract: Backdoors and backbones of Boolean formulas are hidden structural properties. A natural goal, already in part realized, is that solver algorithms seek to obtain substantially better performance by exploiting these structures. However, the present paper is not intended to improve the performance of SAT solvers, but rather is a cautionary paper. In particular, the theme of this paper is that there is a potential chasm between the existence of such structures in the Boolean formula and being able to effectively exploit them. This does not mean that these structures are not useful to solvers. It does mean that one must be very careful not to assume that it is computationally easy to go from the existence of a structure to being able to get one's hands on it and/or being able to exploit the structure. For example, in this paper we show that, under the assumption that P $\neq$ NP, there are easily recognizable families of Boolean formulas with strong backdoors that are easy to find, yet for which it is hard (in fact, NP-complete) to determine whether the formulas are satisfiable. We also show that, also under the assumption P $\neq$ NP, there are easily recognizable sets of Boolean formulas for which it is hard (in fact, NP-complete) to determine whether they have a large backbone.


title: Cosmological constraints from a joint analysis of cosmic microwave background and spectroscopic tracers of the large-scale structure

abstract: The standard model of cosmology, {\Lambda}CDM, is the simplest model that matches the current observations, but it relies on two hypothetical components, to wit, dark matter and dark energy. Future galaxy surveys and cosmic microwave background (CMB) experiments will independently shed light on these components, but a joint analysis that includes cross-correlations will be necessary to extract as much information as possible from the observations. In this paper, we carry out a multi-probe analysis based on pseudo-spectra and test it on publicly available data sets. We use CMB temperature anisotropies and CMB lensing observations from Planck as well as the spectroscopic galaxy and quasar samples of SDSS-III/BOSS, taking advantage of the large areas covered by these surveys. We build a likelihood to simultaneously analyse the auto and cross spectra of CMB lensing and tracer overdensity maps before running Monte-Carlo Markov Chains (MCMC) to assess the constraining power of the combined analysis. We then add the CMB temperature anisotropies likelihood and obtain constraints on cosmological parameters ($H_0$, $\omega_b$, $\omega_c$, ${\ln10^{10}A_s}$, $n_s$ and $z_{re}$) and galaxy biases. We demonstrate that the joint analysis can additionally constrain the total mass of neutrinos ${\Sigma m_{\nu}}$ as well as the dark energy equation of state $w$ at once (for a total of eight cosmological parameters), which is impossible with either of the data sets considered separately. Finally, we discuss limitations of the analysis related to, e.g., the theoretical precision of the models, particularly in the non-linear regime.


title: Electric Currents at Semiconductor Surfaces from the Perspective of Drift-Diffusion Equations

abstract: Surface sensitive electric current measurements are important experimental tools poorly corroborated by theoretical models. We show that the drift-diffusion equations offer a framework for a consistent description of such experiments. The current flow is calculated as a perturbation of an equilibrium solution depicting the space charge layer. We investigate the accumulation and inversion layers in great detail. Relying on numerical findings, we identify the proper length parameter, the relationship of which with the length of the space charge layer is not simple. If the length parameter is large enough, long-ranged modes dominate the Green's function of the current equation, leading to two-dimensional currents. In addition, we demonstrate that the surface behavior of the currents is ruled by only a few parameters. This explains the fact that simplistic conductivity models have proven effective but makes reconstructions of conductance profiles from surface currents rather questionable.


title: High-Order Accurate FDTD Schemes for Dispersive Maxwell's Equations in Second-Order Form Using Recursive Convolutions

abstract: We propose a novel finite-difference time-domain (FDTD) scheme for the solution of the Maxwell's equations in which linear dispersive effects are present. The method uses high-order accurate approximations in space and time for the dispersive Maxwell's equations written as a second-order vector wave equation with a time-history convolution term. The modified equation approach is combined with the recursive convolution (RC) method to develop high-order approximations accurate to any desired order in space and time. High-order-accurate centered approximations of the physical Maxwell interface conditions are derived for the dispersive setting in order to fully restore accuracy at discontinuous material interfaces. Second- and fourth-order accurate versions of the scheme are presented and implemented in two spatial dimensions for the case of the Drude linear dispersion model. The stability of these schemes is analyzed. Finally, our approach is also amenable to curvilinear numerical grids if used with appropriate generalized Laplace operator.


title: Hyperscaling breakdown and Ising Spin Glasses: the Binder cumulant

abstract: Among the Renormalization Group Theory scaling rules relating critical exponents, there are hyperscaling rules involving the dimension of the system. It is well known that in Ising models hyperscaling breaks down above the upper critical dimension. It was shown by M. Schwartz [Europhys. Lett. {\bf 15}, 777 (1991)] that the standard Josephson hyperscaling rule can also break down in Ising systems with quenched random interactions. A related Renormalization Group Theory hyperscaling rule links the critical exponents for the normalized Binder cumulant and the correlation length in the thermodynamic limit. An appropriate scaling approach for analyzing measurements from criticality to infinite temperature is first outlined. Numerical data on the scaling of the normalized correlation length and the normalized Binder cumulant are shown for the canonical Ising ferromagnet model in dimension three where hyperscaling holds, for the Ising ferromagnet in dimension five (so above the upper critical dimension) where hyperscaling breaks down, and then for Ising spin glass models in dimension three where the quenched interactions are random. For the Ising spin glasses there is a breakdown of the normalized Binder cumulant hyperscaling relation in the thermodynamic limit regime, with a return to size independent Binder cumulant values in the finite-size scaling regime around the critical region.


title: Generalized scalar field models with the same energy density and linear stability

abstract: We study how the properties of a Lagrangian density for a single real scalar field in flat spacetime change with inclusion of an overall factor depending only on the field. The focus of the paper is to obtain analytical results. So, we show that even though it is possible to perform a field redefinition to get an equivalent canonical model, it is not always feasible to write the canonical model in terms of elementary functions. Also, we investigate the behavior of the energy density and the linear stability of the solutions. Finally, we show that one can find a class of models that present the same energy density and the same stability potential.


title: Contextual advantage for state discrimination

abstract: Finding quantitative aspects of quantum phenomena which cannot be explained by any classical model has foundational importance for understanding the boundary between classical and quantum theory. It also has practical significance for identifying information processing tasks for which those phenomena provide a quantum advantage. Using the framework of generalized noncontextuality as our notion of classicality, we find one such nonclassical feature within the phenomenology of quantum minimum error state discrimination. Namely, we identify quantitative limits on the success probability for minimum error state discrimination in any experiment described by a noncontextual ontological model. These constraints constitute noncontextuality inequalities that are violated by quantum theory, and this violation implies a quantum advantage for state discrimination relative to noncontextual models. Furthermore, our noncontextuality inequalities are robust to noise and are operationally formulated, so that any experimental violation of the inequalities is a witness of contextuality, independently of the validity of quantum theory. Along the way, we introduce new methods for analyzing noncontextuality scenarios, and demonstrate a tight connection between our minimum error state discrimination scenario and a Bell scenario.


title: Learning without Prejudice: Avoiding Bias in Webly-Supervised Action Recognition

abstract: Webly-supervised learning has recently emerged as an alternative paradigm to traditional supervised learning based on large-scale datasets with manual annotations. The key idea is that models such as CNNs can be learned from the noisy visual data available on the web. In this work we aim to exploit web data for video understanding tasks such as action recognition and detection. One of the main problems in webly-supervised learning is cleaning the noisy labeled data from the web. The state-of-the-art paradigm relies on training a first classifier on noisy data that is then used to clean the remaining dataset. Our key insight is that this procedure biases the second classifier towards samples that the first one understands. Here we train two independent CNNs, a RGB network on web images and video frames and a second network using temporal information from optical flow. We show that training the networks independently is vastly superior to selecting the frames for the flow classifier by using our RGB network. Moreover, we show benefits in enriching the training set with different data sources from heterogeneous public web databases. We demonstrate that our framework outperforms all other webly-supervised methods on two public benchmarks, UCF-101 and Thumos'14.


title: Upper bounds on secret key agreement over lossy thermal bosonic channels

abstract: Upper bounds on the secret-key-agreement capacity of a quantum channel serve as a way to assess the performance of practical quantum-key-distribution protocols conducted over that channel. In particular, if a protocol employs a quantum repeater, achieving secret-key rates exceeding these upper bounds is a witness to having a working quantum repeater. In this paper, we extend a recent advance [Liuzzo-Scorpo et al., arXiv:1705.03017] in the theory of the teleportation simulation of single-mode phase-insensitive Gaussian channels such that it now applies to the relative entropy of entanglement measure. As a consequence of this extension, we find tighter upper bounds on the non-asymptotic secret-key-agreement capacity of the lossy thermal bosonic channel than were previously known. The lossy thermal bosonic channel serves as a more realistic model of communication than the pure-loss bosonic channel, because it can model the effects of eavesdropper tampering and imperfect detectors. An implication of our result is that the previously known upper bounds on the secret-key-agreement capacity of the thermal channel are too pessimistic for the practical finite-size regime in which the channel is used a finite number of times, and so it should now be somewhat easier to witness a working quantum repeater when using secret-key-agreement capacity upper bounds as a benchmark.


title: PMU Assisted Power System Parameter Calibration at Jiangsu Electric Power Company

abstract: An online PMU-assisted Power System Parameter Calibration System (PSPCS) was recently developed and implemented at State Grid Jiangsu Electric Power Company (JEPC). PSPCS leverages high-resolution PMU data and data mining techniques to perform online screening of the EMS and Production Management System (PMS) databases for data cleaning, model validation, and parameter calibration. PSPCS calculates transmission line and generator parameters on a regular real-time basis and compares the results with databases to identify record(s) with significant discrepancy, if any. Once consistent discrepancy is observed, the system will raise a flag and further investigation will be initiated, including a novel density-based spatial clustering procedure for parameter/data calibration. A novel metric is proposed to quantify the credibility of PMU-based parameter identification. This paper discusses the proposed methodologies, challenges, as well as implementation issues identified during the development and deployment of PSPCS.


title: Search for active-sterile neutrino mixing using neutral-current interactions in NOvA

abstract: We report results from the first search for sterile neutrinos mixing with active neutrinos through a reduction in the rate of neutral-current interactions over a baseline of 810\,km between the NOvA detectors. Analyzing a 14-kton detector equivalent exposure of 6.05$\times$10$^{20}$ protons-on-target in the NuMI beam at Fermilab, we observe 95 neutral-current candidates at the Far Detector compared with $83.5 \pm 9.7 \mbox{(stat.)} \pm 9.4 \mbox{(syst.)}$ events predicted assuming mixing only occurs between active neutrino species. No evidence for $\nu_{\mu} \rightarrow \nu_{s}$ transitions is found. Interpreting these results within a 3+1 model, we place constraints on the mixing angles $\theta_{24}<20.8^{\circ}$ and $\theta_{34}<31.2^{\circ}$ at the 90% C.L. for $0.05~eV^2\leq \Delta m^2_{41}\leq 0.5~eV^2$, the range of mass splittings that produce no significant oscillations over the Near Detector baseline.


title: Perturbed moments and a longer mollifier for critical zeros of $\zeta$

abstract: Let $A(s)$ be a general Dirichlet polynomial and $\Phi$ be a smooth function supported in $[1,2]$ with mild bounds on its derivatives. New main terms for the integral $I(\alpha,\beta)=\int_{\mathbb{R}} \zeta(\frac{1}{2}+\alpha+it)\zeta(\frac{1}{2}+\beta+it)|A(\frac{1}{2}+it)|^2 \Phi(\frac{t}{T})dt$ are given. For the error term, we show that the length of the Feng mollifier can be increased from $\theta < \frac{17}{33}$ to $\theta < \frac{6}{11}$ by decomposing the error into Type I and Type II sums and then studying the resulting sums of Kloosterman sums. As an application, we slightly increase the proportion of zeros of $\zeta(s)$ on the critical line.


title: Lyapunov inequality for a boundary value problem involving conformable derivative

abstract: We consider a boundary value problem involving conformable derivative of order $\alpha ,$ $1<\alpha <2$ and Dirichlet conditions. To prove the existence of solutions, we apply the method of upper and lower solutions together with Schauder's fixed-point theorem. \ Futhermore, we give the Lyapunov inequality for the corresponding problem.


title: Security Camera Movie and ERP Data Matching System to Prevent Theft

abstract: In this paper, we propose a SaaS service which prevents shoplifting using image analysis and ERP. In Japan, total damage of shoplifting reaches 450 billion yen. Based on cloud and data analysis technology, we propose a shoplifting prevention service with image analysis of security camera and ERP data check for small shops. We evaluated movie analysis.


title: Sequential Optimal Placement of Distributed Photovoltaics using Downstream Power Index

abstract: The optimization of the size and location of Photovoltaic (PV) Distributed Generation (DG) is a method for reducing distribution networks loss, cutting costs for utilities, and integrating renewable energy into the power grid. However, this optimization problem is a difficult mixed continuous discrete problem that is difficult to solve efficiently. Herein, we propose a novel Sequential PVDG Placement algorithm which utilizes a Downstream Power Index (DPI) for PV allocation. We compare Sequential PVDG Placement results with those of Shuffled Frog Leaping Algorithm (SFLA), a memetic heuristic algorithm. Our algorithm demonstrates similar accuracy and drastically less computation time compared to SFLA. Given its high accuracy and low computation time, sequential PVDG allocation algorithm may be useful for rapid online power dispatch, long-term planning, and microgrid operations under high penetration of renewable sources.


title: Computing equilibrium states of cholesteric liquid crystals in elliptical channels with deflation algorithms

abstract: We study the problem of a cholesteric liquid crystal confined to an elliptical channel. The system is geometrically frustrated because the cholesteric prefers to adopt a uniform rate of twist deformation, but the elliptical domain precludes this. The frustration is resolved by deformation of the layers or introduction of defects, leading to a particularly rich family of equilibrium configurations. To identify the solution set, we adapt and apply a new family of algorithms, known as deflation methods, that iteratively modify the free energy extremisation problem by removing previously known solutions. A second algorithm, deflated continuation, is used to track solution branches as a function of the aspect ratio of the ellipse and preferred pitch of the cholesteric.


title: Discrepancies in simultaneous explanation of Flavor Anomalies and IceCube PeV Events using Leptoquarks

abstract: Leptoquarks have been suggested to solve a variety of discrepancies between the expected and observed phenomenon. In this paper, we show that the scalar doublet Leptoquark with Hypercharge 7/6 can simultaneously explain the recent measurement of $R_{K}$, $R_{K^*}$, the excess in anomalous magnetic moment of muon, and the observed excess in IceCube HESE data. For appropriate choice of couplings, the flavor anomalies are generated at one-loop level and IceCube data is explained via resonant production of the Leptoquark. Several constraints from LHC searches are imposed on the model parameter space.


title: On Calibration of Modern Neural Networks

abstract: Confidence calibration -- the problem of predicting probability estimates representative of the true correctness likelihood -- is important for classification models in many applications. We discover that modern neural networks, unlike those from a decade ago, are poorly calibrated. Through extensive experiments, we observe that depth, width, weight decay, and Batch Normalization are important factors influencing calibration. We evaluate the performance of various post-processing calibration methods on state-of-the-art architectures with image and document classification datasets. Our analysis and experiments not only offer insights into neural network learning, but also provide a simple and straightforward recipe for practical settings: on most datasets, temperature scaling -- a single-parameter variant of Platt Scaling -- is surprisingly effective at calibrating predictions.


title: Large Fermi Arcs in Unconventional Weyl Semimetal RhSi

abstract: The theoretical proposal of chiral fermions in topological semimetals has led to a significant effort towards their experimental realization. In particular, the Fermi surfaces of chiral semimetals carry quantized Chern numbers, making them an attractive platform for the observation of exotic transport and optical phenomena. While the simplest example of a chiral fermion in condensed matter is a conventional $|C|=1$ Weyl fermion, recent theoretical works have proposed a number of unconventional chiral fermions beyond the Standard Model which are protected by unique combinations of topology and crystalline symmetries. However, materials candidates for experimentally probing the transport and response signatures of these unconventional fermions have thus far remained elusive. In this paper, we propose the RhSi family in space group (SG) $\#$198 as the ideal platform for the experimental examination of unconventional chiral fermions. We find that RhSi is a filling-enforced semimetal that features near its Fermi surface a chiral double six-fold-degenerate spin-1 Weyl node at $R$ and a previously uncharacterized four-fold-degenerate chiral fermion at $\Gamma$. Each unconventional fermion displays Chern number $\pm4$ at the Fermi level. We also show that RhSi displays the largest possible momentum separation of compensative chiral fermions, the largest proposed topologically nontrivial energy window, and the longest possible Fermi arcs on its surface. We conclude by proposing signatures of an exotic bulk photogalvanic response in RhSi.


title: Provable benefits of representation learning

abstract: There is general consensus that learning representations is useful for a variety of reasons, e.g. efficient use of labeled data (semi-supervised learning), transfer learning and understanding hidden structure of data. Popular techniques for representation learning include clustering, manifold learning, kernel-learning, autoencoders, Boltzmann machines, etc. To study the relative merits of these techniques, it's essential to formalize the definition and goals of representation learning, so that they are all become instances of the same definition. This paper introduces such a formal framework that also formalizes the utility of learning the representation. It is related to previous Bayesian notions, but with some new twists. We show the usefulness of our framework by exhibiting simple and natural settings -- linear mixture models and loglinear models, where the power of representation learning can be formally shown. In these examples, representation learning can be performed provably and efficiently under plausible assumptions (despite being NP-hard), and furthermore: (i) it greatly reduces the need for labeled data (semi-supervised learning) and (ii) it allows solving classification tasks when simpler approaches like nearest neighbors require too much data (iii) it is more powerful than manifold learning methods.


title: A local to global argument on low dimensional manifolds

abstract: For an oriented manifold $M$ whose dimension is less than $4$, we use the contractibility of certain complexes associated to its submanifolds to cut $M$ into simpler pieces in order to do local to global arguments. In particular, in these dimensions, we give a different proof of a deep theorem of Thurston in foliation theory which says that the natural map between classifying spaces $\mathrm{B}\text{Homeo}^{\delta}(M)\to \mathrm{B}\text{Homeo}(M)$ induces a homology isomorphism where $\text{Homeo}^{\delta}(M)$ denotes the group of homeomorphisms of $M$ made discrete. Our proof shows that in low dimensions, Thurston's theorem can be proved without using foliation theory. Finally, we show that this technique gives a new perspective on the homotopy type of homeomorphism groups in low dimensions. In particular, we give a different proof of Hacher's theorem that the homeomorphism groups of Haken $3$-manifolds with boundary are homotopically discrete without using his disjunction techniques.


title: A change in stripes for cholesteric shells via anchoring in moderation

abstract: Chirality, ubiquitous in complex biological systems, can be controlled and quantified in synthetic materials such as cholesteric liquid crystal (CLC) systems. In this work, we study spherical shells of CLC under weak anchoring conditions. We induce anchoring transitions at the inner and outer boundaries using two independent methods: by changing the surfactant concentration or by raising the temperature close to the clearing point. The shell confinement leads to new states and associated surface structures: a state where large stripes on the shell can be filled with smaller, perpendicular sub-stripes, and a focal conic domain (FCD) state, where thin stripes wrap into at least two, topologically required, double spirals. Focusing on the latter state, we use a Landau-de Gennes model of the CLC to simulate its detailed configurations as a function of anchoring strength. By abruptly changing the topological constraints on the shell, we are able to study the interconversion between director defects and pitch defects, a phenomenon usually restricted by the complexity of the cholesteric phase. This work extends the knowledge of cholesteric patterns, structures that not only have potential for use as intricate, self-assembly blueprints but are pervasive in biological systems.


title: High-rank ternary forms of even degree

abstract: We exhibit, for each even degree, a ternary form of rank strictly greater than the maximum rank of monomials. Together with an earlier result in the odd case, this gives the lower bound \[\operatorname{r_{max}}(3,d)\ge\left\lfloor\frac{d^2+2d+5}4\right\rfloor\] for $d\ge 2$, where $\operatorname{r_{max}}(n,d)$ denotes the maximum rank of degree $d$ forms in $n$ variables with coefficients in an algebraically closed field of characteristic zero.


title: On current contribution to Fronsdal equations

abstract: We explore a local form of second-order Vasiliev equations proposed in [arXiv:1706.03718] and obtain an explicit expression for quadratic corrections to bosonic Fronsdal equations, generated by gauge-invariant higher-spin currents. Our analysis is performed for general phase factor, and for the case of parity-invariant theory we find the agreement with expressions for cubic vertices available in the literature. This provides an additional indication that field redefinition proposed in [arXiv:1706.03718] is the proper one.


title: Nudged elastic band calculations accelerated with Gaussian process regression

abstract: Minimum energy paths for transitions such as atomic and/or spin rearrangements in thermalized systems are the transition paths of largest statistical weight. Such paths are frequently calculated using the nudged elastic band method, where an initial path is iteratively shifted to the nearest minimum energy path. The computational effort can be large, especially when ab initio or electron density functional calculations are used to evaluate the energy and atomic forces. Here, we show how the number of such evaluations can be reduced by an order of magnitude using a Gaussian process regression approach where an approximate energy surface is generated and refined in each iteration. When the goal is to evaluate the transition rate within harmonic transition state theory, the evaluation of the Hessian matrix at the initial and final state minima can be carried out beforehand and used as input in the minimum energy path calculation, thereby improving stability and reducing the number of iterations needed for convergence. A Gaussian process model also provides an uncertainty estimate for the approximate energy surface, and this can be used to focus the calculations on the lesser-known part of the path, thereby reducing the number of needed energy and force evaluations to a half in the present calculations. The methodology is illustrated using the two-dimensional M\"uller-Brown potential surface and performance assessed on an established benchmark involving 13 rearrangement transitions of a heptamer island on a solid surface.

categories: physics.chem-ph physics.atm-clus physics.comp-ph stat.CO stat.ML


title: Boosting the Terahertz Photoconductive Antenna Performance with Optimized Plasmonic Nanostructures

abstract: Advanced nanophotonics penetrates into other areas of science and technology, ranging from applied physics to biology and resulting in many fascinating cross-disciplinary applications. It has been recently demonstrated that suitably engineered light-matter interactions at the nanoscale can overcome the limitations of today's terahertz (THz) photoconductive antennas, making them one step closer to many practical implications. Here we push forward this concept by comprehensive numerical optimization and experimental investigation of a log-periodic THz photoconductive antenna coupled to a silver nanoantenna array. We shed light on the operation principles of the resulting hybrid THz antenna, providing an approach to boost its performance. By tailoring the size of silver nanoantennas and the distance between them, we obtain an enhancement of optical-to-THz conversion efficiency 2-fold larger compared with previously reported results, and the strongest enhancement is around 1 THz, a frequency range barely achievable by other compact THz sources. Moreover, we propose a cost-effective fabrication procedure to realize such hybrid THz antennas with optimized plasmonic nanostructures via thermal dewetting process, which does not require any post processing and makes the proposed solution very attractive for applications.


title: Co-axial monodromy

abstract: For Riemannian metrics of constant positive curvature on a punctured sphere with conic singularities at the punctures and co-axial monodromy of the developing map, possible angles at the singularities are completely described. This completes the recent result of Mondello and Panov. The related problem of describing possible multiplicities of critical points of logarithmic potentials of finitely many charges is also solved.


title: Disturbed Fossil Group Galaxy NGC 1132

abstract: We have analyzed the Chandra archival data of NGC 1132, a well-known fossil group, i.e. a system expected to be old and relaxed long after the giant elliptical galaxy assembly. Instead, the Chandra data reveal that the hot gas morphology is disturbed and asymmetrical, with a cold front following a possible bow shock. We discuss possible origins of the disturbed hot halo, including sloshing by a nearby object, merger, ram pressure by external hotter gas and nuclear outburst. We consider that the first two mechanisms are likely explanations for the disturbed hot halo, with a slight preference for a minor merger with a low impact parameter because of the match with simulations and previous optical observations. In this case, NGC 1132 may be a rare example of unusual late mergers seen in recent simulations. Regardless of the origin of the disturbed hot halo, the paradigm of the fossil system needs to be reconsidered.


title: Transverse contributions to the longitudinal stiffness of the human foot

abstract: Humans rely on the foot's stiffness to withstand the propulsive forces of walking and running. The medial longitudinal arch (MLA) and the transverse tarsal arch (TTA) are key adaptations that increase the human midfoot's longitudinal bending stiffness. The TTA is hypothesized to increase stiffness through a cross-axis coupling between transverse tissue stiffness and longitudinal midfoot stiffness. However, we currently lack in vivo tests of this hypothesis. Such tests are crucial because muscle contraction considerably modulates MLA function and may also affect the TTA's predicted cross-axis coupling. Here we show in vivo evidence for substantial cross-axis coupling in standing human subjects. We found that longitudinal midfoot bending stiffness increased upon increasing transverse forefoot stiffness using an elastic tape (mean$\pm$standard error = $76\% \pm 19\%$, p=0.0002), thus supporting the TTA hypothesis. The cross-axis coupling could be a potential therapeutic target for some foot disorders and help augment foot stiffness for athletics.


title: On the origin of the hard X-Ray excess of high-synchrotron-peaked BL Lac object Mrk 421

abstract: For the first time, Kataoka \& Stawarz reported a clear detection of a hard X-ray excess, above $\gtrsim$20 keV, in the high-synchrotron-peaked BL Lac object Mrk 421. We find that this feature may not be produced by the low-energy part of the same electron population that produced the Fermi/LAT $\gamma$-ray. Because of that it is required that the power-law electron energy go down to $\gamma_{\rm min}\approx19$, which predicts a very strong radio emission (radio flux larger than the observed) even considering the synchrotron self-absorption effect. We investigate the possibility of this excess being produced from the spine/layer jet structure, which has been clearly detected in Mrk 421. We find that (1) similar to one-zone modeling, the spine emissions provide good modeling of the broadband spectral energy distribution, except for the hard X-ray excess; and (2) the hard X-ray excess can be well represented by the synchrotron photons (from the layer) being inverse Compton scattered by the spine electrons.


title: Charged structure constants from modularity

abstract: We derive a universal formula for the average heavy-heavy-light structure constants for 2d CFTs with non-vanishing u(1) charge. The derivation utilizes the modular properties of one-point functions on the torus. Refinements in N=2 SCFTs, show that the resulting Cardy-like formula for the structure constants has precisely the same shifts in the central charge as that of the thermodynamic entropy found earlier. This analysis generalizes the recent results by Kraus and Maloney for CFTs with an additional global u(1) symmetry. Our results at large central charge are also shown to match with computations from the holographic dual, which suggest that the averaged CFT three-point coefficient also serves as an useful probe of detecting black hole hair.


title: HST imaging of the brightest z~8-9 galaxies from UltraVISTA: the extreme bright end of the UV luminosity function

abstract: We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-$K_{S} \sim 25.3-24.8$~mag, $5\sigma$) NIR data from the UltraVISTA DR3 release, deep ground based optical imaging from the CFHTLS and Subaru Suprime Cam programs and Spitzer/IRAC mosaics combining observations from the SMUVS and SPLASH programs. Through the refined spectral energy distributions, which now also include new HyperSuprime Cam g, r, i, z and Y band data, we confirm that 3/5 galaxies have robust $z_{phot}\sim8.0-8.7$, consistent with the initial selection. The remaining 2/5 galaxies have a nominal $z_{phot}\sim2$. However, if we use the HST data alone, these objects have increased probability of being at $z\sim9$. Furthermore, we measure mean UV continuum slopes $\beta=-1.91\pm0.26$ for the three $z\sim8-9$ galaxies, marginally bluer than similarly luminous $z\sim4-6$ in CANDELS but consistent with previous measurements of similarly luminous galaxies at $z\sim7$. The circularized effective radius for our brightest source is $0.9\pm0.2$ kpc, similar to previous measurements for a bright $z\sim11$ galaxy and bright $z\sim7$ galaxies. Finally, enlarging our sample to include the six brightest $z\sim8$ LBGs identified over UltraVISTA (i.e., including three other sources from Labbe et al. 2017, in prep.) we estimate for the first time the volume density of galaxies at the extreme bright ($M_{UV}\sim-22$~mag) end of the $z\sim8$ UV LF. Despite this exceptional result, the still large statistical uncertainties do not allow us to discriminate between a Schechter and a double power-law form.


title: Rotation in [CII]-emitting gas in two galaxies at a redshift of 6.8

abstract: The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang.


title: Quadratic genetic modifications: a streamlined route to cosmological simulations with controlled merger history

abstract: Recent work has studied the interplay between a galaxy's history and its observable properties using "genetically modified" cosmological zoom simulations. The approach systematically generates alternative histories for a halo, while keeping its cosmological environment fixed. Applications to date altered linear properties of the initial conditions such as the mean overdensity of specified regions; we extend the formulation to include quadratic features such as local variance, which determines the overall importance of smooth accretion relative to mergers in a galaxy's history. We introduce an efficient algorithm for this new class of modification and demonstrate its ability to control the variance of a region in a one-dimensional toy model. Outcomes of this work are two-fold: (i) a clarification of the formulation of genetic modifications and (ii) a proof of concept for quadratic modifications leading the way to a forthcoming implementation in cosmological simulations.


title: Gauge Backgrounds and Zero-Mode Counting in F-Theory

abstract: Computing the exact spectrum of charged massless matter is a crucial step towards understanding the effective field theory describing F-theory vacua in four dimensions. In this work we further develop a coherent framework to determine the charged massless matter in F-theory compactified on elliptic fourfolds, and demonstrate its application in a concrete example. The gauge background is represented, via duality with M-theory, by algebraic cycles modulo rational equivalence. Intersection theory within the Chow ring allows us to extract coherent sheaves on the base of the elliptic fibration whose cohomology groups encode the charged zero-mode spectrum. The dimensions of these cohomology groups are computed with the help of modern techniques from algebraic geometry, which we implement in the software gap. We exemplify this approach in models with an Abelian and non-Abelian gauge group and observe jumps in the exact massless spectrum as the complex structure moduli are varied. An extended mathematical appendix gives a self-contained introduction to the algebro-geometric concepts underlying our framework.


title: The Threatening Environment of the TRAPPIST-1 Planets

abstract: Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind and energetic electromagnetic and particle radiation. Here, we use state-of-the-art 3D magnetohydrodynamic models to simulate the wind around TRAPPIST-1 and study the conditions at each planetary orbit. All planets experience a stellar wind pressure between $10^3$ and $10^5$ times the solar wind pressure on Earth. All orbits pass through wind pressure changes of an order of magnitude and most planets spend a large fraction of their orbital period in the sub-Alfv\'enic regime. For plausible planetary magnetic field strengths, all magnetospheres are greatly compressed and undergo much more dynamic change than that of the Earth. The planetary magnetic fields connect with the stellar radial field over much of the planetary surface, allowing direct flow of stellar wind particles onto the planetary atmosphere. These conditions could result in strong atmospheric stripping and evaporation and should be taken into account for any realistic assessment of the evolution and habitability of the TRAPPIST-1 planets.


title: Fantastic Striations and Where to Find Them: Origin of Magnetically Aligned Striations in Interstellar Clouds

abstract: Thin, magnetically aligned striations of relatively moderate contrast with the background are commonly observed in both atomic and molecular clouds. They are also prominent in MHD simulations with turbulent converging shocks. The simulated striations develop within a dense, stagnated sheet in the mid plane of the post-shock region where magnetically-induced converging flows collide. We show analytically that the secondary flows are an inevitable consequence of the jump conditions of oblique MHD shocks. They produce the stagnated, sheet-like sub-layer through a secondary shock when, roughly speaking, the Alfv\'enic speed in the primary converging flows is supersonic, a condition that is relatively easy to satisfy in interstellar clouds. The dense sub-layer is naturally threaded by a strong magnetic field that lies close to the plane of the sub-layer. The substantial magnetic field makes the sheet highly anisotropic, which is the key to the striation formation. Specifically, perturbations of the primary inflow that vary spatially perpendicular to the magnetic field can easily roll up the sheet around the field lines without bending them, creating corrugations that appear as magnetically-aligned striations in column density maps. On the other hand, perturbations that vary spatially along the field lines curve the sub-layer and alter its orientation relative to the magnetic field locally, seeding special locations that become slanted overdense filaments and prestellar cores through enhanced mass accumulation along field lines. In our scenario, the dense sub-layer unique to magnetized oblique shocks is the birthplace for both magnetically-aligned diffuse striations and massive star-forming structures.


title: E_n Jacobi forms and Seiberg-Witten curves

abstract: We discuss Jacobi forms that are invariant under the action of the Weyl group of type E_n (n=6,7,8). For n=6,7 we explicitly construct a full set of generators of the algebra of E_n weak Jacobi forms. We first construct n+1 independent E_n Jacobi forms in terms of Jacobi theta functions and modular forms. By using them we obtain Seiberg-Witten curves of type E_6 and E_7 for the E-string theory. The coefficients of each curve are E_n weak Jacobi forms of particular weights and indices specified by the root system, realizing the generators whose existence was shown some time ago by Wirthm\"uller.


title: No Ly$\alpha$ emitters detected around a QSO at z=6.4: Suppressed by the QSO?

abstract: Understanding how QSO's UV radiation affects galaxy formation is vital to our understanding of reionization era. Using a custom made narrow-band filter, $NB906$, on Subaru/Suprime-Cam, we investigated the number density of Ly$\alpha$ emitters (LAE) around a QSO at z=6.4. To date, this is the highest redshift narrow-band observation, where LAEs around a luminous QSO are investigated. Due to the large field-of-view of Suprime-Cam, our survey area is $\sim$5400~cMpc$^2$, much larger than previously studies at z=5.7 ($\sim$200 cMpc$^2$). In this field, we previously found a factor of 7 overdensity of Lyman break galaxies (LBGs). Based on this, we expected to detect $\sim$100 LAEs down to $NB906$=25 ABmag. However, our 6.4 hour exposure found none. The obtained upper limit on the number density of LAEs is more than an order lower than the blank fields. Furthermore, this lower density of LAEs spans a large scale of 10 $p$Mpc across. A simple argument suggests a strong UV radiation from the QSO can suppress star-formation in halos with $M_{vir}<10^{10}M_{\odot}$ within a $p$Mpc from the QSO, but the deficit at the edge of the field (5 $p$Mpc) remains to be explained.


title: Kinetic theory of transport for inhomogeneous electron fluids

abstract: The interplay between electronic interactions and disorder is neglected in the conventional Boltzmann theory of transport, yet can play an essential role in determining the resistivity of unconventional metals. When quasiparticles are long-lived, one can account for these intertwined effects by solving spatially inhomogeneous Boltzmann equations. Assuming smooth disorder and neglecting umklapp scattering, we solve these inhomogeneous kinetic equations and compute the electrical resistivity across the ballistic-to-hydrodynamic transition. An important consequence of electron-electron interactions is the modification of the momentum relaxation time; this effect is ignored in the conventional theory. We characterize precisely when interactions enhance the momentum scattering rate, and when they decrease it. Our approach unifies existing semiclassical theories of transport and reveals novel transport mechanisms. In particular, we explain how the resistivity can be proportional to the rate of momentum-conserving collisions. We compare this result with existing transport mysteries, including the disorder-independent $T^2$ resistivity of many Fermi liquids, and the linear-in-$T$ "Planckian-limited" resistivity of many strange metals.


title: Short-distance matrix elements for $D^0$-meson mixing for $N_f=2+1$ lattice QCD

abstract: We calculate in three-flavor lattice QCD the short-distance hadronic matrix elements of all five $\Delta C=2$ four-fermion operators that contribute to neutral $D$-meson mixing both in and beyond the Standard Model. We use the MILC Collaboration's $N_f = 2+1$ lattice gauge-field configurations generated with asqtad-improved staggered sea quarks. We also employ the asqtad action for the valence light quarks and use the clover action with the Fermilab interpretation for the charm quark. We analyze a large set of ensembles with pions as light as $M_\pi \approx 180$ MeV and lattice spacings as fine as $a\approx 0.045$ fm, thereby enabling good control over the extrapolation to the physical pion mass and continuum limit. We obtain for the matrix elements in the $\overline{\text{MS}}$-NDR scheme using the choice of evanescent operators proposed by Beneke \emph{et al.}, evaluated at 3 GeV, $\langle D^0|\mathcal{O}_i|\bar{D}^0 \rangle = \{0.0805(55)(16), -0.1561(70)(31), 0.0464(31)(9), 0.2747(129)(55), 0.1035(71)(21)\}~\text{GeV}^4$ ($i=1$--5). The errors shown are from statistics and lattice systematics, and the omission of charmed sea quarks, respectively. To illustrate the utility of our matrix-element results, we place bounds on the scale of CP-violating new physics in $D^0$~mixing, finding lower limits of about 10--50$\times 10^3$ TeV for couplings of $\mathrm{O}(1)$. To enable our results to be employed in more sophisticated or model-specific phenomenological studies, we provide the correlations among our matrix-element results. For convenience, we also present numerical results in the other commonly-used scheme of Buras, Misiak, and Urban.


title: Thermal Jeans fragmentation within 1000 AU in OMC-1S

abstract: We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density of Companions method we find a characteristic spatial scale at ~560 AU, and we use this spatial scale to define the boundary of 19 `cores' in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ~2900 AU, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation towards the southern filament. In addition, the cores of the southern filament are also the densest (within 1100 AU) cores in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 AU seems to be governed by thermal Jeans processes in OMC-1S.


title: The First Scattered Light Image of the Debris Disk around the Sco-Cen target HD 129590

abstract: We present the first scattered light image of the debris disk around HD 129590, a ~1.3 M$_\odot$ G1V member of the Scorpius Centaurus association with age ~10-16 Myr. The debris disk is imaged with the high contrast imaging instrument SPHERE at the Very Large Telescope, and is revealed by both the IRDIS and IFS subsytems, operating in the H and YJ bands respectively. The disk has a high infrared luminosity of $L_{\textrm{IR}}/L_{\textrm{star}}$~5$\times$10$^{-3}$, and has been resolved in other studies using ALMA. We detect a nearly edge on ring, with evidence of an inner clearing. We fit the debris disk using a model characterized by a single bright ring, with radius ~60-70 AU, in broad agreement with previous analysis of the target SED. The disk is vertically thin, and has an inclination angle of ~75$^\circ$. Along with other previously imaged edge-on disks in the Sco-Cen association such as HD 110058, HD 115600, and HD 111520, this disk image will allow of the structure and morphology of very young debris disks, shortly after the epoch of planet formation has ceased.


title: Properties of soliton surfaces associated with integrable $\mathbb{C}P^{N-1}$ sigma models

abstract: We investigate certain properties of $\mathfrak{su}(N)$-valued two-dimensional soliton surfaces associated with the integrable $\mathbb{C}P^{N-1}$ sigma models constructed by the orthogonal rank-one Hermitian projectors, which are defined on the two-dimensional Riemann sphere with finite action functional. Several new properties of the projectors mapping onto one-dimensional subspaces as well as their relations with three mutually different immersion formulas, namely, the generalized Weierstrass, Sym-Tafel and Fokas-Gel'fand have been discussed in detail. Explicit connections among these three surfaces are also established by purely analytical descriptions and, it is demonstrated that the three immersion formulas actually correspond to the single surface parametrized by some specific conditions.


title: Estimation and Mitigation of Channel Non-Reciprocity in Massive MIMO

abstract: Time-division duplex (TDD) based massive MIMO systems rely on the reciprocity of the wireless propagation channels when calculating the downlink precoders based on uplink pilots. However, the effective uplink and downlink channels incorporating the analog radio front-ends of the base station (BS) and user equipments (UEs) exhibit non-reciprocity due to non-identical behavior of the individual transmit and receive chains. When downlink precoder is not aware of such channel non-reciprocity (NRC), system performance can be significantly degraded due to NRC induced interference terms. In this work, we consider a general TDD-based massive MIMO system where frequency-response mismatches at both the BS and UEs, as well as the mutual coupling mismatch at the BS large-array system all coexist and induce channel NRC. Based on the NRC-impaired signal models, we first propose a novel iterative estimation method for acquiring both the BS and UE side NRC matrices and then also propose a novel NRC-aware downlink precoder design which utilizes the obtained estimates. Furthermore, an efficient pilot signaling scheme between the BS and UEs is introduced in order to facilitate executing the proposed estimation method and the NRC-aware precoding technique in practical systems. Comprehensive numerical results indicate substantially improved spectral efficiency performance when the proposed NRC estimation and NRC-aware precoding methods are adopted, compared to the existing state-of-the-art methods.


title: Parity breaking signatures from a Chern-Simons coupling during inflation: the case of non-Gaussian gravitational waves

abstract: Considering high-energy modifications of Einstein gravity during inflation is an interesting issue. We can constrain the strength of the new gravitational terms through observations of inflationary imprints in the actual universe. In this paper we analyze the effects on slow-roll models due to a Chern-Simons term coupled to the inflaton field through a generic coupling function $f(\phi)$. A well known result is the polarization of primordial gravitational waves (PGW) into left and right eigenstates, as a consequence of parity breaking. In such a scenario the modifications to the power spectrum of PGW are suppressed under the conditions that allow to avoid the production of ghost gravitons at a certain energy scale, the so-called Chern-Simons mass $M_{CS}$. In general it has been recently pointed out that there is very little hope to efficiently constrain chirality of PGW on the basis solely of two-point statistics from future CMB data, even in the most optimistic cases. Thus we search if significant parity breaking signatures can arise at least in the bispectrum statistics. We find that the tensor-tensor-scalar bispectra $\langle \gamma \gamma \zeta \rangle$ for each polarization state are the only ones that are not suppressed. Their amplitude, setting the level of parity breaking during inflation, is proportional to the second derivative of the coupling function $f(\phi)$ and they turn out to be maximum in the squeezed limit. We comment on the squeezed-limit consistency relation arising in the case of chiral gravitational waves, and on possible observables to constrain these signatures.


title: Simple and explicit bounds for multi-server queues with $1/(1-\rho)$ scaling

abstract: We consider the FCFS $GI/GI/n$ queue, and prove the first simple and explicit bounds that scale as $\frac{1}{1-\rho}$ under only the assumption that inter-arrival times have finite second moment, and service times have finite $2+\epsilon$ moment for some $\epsilon > 0$. Here $\rho$ denotes the corresponding traffic intensity. Conceptually, our results can be viewed as a multi-server analogue of Kingman's bound. Our main results are bounds for the tail of the steady-state queue length and the steady-state probability of delay. The strength of our bounds (e.g. in the form of tail decay rate) is a function of how many moments of the service distribution are assumed finite. Our bounds scale gracefully even when the number of servers grows large and the traffic intensity converges to unity simultaneously, as in the Halfin-Whitt scaling regime. Some of our bounds scale better than $\frac{1}{1-\rho}$ in certain asymptotic regimes. In these same asymptotic regimes we also prove bounds for the tail of the steady-state number in service. Our main proofs proceed by explicitly analyzing the bounding process which arises in the stochastic comparison bounds of Gamarnik and Goldberg for multi-server queues. Along the way we derive several novel results for suprema of random walks and pooled renewal processes which may be of independent interest. We also prove several additional bounds using drift arguments (which have much smaller pre-factors), and point out a conjecture which would imply further related bounds and generalizations. We also show that when all moments of the service distribution are finite and satisfy a mild growth rate assumption, our bounds can be strengthened to yield explicit tail estimates decaying as $O\big(\exp(-x^{\alpha})\big)$, with $\alpha \in (0,1)$ depending on the growth rate of these moments.


title: Verification Studies for the Noh Problem using Non-ideal Equations of State and Finite Strength Shocks

abstract: The Noh verification test problem is extended beyond the commonly studied ideal gamma-law gas to more realistic equations of state (EOSs) including the stiff gas, the Noble-Abel gas, and the Carnahan-Starling EOS for hard-sphere fluids. Self-similarity methods are used to solve the Euler compressible flow equations, which in combination with the Rankine-Hugoniot jump conditions provide a tractable general solution. This solution can be applied to fluids with EOSs that meet criterion such as it being a convex function and having a corresponding bulk modulus. For the planar case the solution can be applied to shocks of arbitrary strength, but for cylindrical and spherical geometries it is required that the analysis be restricted to strong shocks. The exact solutions are used to perform a variety of quantitative code verification studies of the Los Alamos National Laboratory Lagrangian hydrocode FLAG.


title: Muon Creation in Supernova Matter Facilitates Neutrino-driven Explosions

abstract: Muons can be created in nascent neutron stars (NSs) due to the high electron chemical potentials and the high temperatures. Because of their relatively lower abundance compared to electrons, their role has so far been ignored in numerical simulations of stellar core collapse and NS formation. However, the appearance of muons softens the NS equation of state, triggers faster NS contraction and thus leads to higher luminosities and mean energies of the emitted neutrinos. This strengthens the postshock heating by neutrinos and can facilitate explosions by the neutrino-driven mechanism.


title: A study of sliding motion of a solid body on a rough surface with asymmetric friction

abstract: Recent studies show interest in materials with asymmetric friction forces. We investigate terminal motion of a solid body with circular contact area. We assume that friction forces are asymmetric orthotropic. Two cases of pressure distribution are analyzed: Hertz and Boussinesq laws. Equations for friction force and moment are formulated and solved for these cases. Numer- ical results show significant impact of the asymmetry of friction on the motion. Our results can be used for more accurate prediction of contact behavior of bodies made from new materials with asymmetric surface textures.


title: Stochastic Gradient MCMC Methods for Hidden Markov Models

abstract: Stochastic gradient MCMC (SG-MCMC) algorithms have proven useful in scaling Bayesian inference to large datasets under an assumption of i.i.d data. We instead develop an SG-MCMC algorithm to learn the parameters of hidden Markov models (HMMs) for time-dependent data. There are two challenges to applying SG-MCMC in this setting: The latent discrete states, and needing to break dependencies when considering minibatches. We consider a marginal likelihood representation of the HMM and propose an algorithm that harnesses the inherent memory decay of the process. We demonstrate the effectiveness of our algorithm on synthetic experiments and an ion channel recording data, with runtimes significantly outperforming batch MCMC.


title: Subjects classification from high-dimensional and small-sample size datasets using a strategy based on Clustering Variables around Latent Components (CLV) method

abstract: High-dimensional complex systems can be studied through multivariate analysis, as Principal Component Analysis, however large samples of observations frequently are needed for it. Here it is examined a method for small samples based on clustering variables around latent variables (CLV) to subject classification in two presumed groups. For it, a predictive model was developed to generate datasets with two groups of cases whose variables show randomness features (up to 30% of variables manifest difference between groups, and up to 7% of those are correlated between them). The method recovered the information of the latent factors to classify the subjects with 80 to 95% of agreement, with positive relationship between the classifier precision and the rate [number of variables / number of subjects].


title: A distributed algorithm for average aggregative games with coupling constraints

abstract: We consider the framework of average aggregative games, where the cost function of each agent depends on his own strategy and on the average population strategy. We focus on the case in which the agents are coupled not only via their cost functions, but also via constraints coupling their strategies. We propose a distributed algorithm that achieves an almost Nash equilibrium by requiring only local communications of the agents, as specified by a sparse communication network. The proof of convergence of the algorithm relies on the auxiliary class of network aggregative games and exploits a novel result of parametric convergence of variational inequalities, which is applicable beyond the context of games. We apply our theoretical findings to a multi-market Cournot game with transportation costs and maximum market capacity.


title: Information Potential Auto-Encoders

abstract: In this paper, we suggest a framework to make use of mutual information as a regularization criterion to train Auto-Encoders (AEs). In the proposed framework, AEs are regularized by minimization of the mutual information between input and encoding variables of AEs during the training phase. In order to estimate the entropy of the encoding variables and the mutual information, we propose a non-parametric method. We also give an information theoretic view of Variational AEs (VAEs), which suggests that VAEs can be considered as parametric methods that estimate entropy. Experimental results show that the proposed non-parametric models have more degree of freedom in terms of representation learning of features drawn from complex distributions such as Mixture of Gaussians, compared to methods which estimate entropy using parametric approaches, such as Variational AEs.


title: Minimal Geometry for Valley Filtering in Graphene

abstract: The possibility to effect valley splitting of an electronic current in graphene represents the essential component in the new field of valleytronics in such two-dimensional materials. Based on a symmetry analysis of the scattering matrix, we show that if the spatial distribution of multiple potential scatterers breaks mirror symmetry about the axis of incoming electrons, then a splitting of the current between two valleys is observed. This leads to the appearance of the valley Hall effect. We illustrate the effect of mirror symmetry breaking in a minimal system of two symmetric impurities, demonstrating the splitting between valleys via the differential cross sections and non-vanishing skew parameter. We further discuss the role that these effects may play in transport experiments.


title: Approximating Gains from Trade in Two-sided Markets via Simple Mechanisms

abstract: We design simple mechanisms to approximate the Gains from Trade (GFT) in two-sided markets with multiple unit-supply sellers and multiple unit-demand buyers. A classical impossibility result by Myerson and Satterthwaite showed that even with only one seller and one buyer, no Individually Rational (IR), Bayesian Incentive Compatible (BIC) and Budget-Balanced (BB) mechanism can achieve full GFT (trade whenever buyer's value is higher than the seller's cost). On the other hand, they proposed the "second-best" mechanism that maximizes the GFT subject to IR, BIC and BB constraints, which is unfortunately rather complex for even the single-seller single-buyer case. Our mechanism is simple, IR, BIC and BB, and achieves $\frac{1}{2}$ of the optimal GFT among all IR, BIC and BB mechanisms. Our result holds for arbitrary distributions of the buyers' and sellers' values and can accommodate any downward-closed feasibility constraints over the allocations. The analysis of our mechanism is facilitated by extending the Cai-Weinberg-Devanur duality framework to two-sided markets.


title: Proximal Backpropagation

abstract: We propose proximal backpropagation (ProxProp) as a novel algorithm that takes implicit instead of explicit gradient steps to update the network parameters during neural network training. Our algorithm is motivated by the step size limitation of explicit gradient descent, which poses an impediment for optimization. ProxProp is developed from a general point of view on the backpropagation algorithm, currently the most common technique to train neural networks via stochastic gradient descent and variants thereof. Specifically, we show that backpropagation of a prediction error is equivalent to sequential gradient descent steps on a quadratic penalty energy, which comprises the network activations as variables of the optimization. We further analyze theoretical properties of ProxProp and in particular prove that the algorithm yields a descent direction in parameter space and can therefore be combined with a wide variety of convergent algorithms. Finally, we devise an efficient numerical implementation that integrates well with popular deep learning frameworks. We conclude by demonstrating promising numerical results and show that ProxProp can be effectively combined with common first order optimizers such as Adam.


title: On the origin of generalized uncertainty principle from compactified M5-brane

abstract: In this paper, we demonstrate that compactification in M-theory can lead to a deformation of field theory consistent with the generalized uncertainty principle (GUP).We observe that the matter fields in the M3-brane action contain higher derivative terms. We demonstrate that such terms can also be constructed from a reformulation of the field theory by the GUP. In fact, we will construct the Heisenberg algebra consistent with this deformation, and explicitly demonstrate it to be the Heisenberg algebra obtained from the GUP. Thus, we use compactification in M-theory to motivate for the existence of the GUP.


title: Stochastic Modeling and Simulation of Viral Evolution

abstract: RNA viruses comprise vast populations of closely related, but highly genetically diverse, entities known as quasispecies. Understanding the mechanisms by which this extreme diversity is generated and maintained is fundamental when approaching viral persistence and pathobiology in infected hosts. In this paper, we access quasispecies theory through a mathematical model based on the theory of multitype branching processes, to better understand the roles of mechanisms resulting in viral diversity, persistence and extinction. We accomplish this understanding by a combination of computational simulations and the theoretical analysis of the model. In order to perform the simulations, we have implemented the mathematical model into a computational platform capable of running simulations and presenting the results in a graphical format in real time. Among other things, we show that the establishment of virus populations may display four distinct regimes from its introduction into new hosts until achieving equilibrium or undergoing extinction. Also, we were able to simulate different fitness distributions representing distinct environments within a host which could either be favorable or hostile to the viral success. We addressed the most used mechanisms for explaining the extinction of RNA virus populations called lethal mutagenesis and mutational meltdown. We were able to demonstrate a correspondence between these two mechanisms implying the existence of a unifying principle leading to the extinction of RNA viruses.


title: Pseudo-deterministic Proofs

abstract: We introduce pseudo-deterministic interactive proofs (psdAM): interactive proof systems for search problems where the verifier is guaranteed with high probability to output the same output on different executions. As in the case with classical interactive proofs, the verifier is a probabilistic polynomial time algorithm interacting with an untrusted powerful prover. We view pseudo-deterministic interactive proofs as an extension of the study of pseudo-deterministic randomized polynomial time algorithms: the goal of the latter is to find canonical solutions to search problems whereas the goal of the former is to prove that a solution to a search problem is canonical to a probabilistic polynomial time verifier. Alternatively, one may think of the powerful prover as aiding the probabilistic polynomial time verifier to find canonical solutions to search problems, with high probability over the randomness of the verifier. The challenge is that pseudo-determinism should hold not only with respect to the randomness, but also with respect to the prover: a malicious prover should not be able to cause the verifier to output a solution other than the unique canonical one.


title: Hydrogen Clouds from Comets 266P Christensen and P2008 Y2 (Gibbs) are Candidates for the Source of the 1977 WOW!Signal

abstract: On 1977 August 15, the Ohio State University Radio Observatory detected a strong narrowband signal northwest of the globular star cluster M55 in the constellation Sagittarius. The frequency of the signal, which closely matched the hydrogen line (1420.40575177 MHz), peaked at approximately 2316 EDT. Since then, several investigations into the Wow signal have ruled out the source as terrestrial in origin or other objects such as satellites, planets and asteroids. From 1977 July 27 to 1977 August 15, comets 266P Christensen and P2008 Y2 (Gibbs) were transiting in the neighborhood of the Chi Sagittarii star group. Because the frequency for the Wow! signal fell close to the hydrogen line, and the hydrogen clouds of 266P/Christensen and P/2008 Y2 (Gibbs) were in the proximity of the right ascension and declination values of the Wow signal, the comet(s) and or their hydrogen clouds are strong candidates for the source of the 1977 Wow signal.


title: Bayesian analysis of accumulated damage models in lumber reliability

abstract: Wood products that are subjected to sustained stress over a period of long duration may weaken, and this effect must be considered in models for the long-term reliability of lumber. The damage accumulation approach has been widely used for this purpose to set engineering standards. In this article, we revisit an accumulated damage model and propose a Bayesian framework for analysis. For parameter estimation and uncertainty quantification, we adopt approximation Bayesian computation (ABC) techniques to handle the complexities of the model. We demonstrate the effectiveness of our approach using both simulated and real data, and apply our fitted model to analyze long-term lumber reliability under a stochastic live loading scenario.


title: A characterization of round spheres in space forms

abstract: Let $\mathbb Q^{n+1}_c$ be the complete simply-connected $(n+1)$-dimensional space form of curvature $c$. In this paper we obtain a new characterization of geodesic spheres in $\mathbb Q^{n+1}_c$ in terms of the higher order mean curvatures. In particular, we prove that the geodesic sphere is the only complete bounded immersed hypersurface in $\mathbb Q^{n+1}_c,\;c\leq 0,$ with constant mean curvature and constant scalar curvature. The proof relies on the well known Omori-Yau maximum principle, a formula of Walter for the Laplacian of the $r$-th mean curvature of a hypersurface in a space form, and a classical inequality of G\aa rding for hyperbolic polynomials.


title: Gauging Metallicity of Diffuse Gas Under An Uncertain Ionizing Radiation Field

abstract: Gas metallicity is a key quantity used to determine the physical conditions of gaseous clouds in a wide range of astronomical environments, including interstellar and intergalactic space. In particular, considerable effort in circumgalactic medium (CGM) studies focuses on metallicity measurements, because gas metallicity serves as a critical discriminator for whether the observed heavy ions in the CGM originate in chemically-enriched outflows or in more chemically-pristine gas accreted from the intergalactic medium. However, because the gas is ionized, a necessary first step in determining CGM metallicity is to constrain the ionization state of the gas which, in addition to gas density, depends on the ultraviolet background radiation field (UVB). While it is generally acknowledged that both the intensity and spectral slope of the UVB are uncertain, the impact of an uncertain spectral slope has not been properly addressed in the literature. This Letter shows that adopting a different spectral slope can result in an order of magnitude difference in the inferred CGM metallicity. Specifically, a harder UVB spectrum leads to a higher estimated gas metallicity for a given set of observed ionic column densities . Therefore, such systematic uncertainties must be folded into the error budget for metallicity estimates of ionized gas. An initial study shows that empirical diagnostics are available for discriminating between hard and soft ionizing spectra. Applying these diagnostics helps reduce the systematic uncertainties in CGM metallicity estimates.


title: Differentially Private Learning of Undirected Graphical Models using Collective Graphical Models

abstract: We investigate the problem of learning discrete, undirected graphical models in a differentially private way. We show that the approach of releasing noisy sufficient statistics using the Laplace mechanism achieves a good trade-off between privacy, utility, and practicality. A naive learning algorithm that uses the noisy sufficient statistics "as is" outperforms general-purpose differentially private learning algorithms. However, it has three limitations: it ignores knowledge about the data generating process, rests on uncertain theoretical foundations, and exhibits certain pathologies. We develop a more principled approach that applies the formalism of collective graphical models to perform inference over the true sufficient statistics within an expectation-maximization framework. We show that this learns better models than competing approaches on both synthetic data and on real human mobility data used as a case study.


title: Infrared emission from kilonovae: the case of the nearby short hard burst GRB 160821B

abstract: We present constraints on Ks-band emission from one of the nearest short hard gamma-ray bursts, GRB 160821B, at z=0.16, at three epochs. We detect a reddened relativistic afterglow from the jetted emission in the first epoch but do not detect any excess kilonova emission in the second two epochs. We compare upper limits obtained with Keck I/MOSFIRE to multi-dimensional radiative transfer models of kilonovae, that employ composition-dependent nuclear heating and LTE opacities of heavy elements. We discuss eight models that combine toroidal dynamical ejecta and two types of wind and one model with dynamical ejecta only. We also discuss simple, empirical scaling laws of predicted emission as a function of ejecta mass and ejecta velocity. Our limits for GRB 160821B constrain the ejecta mass to be lower than 0.03 Msun for velocities greater than 0.1c. At the distance sensitivity range of advanced LIGO, similar ground-based observations would be sufficiently sensitive to the full range of predicted model emission including models with only dynamical ejecta. The color evolution of these models shows that I-K color spans 7--16 mag, which suggests that even relatively shallow infrared searches for kilonovae could be as constraining as optical searches.


title: A Novel Non-Perturbative Lattice Regularization of an Anomaly-Free $1 + 1d$ Chiral $SU(2)$ Gauge Theory

abstract: We present a numerical treatment of a novel non-perturbative lattice regularization of a $1+1d$ $SU(2)$ Chiral Gauge Theory. Our approach follows recent proposals that exploit the newly discovered connection between anomalies and topological (or entangled) states to show how to create a lattice regularization of any anomaly-free chiral gauge theory. In comparison to other methods, our regularization enjoys on-site fermions and gauge action, as well as a physically transparent fermion Hilbert space. We follow the `mirror fermion' approach, in which we first create a lattice regularization of both the chiral theory and its mirror conjugate and then introduce interactions that gap out only the mirror theory. The connection between topological states and anomalies shows that such interactions exist if the chiral theory is free of all quantum anomalies. Instead of numerically intractable fermion-fermion interactions, we couple the mirror theory to a Higgs field driven into a symmetry-preserving, disordered, gapped phase.


title: Parabolic Catalan numbers count flagged Schur functions and their appearances as type A Demazure characters (key polynomials)

abstract: Fix an integer partition lambda that has no more than n parts. Let beta be a weakly increasing n-tuple with entries from {1,..,n}. The flagged Schur function indexed by lambda and beta is a polynomial generating function in x_1, .., x_n for certain semistandard tableaux of shape lambda. Let pi be an n-permutation. The type A Demazure character (key polynomial, Demazure polynomial) indexed by lambda and pi is another such polynomial generating function. Reiner and Shimozono and then Postnikov and Stanley studied coincidences between these two families of polynomials. Here their results are sharpened by the specification of unique representatives for the equivalence classes of indexes for both families of polynomials, extended by the consideration of more general beta, and deepened by proving that the polynomial coincidences also hold at the level of the underlying tableau sets. Let R be the set of lengths of columns in the shape of lambda that are less than n. Ordered set partitions of {1,..,n} with block sizes determined by R, called R-permutations, are used to describe the minimal length representatives for the parabolic quotient of the nth symmetric group specified by the set {1,..,n-1}\R. The notion of 312-avoidance is generalized from n-permutations to these set partitions. The R-parabolic Catalan number is defined to be the number of these. Every flagged Schur function arises as a Demazure polynomial. Those Demazure polynomials are precisely indexed by the R-312-avoiding R-permutations. Hence the number of flagged Schur functions that are distinct as polynomials is shown to be the R-parabolic Catalan number. The projecting and lifting processes that relate the notions of 312-avoidance and of R-312-avoidance are described with maps developed for other purposes.


title: Schatten class Hankel and $\overline{\partial}$-Neumann operators on pseudoconvex domains in $\mathbb{C}^n$

abstract: Let $\Omega$ be a $C^2$-smooth bounded pseudoconvex domain in $\mathbb{C}^n$ for $n\geq 2$ and let $\varphi$ be a holomorphic function on $\Omega$ that is $C^2$-smooth on the closure of $\Omega$. We prove that if $H_{\overline{\varphi}}$ is in Schatten $p$-class for $p\leq 2n$ then $\varphi$ is a constant function. As a corollary, we show that the $\overline{\partial}$-Neumann operator on $\Omega$ is not Hilbert-Schmidt.


title: Spatial Regression and the Bayesian Filter

abstract: Regression for spatially dependent outcomes poses many challenges, for inference and for computation. Non-spatial models and traditional spatial mixed-effects models each have their advantages and disadvantages, making it difficult for practitioners to determine how to carry out a spatial regression analysis. We discuss the data-generating mechanisms implicitly assumed by various popular spatial regression models, and discuss the implications of these assumptions. We propose Bayesian spatial filtering as an approximate middle way between non-spatial models and traditional spatial mixed models. We show by simulation that our Bayesian spatial filtering model has several desirable properties and hence may be a useful addition to a spatial statistician's toolkit.


title: Learning a visuomotor controller for real world robotic grasping using simulated depth images

abstract: We want to build robots that are useful in unstructured real world applications, such as doing work in the household. Grasping in particular is an important skill in this domain, yet it remains a challenge. One of the key hurdles is handling unexpected changes or motion in the objects being grasped and kinematic noise or other errors in the robot. This paper proposes an approach to learning a closed-loop controller for robotic grasping that dynamically guides the gripper to the object. We use a wrist-mounted sensor to acquire depth images in front of the gripper and train a convolutional neural network to learn a distance function to true grasps for grasp configurations over an image. The training sensor data is generated in simulation, a major advantage over previous work that uses real robot experience, which is costly to obtain. Despite being trained in simulation, our approach works well on real noisy sensor images. We compare our controller in simulated and real robot experiments to a strong baseline for grasp pose detection, and find that our approach significantly outperforms the baseline in the presence of kinematic noise, perceptual errors and disturbances of the object during grasping.


title: The reactions $\pi\pi\rightarrow\pi\pi$ and $\gamma\gamma\rightarrow\pi\pi$ in $\chi$PT with an isosinglet scalar resonance

abstract: The lowest-lying resonance in the QCD spectrum is the $0^{++}$ isoscalar $\sigma$ meson, also known as the $f_0(500)$. We augment SU(2) chiral perturbation theory ($\chi$PT) by including the $\sigma$ meson as an additional explicit degree of freedom, as proposed by Soto, Talavera, and Tarr\'us and others. In this effective field theory, denoted $\chi$PT$_S$, the $\sigma$ meson's well-established mass and decay width are not sufficient to properly renormalize its self energy. At $\mathcal{O}(p^4)$ another low-energy constant appears in the dressed $\sigma$-meson propagator; we adjust it so that the isoscalar pion-pion scattering length is also reproduced. We compare the resulting amplitudes for the $\pi\pi\rightarrow\pi\pi$ and $\gamma\gamma\rightarrow\pi\pi$ reactions to data from threshold through the energies at which the $\sigma$-meson resonance affects observables. The leading-order (LO) $\pi \pi$ amplitude reproduces the $\sigma$-meson pole position, the isoscalar $\pi \pi$ scattering lengths and $\pi \pi$ scattering and $\gamma \gamma \rightarrow \pi \pi$ data up to $\sqrt{s} \approx 0.5$ GeV. It also yields a $\gamma\gamma\rightarrow\pi\pi$ amplitude that obeys the Ward identity. The value obtained for the $\pi^0$ polarizability is, however, only slightly larger than that obtained in standard $\chi$PT.


title: Experimental observation of magnetic bobbers for a new concept of magnetic solid-state memory

abstract: The use of chiral skyrmions, which are nanoscale vortex-like spin textures, as movable data bit carriers forms the basis of a recently proposed concept for magnetic solid-state memory. In this concept, skyrmions are considered to be unique localized spin textures, which are used to encode data through the quantization of different distances between identical skyrmions on a guiding nanostripe. However, the conservation of distances between highly mobile and interacting skyrmions is difficult to implement in practice. Here, we report the direct observation of another type of theoretically-predicted localized magnetic state, which is referred to as a chiral bobber (ChB), using quantitative off-axis electron holography. We show that ChBs can coexist together with skyrmions. Our results suggest a novel approach for data encoding, whereby a stream of binary data representing a sequence of ones and zeros can be encoded via a sequence of skyrmions and bobbers. The need to maintain defined distances between data bit carriers is then not required. The proposed concept of data encoding promises to expedite the realization of a new generation of magnetic solid-state memory.


title: Particle Detection Algorithms for Complex Plasmas

abstract: In complex plasmas, the behavior of freely floating micrometer sized particles is studied. The particles can be directly visualized and recorded by digital video cameras. To analyze the dynamics of single particles, reliable algorithms are required to accurately determine their positions to sub-pixel accuracy from the recorded images. Typically, straightforward algorithms are used for this task. Here, we combine the algorithms with common techniques for image processing. We study several algorithms and pre- and post-processing methods, and we investigate the impact of the choice of threshold parameters, including an automatic threshold detection. The results quantitatively show that each algorithm and method has its own advantage, often depending on the problem at hand. This knowledge is applicable not only to complex plasmas, but useful for any kind of comparable image-based particle tracking, e.g. in the field of colloids or granular matter.


title: Scattering Cross Section in a Cylindrical anisotropic layered metamaterial

abstract: To design a uniaxial anisotropic metamaterial a layered cylindrical metamaterial is introduced for TE polarization. Unlike to the previous work, which the layers were in radial direction, here the layers are in azimuthal direction. Scattering efficiency for this metamaterial in different frequency is analyzed with solving Maxwell's wave equation. It is observed that in some frequencies when the effective permittivity of the structure goes to zero the scattering efficiency would be negligible. This result approves the previous predictions. It is also found out that the scattering cancellation depends on the relative permittivity of the environmental medium for the cylinder. The finite element simulations are also confirmed the results.


title: Proper motions of collimated jets from intermediate-mass protostars in the Carina Nebula

abstract: We present proper motion measurements of 37 jets and HH objects in the Carina Nebula measured in two epochs of H$\alpha$ images obtained $\sim 10$ yrs apart with HST/ACS. Transverse velocities in all but one jet are faster than $\gtrsim 25$ km s$^{-1}$, confirming that the jet-like H$\alpha$ features identified by Smith et al. (2010) trace outflowing gas. Proper motions constrain the location of the jet-driving source and provide kinematic confirmation of the intermediate-mass protostars that we identify for 20/37 jets. Jet velocities do not correlate with the estimated protostar mass and embedded driving sources do not have slower jets. Instead, transverse velocities (median $\sim 75$ km s$^{-1}$) are similar to those in jets from low-mass stars. Assuming a constant velocity since launch, we compute jet dynamical ages (median $\sim 10^4$ yr). If continuous emission from inner jets traces the duration of the most recent accretion bursts, then these episodes are sustained longer (median $\sim 700$ yr) than the typical decay time of an FU Orionis outburst. These jets can carry appreciable momentum that may be injected into the surrounding environment. The resulting outflow force, $dP/dt$, lies between that measured in low- and high-mass sources, despite the very different observational tracers used. Smooth scaling of the outflow force argues for a common physical process underlying outflows from protostars of all masses. This latest kinematic result adds to a growing body of evidence that intermediate-mass star formation proceeds like a scaled-up version of the formation of low-mass stars.


title: Aa Unbalanced Feed Design for Wideband Phased Arrays

abstract: In this work, a planar phased array antenna was engineered with ultra-wideband (UWB) performance that covers portions of the L-, S-, and C-bands. The unit cell design contains a driven layer that is directly connected to a single coaxial feed and a parasitic layer located above the driven layer separated by free space. This design does not require a balun or any vias between the antenna ground plane and the driven layer, resulting in a simple antenna stack-up consisting of only planar layers, aside from the feed structure. As a consequence, the complexity, and potentially sensitivity to mechanical tolerances, is reduced. The simulated results of this unbalanced phased array are discussed and experimentally validated.


title: On the radius of spatial analyticity for cubic nonlinear Schr\"{o}dinger equation

abstract: It is shown that the uniform radius of spatial analyticity $\sigma(t)$ of solutions at time $t$ to the 1d, 2d and 3d cubic nonlinear Schr\"{o}dinger equations cannot decay faster than $1/|t|$ as $|t| \to \infty$, given initial data that is analytic with fixed radius $\sigma_0$.


title: On Counting Triangles through Edge Sampling in Large Dynamic Graphs

abstract: Traditional frameworks for dynamic graphs have relied on processing only the stream of edges added into or deleted from an evolving graph, but not any additional related information such as the degrees or neighbor lists of nodes incident to the edges. In this paper, we propose a new edge sampling framework for big-graph analytics in dynamic graphs which enhances the traditional model by enabling the use of additional related information. To demonstrate the advantages of this framework, we present a new sampling algorithm, called Edge Sample and Discard (ESD). It generates an unbiased estimate of the total number of triangles, which can be continuously updated in response to both edge additions and deletions. We provide a comparative analysis of the performance of ESD against two current state-of-the-art algorithms in terms of accuracy and complexity. The results of the experiments performed on real graphs show that, with the help of the neighborhood information of the sampled edges, the accuracy achieved by our algorithm is substantially better. We also characterize the impact of properties of the graph on the performance of our algorithm by testing on several Barabasi-Albert graphs.


title: Engineering the breaking of time-reversal symmetry in gate-tunable hybrid ferromagnet/topological insulator heterostructures

abstract: Studying the influence of breaking time-reversal symmetry on topological insulator surface states is an important problem of current interest in condensed matter physics and could provide a route toward proof-of-concept spintronic devices that exploit spin-textured surface states. Here, we develop a new model system for studying the effect of breaking time-reversal symmetry: a hybrid heterostructure wherein a ferromagnetic semiconductor Ga1-xMnxAs, with an out-of-plane component of magnetization, is cleanly interfaced with a three-dimensional topological insulator (Bi,Sb)2(Te,Se)3 by molecular beam epitaxy. Lateral electrical transport in this bilayer is dominated by conduction through the topological insulator whose conductivity is a few orders of magnitude higher than that of the highly resistive ferromagnetic semiconductor with a low Mn concentration. Electrical transport measurements of a top-gated heterostructure device reveal a crossover from weak anti-localization (negative magneto-conductance) to weak localization (positive magneto-conductance) as the temperature is lowered or as the chemical potential approaches the Dirac point. This is accompanied by a systematic emergence of an anomalous Hall effect. These results are interpreted in terms of the opening of a gap at the Dirac point as a result of the exchange coupling between the topological insulator surface state and the ferromagnetic ordering in the Ga1-xMnxAs layer. Our study shows that this hybrid system is well suited to explore topological quantum phenomena and to realize proof-of-concept demonstrations of topological spintronic devices at cryogenic temperatures.


title: Polariton chimeras: Bose-Einstein condensates with intrinsic chaoticity and spontaneous long-range ordering

abstract: The system of cavity polaritons driven by a plane electromagnetic wave is found to undergo the spontaneous breaking of spatial symmetry, which results in a lifted phase locking with respect to the driving field and, consequently, in the possibility of internal ordering. In particular, periodic spin and intensity patterns arise in polariton wires; they exhibit strong long-range order and can serve as media for signal transmission. Such patterns have the properties of dynamical chimeras: they are formed spontaneously in perfectly homogeneous media and can be partially chaotic. The reported new mechanism of chimera formation requires neither time-delayed feedback loops nor non-local interactions.