title: Complexity Equals Action

abstract: We conjecture that the quantum complexity of a holographic state is dual to the action of a certain spacetime region that we call a Wheeler-DeWitt patch. We illustrate and test the conjecture in the context of neutral, charged, and rotating black holes in AdS, as well as black holes perturbed with static shells and with shock waves. This conjecture evolved from a previous conjecture that complexity is dual to spatial volume, but appears to be a major improvement over the original. In light of our results, we discuss the hypothesis that black holes are the fastest computers in nature.


title: A Constrained-Gradient Method to Control Divergence Errors in Numerical MHD

abstract: In numerical magnetohydrodynamics (MHD), a major challenge is maintaining zero magnetic field-divergence (div-B). Constrained transport (CT) schemes can achieve this at high accuracy, but have generally been restricted to very specific methods. For more general (meshless, moving-mesh, or ALE) methods, 'divergence-cleaning' schemes reduce the div-B errors, however they can still be significant, especially at discontinuities, and can lead to systematic deviations from correct solutions which converge away very slowly. Here we propose a new constrained gradient (CG) scheme which augments these with a hybrid projection step, and can be applied to any numerical scheme with a reconstruction. This iteratively approximates the least-squares minimizing, globally divergence-free reconstruction of the fluid. We emphasize that, unlike 'locally divergence free' methods, this actually minimizes the numerically unstable div-B terms, without affecting the convergence order of the method. We implement this in the mesh-free code GIZMO and compare a wide range of test problems. Compared to state-of-the-art cleaning schemes, our CG method reduces the maximum div-B errors in each problem by 1-3 orders of magnitude (2-5 dex below the typical errors if no div-B cleaning is used). By preventing large div-B even at unresolved discontinuities, the method eliminates systematic errors at jumps. In every problem, the accuracy of our CG results is comparable to CT methods. The cost is modest, ~30% of the hydro algorithm, and the CG correction can be easily implemented in a wide range of different numerical MHD methods. While for many problems, we find Dedner-type cleaning schemes are sufficient for good results, we identify a wide range of problems where using only the simplest Powell or '8-wave' cleaning can produce systematic, order-of-magnitude errors.

categories: astro-ph.IM astro-ph.CO astro-ph.GA astro-ph.SR physics.flu-dyn


title: Outer-disk reddening and gas-phase metallicities: The CALIFA connection

abstract: We study, for the first time in a statistically significant and well-defined sample, the relation between the outer-disk ionized-gas metallicity gradients and the presence of breaks in the surface brightness profiles of disk galaxies. SDSS g'- and r'-band surface brightness, (g'- r') color, and ionized-gas oxygen abundance profiles for 324 galaxies within the CALIFA survey are used for this purpose. We perform a detailed light-profile classification finding that 84% of our disks show down- or up-bending profiles (Type II and Type III, respectively) while the remaining 16% are well fitted by one single exponential (Type I). The analysis of the color gradients at both sides of this break shows a U-shaped profile for most Type II galaxies with an average minimum (g'- r') color of ~0.5 mag and a ionized-gas metallicity flattening associated to it only in the case of low-mass galaxies. More massive systems show a rather uniform negative metallicity gradient. The correlation between metallicity flattening and stellar mass results in p-values as low as 0.01. Independently of the mechanism having shaped the outer light profiles of these galaxies, stellar migration or a previous episode of star formation in a shrinking star-forming disk, it is clear that the imprint in their ionized-gas metallicity was different for low- and high-mass Type II galaxies. In the case of Type III disks, a positive correlation between the change in color and abundance gradient is found (the null hypothesis is ruled out with a p-value of 0.02), with the outer disks of Type III galaxies with masses $\leq$10$^{10}$ M$_{\odot}$ showing a weak color reddening or even a bluing. This is interpreted as primarily due to a mass down-sizing effect on the population of Type III galaxies having recently experienced an enhanced inside-out growth.


title: Van der Waals heterojunction devices based on organohalide perovskites and two-dimensional materials

abstract: The recently emerged organohalide perovskites (e.g., CH3NH3PbI3) have drawn intense attention for high efficiency solar cells. However, with a considerable solubility in many solvents, these perovskites are not typically compatible with conventional lithography processes for more complicated device fabrications that are important for both fundamental studies and technological applications. Here we report the creation of novel heterojunction devices based on perovskites and two-dimensional (2D) crystals by taking advantage of the layered characteristic of lead iodide (PbI2) and vapor phase intercalation. We show a graphene/perovskite/graphene vertical stack can deliver a highest photoresponsivity of ~950 A/W and photoconductive gain of ~2200, and a graphene/WSe2/perovskite/graphene heterojunction can display a high on/off ratio (~10^6) transistor behavior with distinct gate-tunable diode characteristics and open-circuit voltages. Such unique perovskite-2D heterostructures have significant potential for future optoelectronic research and can enable broad possibilities with compositional tunability of organohalide perovskites and the versatility offered by diverse 2D materials.


title: Archival Legacy Investigations of Circumstellar Environments (ALICE): Statistical assessment of point source detections

abstract: The ALICE program, for Archival Legacy Investigation of Circumstellar Environment, is currently conducting a virtual survey of about 400 stars, by re-analyzing the HST-NICMOS coronagraphic archive with advanced post-processing techniques. We present here the strategy that we adopted to identify detections and potential candidates for follow-up observations, and we give a preliminary overview of our detections. We present a statistical analysis conducted to evaluate the confidence level on these detection and the completeness of our candidate search.


title: Thermodynamic and Kinetic Anisotropies in Octane Thin Films

abstract: Confinement breaks the translational symmetry of materials. Such symmetry breaking can be used to obtain configurations that are not otherwise accessible in the bulk. Here, we explore the effect of substrate-liquid interactions on the induced thermodynamic and kinetic anisotropies. We consider n-octane nanofilms that are in contact with substrates with varying degrees of attraction. Complete freezing of octane nanofilms is observed at low temperatures, while at intermediate temperatures, a frozen monolayer emerges at both interfaces. By carefully inspecting the profiles of translational and orientational relaxation times, we confirm that the translational and orientational degrees of freedom are decoupled at these frozen monolayers. At sufficiently high temperatures, however, free interfaces and solid-liquid interfaces close to loose substrates undergo pre-freezing, characterized by mild peaks in several thermodynamic quantities. Two distinct dynamic regimes are observed. The dynamics is accelerated in the vicinity of loose substrates, while sticky substrates decelerate dynamics, sometimes by as much as two orders of magnitude. These two distinct dynamical regimes have been previously by us [JCP 141: 024506, 2014] for a model atomic glass-forming liquid. We also confirm the existence of two correlations proposed in the above-mentioned work in solid-liquid subsurface regions of octane films, i.e., a correlation between density and normal stress, and between atomic translational relaxation time and lateral stress. Finally, we inspect the ability of different regions of a film to explore the potential energy landscape, and observe no noticeable difference between the free surface and the bulk. This is unlike the films of model atomic glass formers that tend to sample their respective landscape more efficiently at free surfaces.


title: Sliding droplets of Xanthan solutions: a joint experimental and numerical study

abstract: We have investigated the sliding of droplets made of solutions of Xanthan, a stiff rodlike polysaccharide exhibiting a non-newtonian behavior, notably characterized by a shear-rate dependence of the viscosity. These experimental results are quantitatively compared with those of newtonian fluids (water). The impact of the non-newtonian behavior on the sliding process was shown through the relation between the average dimensionless velocity (i.e. the Capillary number) and the dimensionless volume forces (i.e. the Bond number). To this aim, it is needed to define operative strategies to compute the Capillary number for the shear thinning fluids and compare with the corresponding newtonian case. Results from experiments were complemented with lattice Boltzmann numerical simulations of sliding droplets, aimed to disentangle the influence that non-newtonian flow properties have on the sliding.


title: Determinantal representations of W-weighted Drazin inverse solutions of some quaternion matrix equations

abstract: By using determinantal representations of the W-weighted Drazin inverse previously obtained by the author within the framework of the theory of the column-row determinants, we get explicit formulas for determinantal representations of the W-weighted Drazin inverse solutions (analogs of Cramer's rule) of the quaternion matrix equations $ {\bf W}{\bf A}{\bf W}{\bf X}={\bf D}$, $ {\bf X}{\bf W}{\bf A}{\bf W}={\bf D} $, and ${\bf W}_{1}{\bf A}{\bf W}_{1}{\bf X}{\bf W}_{2}{\bf B}{\bf W}_{2}={\bf D} $.


title: Berline-Vergne valuation and generalized permutohedra

abstract: Generalizing a conjecture by De Loera et al., we conjecture that integral generalized permutohedra all have positive Ehrhart coefficients. Berline and Vergne construct a valuation that assigns values to faces of polytopes, which provides a way to write Ehrhart coefficients of a polytope as positive sums of these values. Based on available results, we pose a stronger conjecture: Berline-Vergne's valuation is always positive on permutohedra, which implies our first conjecture. This article proves that our strong conjecture on Berline-Vergne's valuation is true for dimension up to 6, and is true if we restrict to faces of codimension up to 3. In addition to investigating the positivity conjectures, we study the Berline-Vergne's valuation, and show that it is the unique construction for McMullen's formula used to describe number of lattice points in permutohedra under certain symmetry constraints. We also give an equivalent statement to the strong conjecture in terms of mixed valuations.


title: Strong subadditivity and holography

abstract: We study in detail the relationship between strong subadditivity for a boundary field theory and energy conditions for its bulk dual in 2+1 dimensions. We provide a discussion of known facts and new results organized from the simplest case of a static system with collinear intervals to a time dependent one in a generic configuration, with particular focus on the holographic geometric description.


title: Systematic Renormalization of the Effective Theory of Large Scale Structure

abstract: A perturbative description of Large Scale Structure is a cornerstone of our understanding of the observed distribution of matter in the universe. Renormalization is an essential and defining step to make this description physical and predictive. Here we introduce a systematic renormalization procedure, which neatly associates counterterms to the UV-sensitive diagrams order by order, as it is commonly done in quantum field theory. As a concrete example, we renormalize the one-loop power spectrum and bispectrum of both density and velocity. In addition, we present a series of results that are valid to all orders in perturbation theory. First, we show that while systematic renormalization requires temporally non-local counterterms, in practice one can use an equivalent basis made of local operators. We give an explicit prescription to generate all counterterms allowed by the symmetries. Second, we present a formal proof of the well-known general argument that the contribution of short distance perturbations to large scale density contrast $\delta$ and momentum density $\mathbf\pi(\mathbf k)$ scale as $k^2$ and $k$, respectively. Third, we demonstrate that the common practice of introducing counterterms only in the Euler equation when one is interested in correlators of $ \delta$ is indeed valid to all orders.


title: Development of a coarse-grained water forcefield via multistate iterative Boltzmann inversion

abstract: A coarse-grained water model is developed using multistate iterative Boltzmann inversion. Following previous work, the k-means algorithm is used to dynamically map multiple water molecules to a single coarse-grained bead, allowing the use of structure-based coarse-graining methods. The model is derived to match the bulk and interfacial properties of liquid water and improves upon previous work that used single state iterative Boltzmann inversion. The model accurately reproduces the density and structural correlations of water at 305 K and 1.0 atm, stability of a liquid droplet at 305 K, and shows little tendency to crystallize at physiological conditions. This work also illustrates several advantages of using multistate iterative Boltzmann inversion for deriving generally applicable coarse-grained forcefields.


title: Synchronization of Distant Optical Clocks at the Femtosecond Level

abstract: The use of optical clocks/oscillators in future ultra-precise navigation, gravitational sensing, coherent arrays, and relativity experiments will require time comparison and synchronization over terrestrial or satellite free-space links. Here we demonstrate full unambiguous synchronization of two optical timescales across a free-space link. The time deviation between synchronized timescales is below 1 fs over durations from 0.1 s to 6500 s, despite atmospheric turbulence and kilometer-scale path length variations. Over several days, the time wander is 40 fs peak-to-peak. Our approach relies on the two-way reciprocity of a single-spatial-mode optical link, valid to below 225 attoseconds across a turbulent 4-km path. This femtosecond level of time-frequency transfer should enable optical networks using state-of-the-art optical clocks/oscillators.


title: Isotropic Negative Thermal Expansion Metamaterials

abstract: Negative thermal expansion materials are important and desirable in science and engineering applications. However, natural materials with isotropic negative thermal expansion are rare and usually unsatisfied in performance. Here, we propose a novel method to achieve negative thermal expansion via a metamaterial approach. The metamaterial is constructed with unit cells that combine bi-material strips and anti-chiral structures. Both experimental and simulation results display isotropic negative thermal expansion properties. The coefficient of negative thermal expansion of the metamaterials is demonstrated to be dependent on the difference between the thermal expansion coefficients of two component materials in the bi-material strips, as well as on the circular node radius and the ligament length in the anti-chiral structures. The measured value of the linear negative thermal expansion coefficient reaches -68.1X10-6 1/K in an operating temperature range from 303.15 K to 773.15 K, which is among the largest achieved in experiments to date. Our findings provide a novel and practical approach to obtaining materials with tunable isotropic negative thermal expansion on any scale.


title: Effects of Sliding Speed on the Intensity of Triboluminescence in Slide contact: Experimental Measurements and Theoretical Analyses

abstract: Triboluminescence (TL) is the emission of light produced by rubbing or striking two materials together. Here, the light emission has been observed from the sliding contact between two disks under dry condition. The effects of the sliding speed on the intensity of TL have been experimentally investigated. The results show that the intensity of the emission light increases significantly with the sliding speed. A theoretical model is also proposed and an analytical expression is deduced for the intensity of TL in the slide contact. The theoretical prediction is found consistent with the experimental results. The present work may be helpful to the understanding of the mechanism of light emission when friction.


title: Thin-disk laser pump schemes for large number of passes and moderate pump source quality

abstract: Novel thin-disk laser pump layouts are proposed yielding an increased number of passes for a given pump module size and pump source quality. These novel layouts result from a general scheme which bases on merging two simpler pump optics arrangements. Some peculiar examples can be realized by adapting standard commercially available pump optics simply by intro ducing an additional mirror-pair. More pump passes yield better efficiency, opening the way for usage of active materials with low absorption. In a standard multi-pass pump design, scaling of the number of beam passes brings ab out an increase of the overall size of the optical arrangement or an increase of the pump source quality requirements. Such increases are minimized in our scheme, making them eligible for industrial applications


title: Evasion and Hardening of Tree Ensemble Classifiers

abstract: Classifier evasion consists in finding for a given instance $x$ the nearest instance $x'$ such that the classifier predictions of $x$ and $x'$ are different. We present two novel algorithms for systematically computing evasions for tree ensembles such as boosted trees and random forests. Our first algorithm uses a Mixed Integer Linear Program solver and finds the optimal evading instance under an expressive set of constraints. Our second algorithm trades off optimality for speed by using symbolic prediction, a novel algorithm for fast finite differences on tree ensembles. On a digit recognition task, we demonstrate that both gradient boosted trees and random forests are extremely susceptible to evasions. Finally, we harden a boosted tree model without loss of predictive accuracy by augmenting the training set of each boosting round with evading instances, a technique we call adversarial boosting.


title: The role of the Havriliak-Negami relaxation in the description of local structure of Kohlrausch's function in the frequency domain. Part I

abstract: An improved approximation via Havriliak-Negami (HN) functions to the Fourier Transform (FT) of certain Weibull distributions, -\psi_{\beta}, (the time derivative of the Kohlrausch-Williams-Watts function), is given for a large interval of frequencies: \omega/2\pi\in[0,10^{12}] if 0<\beta\leq1 and \omega/2\pi\in[0,10^{7}] if 1<\beta\leq2. The model is free from the usual numerical distortions, or restrictions associated to sampling step and finite size, present in similar adjustments to complex relaxation functions. Further indicates that the identification of (FT) Weibull data with a double HN approximant, \psi_{\beta}\simeq\mathcal{A}p_{2}HN, is quite exact locally even though the parameters involved should vary adiabatically with the frequency, i.e. \{\alpha_{1,2},\gamma_{1,2},\tau_{1,2},\lambda\}(\omega). This fact is the base for the high sensibility of the parameter models, as functions of \beta, to the available data and sampling associated errors. Presumably is also the root of the different proposals for asymptotic traits of \alpha\cdot\gamma(\beta)_{1,2}|_{\omega\rightarrow\infty} functions found in scientific literature.


title: Light Leptonic New Physics at the Precision Frontier

abstract: Precision probes of new physics are often interpreted through their indirect sensitivity to short-distance scales. In this proceedings contribution, we focus on the question of which precision observables, at current sensitivity levels, allow for an interpretation via either short-distance new physics or consistent models of long-distance new physics, weakly coupled to the Standard Model. The electroweak scale is chosen to set the dividing line between these scenarios. In particular, we find that inverse see-saw models of neutrino mass allow for light new physics interpretations of most precision leptonic observables, such as lepton universality, lepton flavor violation, but not for the electron EDM.


title: Hadronic resonance production and interaction in partonic and hadronic matter in EPOS3 with and without the hadronic afterburner UrQMD

abstract: We study the production of hadronic resonances and their interaction in the partonic and hadronic medium using the EPOS3 model, which employs the UrQMD model for the description of the hadronic phase. We investigate the centrality dependence of the yields and momentum distributions for various resonances (rho(770),K(892),phi(1020),Delta(1232),Sigma(1385),Lambda(1520),Xi(1530) and their antiparticles) in Pb-Pb collisions at sNN= 2.76 TeV. The predictions for K(892) and phi(1020) will be compared with the experimental data from the ALICE collaboration. The observed signal suppression of the K(892) with increasing centrality will be discussed with respect to the resonance interaction in the hadronic medium. The mean transverse momentum and other particle ratios such as phi(1020)/p and Omega/phi(1020) will be discussed with respect to additional contributions from the hadronic medium interactions.


title: Asymmetry energy of nuclear matter: Temperature and density dependence, and validity of semi-empirical formula

abstract: In this work, we have done a completely microscopic calculation using a many-body variational method based on the cluster expansion of energy to compute the asymmetry energy of nuclear matter. In our calculations, we have employed the $AV_{18}$ nuclear potential. We have also investigated the temperature and density dependence of asymmetry energy. Our results show that the asymmetry energy of nuclear matter depends on both density and temperature. We have also studied the effects of different terms in the asymmetry energy of nuclear matter. These investigations indicate that at different densities and temperatures, the contribution of parabolic term is very substantial with respect to the other terms. Therefore, we can conclude that the parabolic approximation is a relatively good estimation, and our calculated binding energy of asymmetric nuclear matter is in a relatively good agreement with that of semi-empirical mass formula. However, for the accurate calculations, it is better to consider the effects of other terms.


title: Quantum Look at two Common Logics: the Logic of Primitive Thinking and the Logic of Everyday Human Reasoning

abstract: Based on ideas of quantum theory of open systems and psychological dual system theory we propose two novel versions of Non-Boolean logic. The first version can be interpreted in our opinion as simplified description of primitive (mythological) thinking and the second one as the toy model of everyday human reasoning in which aside from logical deduction, heuristic elements and beliefs also play the considerable role. Several arguments in favor of the interpretations proposed are adduced and discussed in the paper as well.


title: Fibonacci and Lucas Numbers Associated with Brocard-Ramanujan Equation

abstract: We explicitly solve the diophantine equations of the form $$ A_{n_1}A_{n_2}\cdots A_{n_k}\pm 1 = B_m^2 $$ where $(A_n)_{n\geq 0}$ and $(B_m)_{m\geq 0}$ are either the Fibonacci sequence or Lucas sequence. This extends the result of D. Marques (2011) and L. Szalay (2012) concerning a variant of Brocard-Ramanujan equation. This is a manuscript of the article published in Communications of the Korean Mathematical Society, 32(3) (2017), pp. 511-522


title: Cultural Barriers to Software Productivity Practices at Los Alamos

abstract: In recent years, code projects in the nuclear weapons program at Los Alamos National Laboratory (LANL) have given increased attention to modern software productivity practices. We found that some of the biggest barriers to adoption of new practices were not technical but cultural. This paper describes several of the cultural differences between the physics and computer science communities at LANL.


title: Bayesian sequential parameter estimation with a Laplace type approximation

abstract: A method for sequential inference of the fixed parameters of a dynamic latent Gaussian models is proposed and evaluated that is based on the iterated Laplace approximation. The method provides a useful trade-off between computational performance and the accuracy of the approximation to the true posterior distribution. Approximation corrections are shown to improve the accuracy of the approximation in simulation studies. A population-based approach is also shown to provide a more robust inference method.


title: Simulations of Solar Jets Confined by Coronal Loops

abstract: Coronal jets are collimated, dynamic events that occur over a broad range of spatial scales in the solar corona. In the open magnetic field of coronal holes, jets form quasi-radial spires that can extend far out into the heliosphere, while in closed-field regions the jet outflows are confined to the corona. We explore the application of the embedded-bipole model to jets occurring in closed coronal loops. In this model, magnetic free energy is injected slowly by footpoint motions that introduce twist within the closed dome of the jet source region, and is released rapidly by the onset of an ideal kink-like instability. Two length scales characterize the system: the width (N) of the jet source region and the footpoint separation (L) of the coronal loop that envelops the jet source. We find that the jet characteristics are highly sensitive to the ratio L/N, in both the conditions for initiation and the subsequent dynamics. The longest-lasting and most energetic jets occur along long coronal loops with large L/N ratios, and share many features of open-field jets, while smaller L/N ratios produce shorter-duration, less energetic jets that are affected by reflections from the far-loop footpoint. We quantify the transition between these behaviours and show that our model replicates key qualitative and quantitative aspects of both quiet-Sun and active-region loop jets. We also find that the reconnection between the closed dome and surrounding coronal loop is very extensive: the cumulative reconnected flux at least matches the total flux beneath the dome for small L/N, and is more than double that value for large L/N.


title: Enabling Arbitrary Wavelength Optical Frequency Combs on Chip

abstract: A necessary condition for generation of bright soliton Kerr frequency combs in microresonators is to achieve anomalous group velocity dispersion (GVD) for the resonator modes. This condition is hard to implement in visible as well as ultraviolet since the majority of optical materials are characterized with large normal GVD in these wavelength regions. We overcome this challenge by borrowing ideas from strongly dispersive coupled systems in solid state physics and optics. We show that photonic compound ring resonators can possess large anomalous GVD at any desirable wavelength, even if each individual resonator is characterized with normal GVD. Based on this concept we design a mode locked frequency comb with thin-film silicon nitride compound ring resonators in the vicinity of Rubidium D1 line (794.6nm) and propose to use this optical comb as a flywheel for chip-scale optical clocks.


title: Influence of spin polarizability on liquid gas phase transition in the nuclear matter

abstract: In this paper, we investigate the liquid gas phase transition for the spin polarized nuclear matter. Applying the lowest order constrained variational (LOCV) method, and using two microscopic potentials, $AV_{18}$ and $UV_{14}$+TNI, we calculate the free energy, equation of state, order parameter, entropy, heat capacity and compressibility to derive the critical properties of spin polarized nuclear matter. Our results indicate that for the spin polarized nuclear matter, the second order phase transition takes place at lower temperatures with respect to the unpolarized one. It is also shown that the critical temperature of our spin polarized nuclear matter with a specific value of spin polarization parameter is in good agreement with the experimental result.


title: Fragmentation contributions to hadroproduction of prompt $J/\psi$, $\chi_{cJ}$, and $\psi(2S)$ states

abstract: We compute fragmentation corrections to hadroproduction of the quarkonium states $J/\psi$, $\chi_{cJ}$, and $\psi(2S)$ at leading power in $m_c^2/p_T^2$, where $m_c$ is the charm-quark mass and $p_T$ is the quarkonium transverse momentum. The computation is carried out in the framework of nonrelativistic QCD. We include corrections to the parton-production cross sections through next-to-leading order in the strong coupling $\alpha_s$ and corrections to the fragmentation functions through second order in $\alpha_s$. We also sum leading logarithms of $p_T^2/m_c^2$ to all orders in perturbation theory. We find that, when we combine these leading-power fragmentation corrections with fixed-order calculations through next-to-leading order in $\alpha_s$, we are able to obtain good fits for $p_T\geq 10$ GeV to hadroproduction cross sections that were measured at the Tevatron and the LHC. Using values for the nonperturbative long-distance matrix elements that we extract from the cross-section fits, we make predictions for the polarizations of the quarkonium states. We obtain good agreement with measurements of the polarizations, with the exception of the CDF Run II measurement of the prompt $J/\psi$ polarization, for which the agreement is only fair. In the predictions for the prompt-$J/\psi$ cross sections and polarizations, we take into account feeddown from the $\chi_{cJ}$ and $\psi(2S)$ states.


title: Fast Quantum Control for Weakly Nonlinear Qubits: On Two-Quadrature Adiabatic Gates

abstract: Adiabatic or slowly varying gate operations are typically required in order to remain within the qubit subspace in an anharmonic oscillator. However significant speed ups are possible by using the two quadrature derivative-removal-by-adiabatic-gate(DRAG) technique, where a second time derivative pulse component burns a spectral hole near an unwanted transition. It is shown here, that simultaneous optimization of the detuning and the pulse norm in addition, further reduces leakage errors and significantly improve gate fidelities. However, with this optimization accounting for the AC Stark shift, there is a low spectral weight pulse envelope regime, where DRAG is almost not needed and where the two state error fidelities are stable against pulse jitter. Explicit time evolution calculations are carried out in the lab frame for truncated multi-level Transmon qubit models obtained from a tight-binding model.


title: Reconciling Muon g-2, 125 GeV Higgs and Dark Matter in Gauge Mediation Models

abstract: We present a class of models in the framework of gauge mediation supersymmetry breaking where the standard model is supplemented by additional U(1) symmetry which acts only on the third generation fermions. The messenger fields carry non-trivial U(1) charge and are vector-like particles under this symmetry. This leads to additional contribution to the soft supersymmetry breaking mass terms for the third generation squarks and sleptons. In this framework we show that the muon g-2 anomaly, the observed 125 GeV Higgs boson mass and the detected relic dark matter abundance (gravitino in our case) can be simultaneously accommodated. The resolution of the muon g-2 anomaly, in particular, yields the result that the first two generation squark masses, as well the gluino mass, should be <~ 2.5 TeV, which will be tested at LHC14.


title: Probing the Orbital Origin of Conductance Oscillations in Atomic Chains

abstract: We investigate periodical oscillations in the conductance of suspended Au and Pt atomic chains during elongation under mechanical stress. Analysis of conductance and shot noise measurements reveals that the oscillations are mainly related to variations in a specific conduction channel as the chain undergoes transitions between zigzag and linear atomic configurations. The calculated local electronic structure shows that the oscillations originate from varying degrees of hybridization between the atomic orbitals along the chain as a function of the zigzag angle. These variations are highly dependent on the directionality and symmetry of the relevant orbitals, in agreement with the order-of-magnitude difference between the Pt and Au oscillation amplitudes observed in experiment. Our results demonstrate that the sensitivity of conductance to structural variations can be controlled by designing atomic-scale conductors in view of the directional interactions between atomic orbitals.


title: Helly-type theorems for the diameter

abstract: We study versions of Helly's theorem that guarantee that the intersection of a family of convex sets in $R^d$ has a large diameter. This includes colourful, fractional and $(p,q)$ versions of Helly's theorem. In particular, the fractional and $(p,q)$ versions work with conditions where the corresponding Helly theorem does not. We also include variants of Tverberg's theorem, B\'ar\'any's point selection theorem and the existence of weak epsilon-nets for convex sets with diameter estimates.


title: Proposal of room-temperature diamond maser

abstract: Lasers have revolutionized optical science and technology, but their microwave counterpart, maser, has not realized its great potential due to its demanding work conditions (high-vacuum for gas maser and liquid-helium temperature for solid-state maser). Room-temperature solid-state maser is highly desirable, but under such conditions the lifetimes of emitters (usually electron spins) are usually too short (~ns) for population inversion. The only room-temperature solid-state maser is based on a pentacene-doped p-terphenyl crystal, which has long spin lifetime (~0.1 ms). This maser, however, operates only in the pulse mode and the material is unstable. Here we propose room-temperature maser based on nitrogen-vacancy (NV) centres in diamond, which feature long spin lifetimes at room temperature (~10 ms), high optical pump efficiency, and material stability. We demonstrate that under readily accessible conditions, room-temperature diamond maser is feasible. Room-temperature diamond maser may facilitate a broad range of microwave technologies.


title: Computable randomness and monotonicity

abstract: We show that $z\in\R^n$ is computably random if and only if every computable monotone function on $\R^n$ is differentiable at $z$.


title: Nonlinear MHD waves in a Prominence Foot

abstract: We study nonlinear waves in a prominence foot using 2.5D MHD model motivated by recent high-resolution observations with Hinode/SOT in Ca~II emission of a prominence on October 10, 2012 showing highly dynamic small-scale motions in the prominence material. Observations of H$\alpha$ intensities and of Doppler shifts show similar propagating fluctuations. However the optically thick nature of the emission lines inhibits unique quantitative interpretation in terms of density. Nevertheless, we find evidence of nonlinear wave activity in the prominence foot by examining the relative magnitude of the fluctuation intensity ($\delta I/I\sim \delta n/n$). The waves are evident as significant density fluctuations that vary with height, and apparently travel upward from the chromosphere into the prominence material with quasi-periodic fluctuations with typical period in the range of 5-11 minutes, and wavelengths $\sim <$2000 km. Recent Doppler shift observations show the transverse displacement of the propagating waves. The magnetic field was measured with THEMIS instrument and was found to be 5-14 G. For the typical prominence density the corresponding fast magnetosonic speed is $\sim$20 km s$^{-1}$, in qualitative agreement with the propagation speed of the detected waves. The 2.5D MHD numerical model is constrained with the typical parameters of the prominence waves seen in observations. Our numerical results reproduce the nonlinear fast magnetosonic waves and provide strong support for the presence of these waves in the prominence foot. We also explore gravitational MHD oscillations of the heavy prominence foot material supported by dipped magnetic field structure.


title: A Six-Planet System Orbiting HD 219134

abstract: We present new, high-precision Doppler radial velocity (RV) data sets for the nearby K3V star HD 219134. The data include 175 velocities obtained with the HIRES Spectrograph at the Keck I Telescope, and 101 velocities obtained with the Levy Spectrograph at the Automated Planet Finder Telescope (APF) at Lick Observatory. Our observations reveal six new planetary candidates, with orbital periods of P=3.1, 6.8, 22.8, 46.7, 94.2 and 2247 days, spanning masses of msini=3.8, 3.5, 8.9, 21.3, 10.8 and 108 M_earth respectively. Our analysis indicates that the outermost signal is unlikely to be an artifact induced by stellar activity. In addition, several years of precision photometry with the T10 0.8~m automatic photometric telescope (APT) at Fairborn Observatory demonstrated a lack of brightness variability to a limit of ~0.0002 mag, providing strong support for planetary-reflex motion as the source of the radial velocity variations. The HD 219134 system, with its bright (V=5.6) primary provides an excellent opportunity to obtain detailed orbital characterization (and potentially follow-up observations) of a planetary system that resembles many of the multiple-planet systems detected by Kepler, and which are expected to be detected by NASA's forthcoming TESS Mission and by ESA's forthcoming PLATO Mission.


title: Clinical trial design enabling {\epsilon}-optimal treatment rules

abstract: Medical research has evolved conventions for choosing sample size in randomized clinical trials that rest on the theory of hypothesis testing. Bayesians have argued that trials should be designed to maximize subjective expected utility in settings of clinical interest. This perspective is compelling given a credible prior distribution on treatment response, but Bayesians have struggled to provide guidance on specification of priors. We use the frequentist statistical decision theory of Wald (1950) to study design of trials under ambiguity. We show that {\epsilon}-optimal rules exist when trials have large enough sample size. An {\epsilon}-optimal rule has expected welfare within {\epsilon} of the welfare of the best treatment in every state of nature. Equivalently, it has maximum regret no larger than {\epsilon}. We consider trials that draw predetermined numbers of subjects at random within groups stratified by covariates and treatments. The principal analytical findings are simple sufficient conditions on sample sizes that ensure existence of {\epsilon}-optimal treatment rules when outcomes are bounded. These conditions are obtained by application of Hoeffding (1963) large deviations inequalities to evaluate the performance of empirical success rules.


title: Uo-convergence and its applications to Ces\`aro means in Banach lattices

abstract: A net $(x_\alpha)$ in a vector lattice $X$ is said to uo-converge to $x$ if $|x_\alpha-x|\wedge u\xrightarrow{\rm o}0$ for every $u\ge 0$. In the first part of this paper, we study some functional-analytic aspects of uo-convergence. We prove that uo-convergence is stable under passing to and from regular sublattices. This fact leads to numerous applications presented throughout the paper. In particular, it allows us to improve several results in [26,27]. In the second part, we use uo-convergence to study convergence of Ces\`aro means in Banach lattices. In particular, we establish an intrinsic version of Koml\'os' Theorem, which extends the main results of [35,16,31] in a uniform way. We also develop a new and unified approach to Banach-Saks properties and Banach-Saks operators based on uo-convergence. This approach yields, in particular, short direct proofs of several results in [21,24,25].


title: Free and based path groupoids

abstract: We give an explicit description of the free path and loop groupoids in the Morita bicategory of translation topological groupoids. We prove that the free path groupoid of a discrete group acting on a topological space $X$ is a translation groupoid given by the same group acting on the topological path space $X^I$. We give a detailed description of based path and loop groupoids and show that both are equivalent to topological spaces. We also establish the notion of homotopy and fibration in this context.


title: Gauge invariant fluctuations of the metric during inflation from new scalar-tensor Weyl-Integrable gravity model

abstract: We investigate gauge invariant scalar fluctuations of the metric during inflation in a non-perturbative formalism in the framework of a recently introduced scalar-tensor theory of gravity formulated on a Weyl-Integrable geometry. We found that the Weyl scalar field can play the role of the inflaton field in this theory. As an application we study the case of a power law inflation. In this case the quasi-scale invariance of the spectrum for scalar fluctuations of the metric is achieved for determined values of the $\omega$ parameter of the scalar-tensor theory. In our formalism the physical inflaton field has a geometrical origin.


title: The SDSS-III BOSS quasar lens survey: discovery of thirteen gravitationally lensed quasars

abstract: We report the discovery of 13 confirmed two-image quasar lenses from a systematic search for gravitationally lensed quasars in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We adopted a methodology similar to that used in the SDSS Quasar Lens Search (SQLS). In addition to the confirmed lenses, we report 11 quasar pairs with small angular separations ($\lesssim$2") confirmed from our spectroscopy, which are either projected pairs, physical binaries, or possibly quasar lens systems whose lens galaxies have not yet been detected. The newly discovered quasar lens system, SDSS J1452+4224 at zs$\approx$4.8 is one of the highest redshift multiply imaged quasars found to date. Furthermore, we have over 50 good lens candidates yet to be followed up. Owing to the heterogeneous selection of BOSS quasars, the lens sample presented here does not have a well-defined selection function.


title: Algorithmic Methods to Infer the Evolutionary Trajectories in Cancer Progression

abstract: The genomic evolution inherent to cancer relates directly to a renewed focus on the voluminous next generation sequencing (NGS) data, and machine learning for the inference of explanatory models of how the (epi)genomic events are choreographed in cancer initiation and development. However, despite the increasing availability of multiple additional -omics data, this quest has been frustrated by various theoretical and technical hurdles, mostly stemming from the dramatic heterogeneity of the disease. In this paper, we build on our recent works on "selective advantage" relation among driver mutations in cancer progression and investigate its applicability to the modeling problem at the population level. Here, we introduce PiCnIc (Pipeline for Cancer Inference), a versatile, modular and customizable pipeline to extract ensemble-level progression models from cross-sectional sequenced cancer genomes. The pipeline has many translational implications as it combines state-of-the-art techniques for sample stratification, driver selection, identification of fitness-equivalent exclusive alterations and progression model inference. We demonstrate PiCnIc's ability to reproduce much of the current knowledge on colorectal cancer progression, as well as to suggest novel experimentally verifiable hypotheses.


title: Analysis of A Splitting Approach for the Parallel Solution of Linear Systems on GPU Cards

abstract: We discuss an approach for solving sparse or dense banded linear systems ${\bf A} {\bf x} = {\bf b}$ on a Graphics Processing Unit (GPU) card. The matrix ${\bf A} \in {\mathbb{R}}^{N \times N}$ is possibly nonsymmetric and moderately large; i.e., $10000 \leq N \leq 500000$. The ${\it split\ and\ parallelize}$ (${\tt SaP}$) approach seeks to partition the matrix ${\bf A}$ into diagonal sub-blocks ${\bf A}_i$, $i=1,\ldots,P$, which are independently factored in parallel. The solution may choose to consider or to ignore the matrices that couple the diagonal sub-blocks ${\bf A}_i$. This approach, along with the Krylov subspace-based iterative method that it preconditions, are implemented in a solver called ${\tt SaP::GPU}$, which is compared in terms of efficiency with three commonly used sparse direct solvers: ${\tt PARDISO}$, ${\tt SuperLU}$, and ${\tt MUMPS}$. ${\tt SaP::GPU}$, which runs entirely on the GPU except several stages involved in preliminary row-column permutations, is robust and compares well in terms of efficiency with the aforementioned direct solvers. In a comparison against Intel's ${\tt MKL}$, ${\tt SaP::GPU}$ also fares well when used to solve dense banded systems that are close to being diagonally dominant. ${\tt SaP::GPU}$ is publicly available and distributed as open source under a permissive BSD3 license.


title: Finitely based sets of 2-limited block-2-simple words

abstract: Let $\mathfrak A$ be an alphabet and $W$ be a set of words in the free monoid ${\mathfrak A}^*$. Let $S(W)$ denote the Rees quotient over the ideal of ${\mathfrak A}^*$ consisting of all words that are not subwords of words in $W$. A set of words $W$ is called {\em finitely based} if the monoid $S(W)$ is finitely based. A word $\bf u$ is called 2-limited if each variable occurs in $\bf u$ at most twice. A {\em block} of a word $\bf u$ is a maximal subword of $\bf u$ that does not contain any linear variables. We say that a word $\bf u$ is {\em block-2-simple} if each block of $\bf u$ involves at most two distinct variables. We provide an algorithm that recognizes finitely based sets of words among sets of 2-limited block-2-simple words. We also present new sufficient conditions under which a set of words is non-finitely based.


title: Estimation of direct laser acceleration in laser wakefield accelerators using particle-in-cell simulations

abstract: Many current laser wakefield acceleration (LWFA) experiments are carried out in a regime where the laser pulse length is on the order of or longer than the wake wavelength and where ionization injection is employed to inject electrons into the wake. In these experiments, the trapped electrons will co-propagate with the longitudinal wakefield and the transverse laser field. In this scenario, the electrons can gain a significant amount of energy from both the direct laser acceleration (DLA) mechanism as well as the usual LWFA mechanism. Particle-in-cell (PIC) codes are frequently used to discern the relative contribution of these two mechanisms. However, if the longitudinal resolution used in the PIC simulations is inadequate, it can produce numerical heating that can overestimate the transverse motion, which is important in determining the energy gain due to DLA. We have therefore carried out a systematic study of this LWFA regime by varying the longitudinal resolution of PIC simulations from the standard, best-practice resolution of 30 points per laser wavelength to four times that value and then examining the energy gain characteristics of both the highest-energy electrons and the bulk electrons. By calculating the contribution of DLA to the final energies of the electrons produced from the LWFA, we find that although the transverse momentum and oscillation radii are over-estimated in the lower-resolution simulations, this over-estimation does not lead to artificial energy gain by DLA. Rather, the DLA contribution to the highest-energy electrons is larger in the higher-resolution cases because the DLA resonance is better maintained. Thus, even at the highest longitudinal resolutions, DLA contributes a significant portion of the energy gained by the highest-energy electrons and also contributes to accelerating the bulk of the charge in the electron beam produced by the LWFA.


title: Suboptimal Stabilizing Controllers for Linearly Solvable System

abstract: This paper presents a novel method to synthesize stochastic control Lyapunov functions for a class of nonlinear, stochastic control systems. In this work, the classical nonlinear Hamilton-Jacobi-Bellman partial differential equation is transformed into a linear partial differential equation for a class of systems with a particular constraint on the stochastic disturbance. It is shown that this linear partial differential equation can be relaxed to a linear differential inclusion, allowing for approximating polynomial solutions to be generated using sum of squares programming. It is shown that the resulting solutions are stochastic control Lyapunov functions with a number of compelling properties. In particular, a-priori bounds on trajectory suboptimality are shown for these approximate value functions. The result is a technique whereby approximate solutions may be computed with non-increasing error via a hierarchy of semidefinite optimization problems.


title: Bilinear quadratures for inner products

abstract: A bilinear quadrature numerically evaluates a continuous bilinear map, such as the $L^2$ inner product, on continuous $f$ and $g$ belonging to known finite-dimensional function spaces. Such maps arise in Galerkin methods for differential and integral equations. The construction of bilinear quadratures over arbitrary domains in $\mathbb{R}^d$ is presented. In one dimension, integration rules of this type include Gaussian quadrature for polynomials and the trapezoidal rule for trigonometric polynomials as special cases. A numerical procedure for constructing bilinear quadratures is developed and validated.


title: Testing the validity of the Kirkwood approximation using an extended Sznajd model

abstract: We revisit the deduction of the exit probability of the one dimensional Sznajd model through the Kirkwood approximation [F. Slanina et al., Europhys. Lett. 82, 18006 (2008)]. This approximation is peculiar in that in spite of the agreement with simulation results [F. Slanina et al., Europhys. Lett. 82, 18006 (2008), R. Lambiotte and S. Redner, Europhys. Lett. 82, 18007 (2008), A. M. Timpanaro and C. P. C. Prado, 89, 052808 (2014)] the hypothesis about the correlation lenghts behind it are inconsistent and fixing these inconsistencies leads to the same results as a simple mean field. We use an extended version of the Sznajd model to test the Kirkwood approximation in a wider context. This model includes the voter, Sznajd and "United we Stand, Divided we Fall" (USDF) models [R. A. Holley and T. M. Liggett, Ann. Prob. 3, 643 (1975), K. Sznajd-Weron and J. Sznajd, Int. Journ. Mod. Phys. C 11, 1157 (2000)] as different parameter combinations, meaning that some analytical results from these models can be used to evaluate the performance of the Kirkwood approximation. We also compare the predicted exit probability with simulation results for networks with $10^3$ sites. The results show clearly the regions in parameter space where the approximation gives accurate predictions, as well as where it starts failing, leading to a better understanding of its reliability.


title: Multidimensional Butterfly Factorization

abstract: This paper introduces the multidimensional butterfly factorization as a data-sparse representation of multidimensional kernel matrices that satisfy the complementary low-rank property. This factorization approximates such a kernel matrix of size $N\times N$ with a product of $\O(\log N)$ sparse matrices, each of which contains $\O(N)$ nonzero entries. We also propose efficient algorithms for constructing this factorization when either (i) a fast algorithm for applying the kernel matrix and its adjoint is available or (ii) every entry of the kernel matrix can be evaluated in $\O(1)$ operations. For the kernel matrices of multidimensional Fourier integral operators, for which the complementary low-rank property is not satisfied due to a singularity at the origin, we extend this factorization by combining it with either a polar coordinate transformation or a multiscale decomposition of the integration domain to overcome the singularity. Numerical results are provided to demonstrate the efficiency of the proposed algorithms.


title: Spin relaxation near a ferromagnetic transition

abstract: We study spin relaxation in dilute magnetic semiconductors near a ferromagnetic transition, where spin fluctuations become strong. An enhancement in the scattering rate of itinerant carriers from the spin fluctuations of localized impurities leads to a change in the dominant spin relaxation mechanism from Dyakonov-Perel to spin flips in scattering. On the ferromagnetic side of the transition, we show that due to the presence of two magnetic components -- the itinerant carriers and the magnetic impurities -- with different gyromagnetic ratios, the relaxation rate of the total magnetization can be quite different from the relaxation rate of the spin. Following a disturbance of the equilibrium magnetization, the spin is initially redistributed between the two components to restore the equilibrium magnetization. It is only on a longer time scale, controlled by the spin-orbit interaction, that the total spin itself relaxes to its equilibrium state.


title: Algorithms for Linear Bandits on Polyhedral Sets

abstract: We study stochastic linear optimization problem with bandit feedback. The set of arms take values in an $N$-dimensional space and belong to a bounded polyhedron described by finitely many linear inequalities. We provide a lower bound for the expected regret that scales as $\Omega(N\log T)$. We then provide a nearly optimal algorithm and show that its expected regret scales as $O(N\log^{1+\epsilon}(T))$ for an arbitrary small $\epsilon >0$. The algorithm alternates between exploration and exploitation intervals sequentially where deterministic set of arms are played in the exploration intervals and greedily selected arm is played in the exploitation intervals. We also develop an algorithm that achieves the optimal regret when sub-Gaussianity parameter of the noise term is known. Our key insight is that for a polyhedron the optimal arm is robust to small perturbations in the reward function. Consequently, a greedily selected arm is guaranteed to be optimal when the estimation error falls below some suitable threshold. Our solution resolves a question posed by Rusmevichientong and Tsitsiklis (2011) that left open the possibility of efficient algorithms with asymptotic logarithmic regret bounds. We also show that the regret upper bounds hold with probability $1$. Our numerical investigations show that while theoretical results are asymptotic the performance of our algorithms compares favorably to state-of-the-art algorithms in finite time as well.


title: Quantized Massive MU-MIMO-OFDM Uplink

abstract: Coarse quantization at the base station (BS) of a massive multi-user (MU) multiple-input multiple-output (MIMO) wireless system promises significant power and cost savings. Coarse quantization also enables significant reductions of the raw analog-to-digital converter (ADC) data that must be transferred from a spatially-separated antenna array to the baseband processing unit. The theoretical limits as well as practical transceiver algorithms for such quantized MU-MIMO systems operating over frequency-flat, narrowband channels have been studied extensively. However, the practically relevant scenario where such communication systems operate over frequency-selective, wideband channels is less well understood. This paper investigates the uplink performance of a quantized massive MU-MIMO system that deploys orthogonal frequency-division multiplexing (OFDM) for wideband communication. We propose new algorithms for quantized maximum a-posteriori (MAP) channel estimation and data detection, and we study the associated performance/quantization trade-offs. Our results demonstrate that coarse quantization (e.g., four to six bits, depending on the ratio between the number of BS antennas and the number of users) in massive MU-MIMO-OFDM systems entails virtually no performance loss compared to the infinite-precision case at no additional cost in terms of baseband processing complexity.


title: Perturbative results for two and three particle threshold energies in finite volume

abstract: We calculate the energy of the state closest to threshold for two and three identical, spinless particles confined to a cubic spatial volume with periodic boundary conditions and with zero total momentum in the finite-volume frame. The calculation is performed in relativistic quantum field theory with particles coupled via a $\lambda \phi^4$ interaction, and we work through order $\lambda^3$. The energy shifts begin at ${\cal O}(1/L^3)$, and we keep subleading terms proportional to $1/L^4$, $1/L^5$ and $1/L^6$. These terms allow a non-trivial check of the results obtained from quantization conditions that hold for arbitrary interactions, namely that of L\"uscher for two particles and our recently developed formalism for three particles. We also compare to previously obtained results based on non-relativistic quantum mechanics.


title: Spin-Orbit Misalignment of Two-Planet-System KOI-89 Via Gravity Darkening

abstract: We constrain the true spin-orbit alignment of the KOI-89 system by numerically fitting the two \emph{Kepler} photometric lightcurves produced by transiting planets KOI-89.01 and KOI-89.02. The two planets have periods of 84.69 days and 207.58 days, respectively. We find that the two bodies are low-density giant planets with radii $0.45 \pm 0.03~\mathrm{R_{jup}}$ and $0.43 \pm 0.05~\mathrm{R_{jup}}$ and spin-orbit misalignments $72^{\circ} \pm 3^\circ$ and $73^{\circ+11}_{-5}$, respectively. Via dynamic stability tests we demonstrate the general trend of higher system stability with the two planets close to mutual alignment and estimate their coalignment angle to $20^\circ \pm 20^\circ$ -- i.e. the planets are misaligned with the star but may be aligned with each other. From these results, we limit KOI-89's misalignment mechanisms to star-disk-binary interactions, disk warping via planet-disk interactions, planet-planet scattering, Kozai resonance, or internal gravity waves.


title: Scattering of entangled two-photon states

abstract: We consider the scattering of entangled two-photon states from collections of small particles. We also study the related Mie problem of scattering from a sphere. In both cases, we calculate the entropy of entanglement and investigate the influence of the entanglement of the incident field on the entanglement of the scattered field.


title: The stable homotopy classification of $(n-1)$-connected $(n+4)$-dimensional polyhedra with 2 torsion free homology

abstract: In this paper, we study the stable homotopy types of $\mathbf{F}^4_{n(2)}$-polyhedra, i.e., $(n-1)$-connected, at most $(n+4)$-dimensional polyhedra with 2-torsion free homologies. We are able to classify the indecomposable $\mathbf{F}^4_{n(2)}$-polyhedra. The proof relies on the matrix problem technique which was developed in the classification of representaions of algebras and applied to homotopy theory by Baues and Drozd.


title: Beam specific planning target volumes incorporating 4DCT for pencil beam scanning proton therapy of thoracic tumors

abstract: The purpose of this study is to determine whether organ sparing and target coverage can be simultaneously maintained for pencil beam scanning (PBS) proton therapy treatment of thoracic tumors in the presence of motion, stopping power uncertainties and patient setup variations. Ten consecutive patients that were previously treated with proton therapy to 66.6/1.8 Gy (RBE) using double scattering (DS) were replanned with PBS. Minimum and maximum intensity images from 4DCT were used to introduce flexible smearing in the determination of the beam specific PTV (BSPTV). Datasets from eight 4DCT phases, using +-3% uncertainty in stopping power, and +-3 mm uncertainty in patient setup in each direction were used to create 8X12X10=960 PBS plans for the evaluation of ten patients. Plans were normalized to provide identical coverage between DS and PBS. The average lung V20, V5, and mean doses were reduced from 29.0%, 35.0%, and 16.4 Gy with DS to 24.6%, 30.6%, and 14.1 Gy with PBS, respectively. The average heart V30 and V45 were reduced from 10.4% and 7.5% in DS to 8.1% and 5.4% for PBS, respectively. Furthermore, the maximum spinal cord, esophagus and heart dose were decreased from 37.1 Gy, 71.7 Gy and 69.2 Gy with DS to 31.3 Gy, 67.9 Gy and 64.6 Gy with PBS. The conformity index (CI), homogeneity index (HI), and global maximal dose were improved from 3.2, 0.08, 77.4 Gy with DS to 2.8, 0.04 and 72.1 Gy with PBS. All differences are statistically significant, with p values <0.05, with the exception of the heart V45 (p= 0.146). PBS with BSPTV achieves better organ sparing and improves target coverage using a repainting method for the treatment of thoracic tumors. Incorporating motion-related uncertainties is essential in maintaining marginal coverage and homogenous dose of treatment targets.


title: Particles Acceleration in Converged Two Shocks

abstract: Observations show that there is a proton spectral "break" with E$_{break}$ at 1-10MeV in some large CME-driven shocks. Theoretical model usually attribute this phenomenon to a diffusive shock acceleration. However, the underlying physics of the shock acceleration still remains uncertain. Although previous numerical models can hardly predict this "break" due to either high computational expense or shortcomings of current models, the present paper focuses on simulating this energy spectrum in converged two shocks by Monte Carlo numerical method. Considering the Dec 13 2006 CME-driven shock interaction with an Earth bow shock, we examine whether the energy spectral "break" could occur on an interaction between two shocks. As result, we indeed obtain the maximum proton energy up to 10MeV, which is the premise to investigate the existence of the energy spectral "break". Unexpectedly, we further find a proton spectral "break" appears distinctly at the energy $\sim$5MeV.


title: A note on the dynamic dominant resource fairness mechanism

abstract: Multi-resource fair allocation has beena hot topic of resource allocation. Most recently, a dynamic dominant resource fairness (DRF) mechanism is proposed for dynamic multi-resource fair allocation. In this paper, we prove that the competitive ratio of the dynamic DRF mechanism is the reciprocal of the number of resource types, for two different objectives. Moreover, we develop a linear-time algorithm to find a dynamic DRF solution at each step.


title: Observation of quantum dynamical oscillations of ultracold atoms in the F and D bands of an optical lattice

abstract: We report the observation of quantum dynamical oscillations of ultracold atomic gases in the F and D bands of a single-well optical lattice. We are able to control the Bragg reflections at the Brillouin zone edge up to the third order. As a result, we can switch the quantum dynamics from oscillations across both the F and D bands to oscillations only within the F-band. Our capability to observe these remarkable oscillations comes from the innovative non-adiabatic technique which allows us to load ultracold atoms efficiently to the G-band of an optical lattice.


title: Destroying coherence in high temperature superconductors with current flow

abstract: The loss of single-particle coherence going from the superconducting state to the normal state in underdoped cuprates is a dramatic effect that has yet to be understood. Here, we address this issue by performing angle resolved photoemission spectroscopy (ARPES) measurements in the presence of a transport current. We find that the loss of coherence is associated with the development of an onset in the resistance, in that well before the midpoint of the transition is reached, the sharp peaks in the ARPES spectra are completely suppressed. Since the resistance onset is a signature of phase fluctuations, this implies that the loss of single-particle coherence is connected with the loss of long-range phase coherence.


title: Spin one matter fields

abstract: Spin-one matter fields are relevant both for the description of hadronic states and as potential extensions of the Standard Model. In this work we present a formalism for the description of massive spin-one fields transforming in the $(1,0)\oplus(0,1)$ representation of the Lorentz group, based on the covariant projection onto parity eigenspaces and Poincar\'e orbits. The formalism yields a constrained dynamics. We solve the constraints and perform the canonical quantization accordingly. This formulation uses the recent construction of a parity-based covariant basis for matrix operators acting on the $(j,0)\oplus(0,j) $ representations. The algebraic properties of the covariant basis play an important role in solving the constraints and allowing the canonical quantization of the theory. We study the chiral structure of the theory and conclude that it is not chirally symmetric in the massless limit, hence it is not possible to have chiral gauge interactions. However, spin-one matter fields can have vector gauge interactions. Also, the dimension of the field makes self-interactions naively renormalizable. Using the covariant basis, we classify all possible naively renormalizable self-interaction terms.


title: Novel Collective Phenomena in High-Energy Proton-Proton and Proton-Nucleus Collisions

abstract: The observation of long-range collective correlations for particles emitted in high-multiplicity pp and pPb collisions has opened up new opportunities of investigating novel high-density QCD phenomena in small colliding systems. We review experimental results related to the studies of collective phenomena in small systems from RHIC and the LHC over the past several years. Latest development in theoretical interpretations motivated by different frameworks are also reviewed, and confronted with the experimental data. Perspectives on possible future directions are discussed, with the aim of further exploring the rich emergent QCD phenomena.


title: Kosambi-Cartan-Chern (KCC) theory for higher order dynamical systems

abstract: The Kosambi-Cartan-Chern (KCC) theory represents a powerful mathematical method for the investigation of the properties of dynamical systems. The KCC theory introduces a geometric description of the time evolution of a dynamical system, with the solution curves of the dynamical system described by methods inspired by the theory of geodesics in a Finsler spaces. The evolution of a dynamical system is geometrized by introducing a non-linear connection, which allows the construction of the KCC covariant derivative, and of the deviation curvature tensor. In the KCC theory the properties of any dynamical system are described in terms of five geometrical invariants, with the second one giving the Jacobi stability of the system. Usually, the KCC theory is formulated by reducing the dynamical evolution equations to a set of second order differential equations. In the present paper we introduce and develop the KCC approach for dynamical systems described by systems of arbitrary $n$-dimensional first order differential equations. We investigate in detail the properties of the $n$-dimensional autonomous dynamical systems, as well as the relationship between the linear stability and the Jacobi stability. As a main result we find that only even-dimensional dynamical systems can exhibit both Jacobi stability and instability behaviors, while odd-dimensional dynamical systems are always Jacobi unstable, no matter their Lyapunov stability. As applications of the developed formalism we consider the geometrization and the study of the Jacobi stability of the complex dynamical networks, and of the $\Lambda $-Cold Dark Matter cosmological models, respectively.


title: Pattern Avoidance for Random Permutations

abstract: Using techniques from Poisson approximation, we prove explicit error bounds on the number of permutations that avoid any pattern. Most generally, we bound the total variation distance between the joint distribution of pattern occurrences and a corresponding joint distribution of independent Bernoulli random variables, which as a corollary yields a Poisson approximation for the distribution of the number of occurrences of any pattern. We also investigate occurrences of consecutive patterns in random Mallows permutations, of which uniform random permutations are a special case. These bounds allow us to estimate the probability that a pattern occurs any number of times and, in particular, the probability that a random permutation avoids a given pattern.


title: Bayesian Estimators in Uncertain Nested Error Regression Models

abstract: Nested error regression models are useful tools for analysis of grouped data, especially in the case of small area estimation. This paper suggests a nested error regression model using uncertain random effects in which the random effect in each area is expressed as a mixture of a normal distribution and a positive mass at $0$. For estimation of the model parameters and prediction of the random effects, an objective Bayesian inference is proposed by setting non-informative prior distributions on the model parameters. Under mild sufficient conditions, it is shown that the posterior distribution is proper and the posterior variances are finite, confirming the validity of posterior inference. To generate samples from the posterior distribution, we provide the Gibbs sampling method with familiar forms for all the full conditional distributions. This paper also addresses the problem of predicting finite population means, and a sampling-based method is suggested to tackle this issue. Finally, the proposed model is compared with the conventional nested error regression model through simulation and empirical studies.


title: Super-Resolution Off the Grid

abstract: Super-resolution is the problem of recovering a superposition of point sources using bandlimited measurements, which may be corrupted with noise. This signal processing problem arises in numerous imaging problems, ranging from astronomy to biology to spectroscopy, where it is common to take (coarse) Fourier measurements of an object. Of particular interest is in obtaining estimation procedures which are robust to noise, with the following desirable statistical and computational properties: we seek to use coarse Fourier measurements (bounded by some cutoff frequency); we hope to take a (quantifiably) small number of measurements; we desire our algorithm to run quickly. Suppose we have k point sources in d dimensions, where the points are separated by at least \Delta from each other (in Euclidean distance). This work provides an algorithm with the following favorable guarantees: - The algorithm uses Fourier measurements, whose frequencies are bounded by O(1/\Delta) (up to log factors). Previous algorithms require a cutoff frequency which may be as large as {\Omega}( d/\Delta). - The number of measurements taken by and the computational complexity of our algorithm are bounded by a polynomial in both the number of points k and the dimension d, with no dependence on the separation \Delta. In contrast, previous algorithms depended inverse polynomially on the minimal separation and exponentially on the dimension for both of these quantities. Our estimation procedure itself is simple: we take random bandlimited measurements (as opposed to taking an exponential number of measurements on the hyper-grid). Furthermore, our analysis and algorithm are elementary (based on concentration bounds for sampling and the singular value decomposition).


title: Unit-Regularity of Regular Nilpotent Elements

abstract: Let $a$ be a regular element of a ring $R$. If either $K:=\rm{r}_R(a)$ has the exchange property or every power of $a$ is regular, then we prove that for every positive integer $n$ there exist decompositions $$ R_R = K \oplus X_n \oplus Y_n = E_n \oplus X_n \oplus aY_n,$$ where $Y_n \subseteq a^nR$ and $E_n \cong R/aR$. As applications we get easier proofs of the results that a strongly $\pi$-regular ring has stable range one and also that a strongly $\pi$-regular element whose every power is regular is unit-regular.


title: Representation-tame algebras need not be homologically tame

abstract: We show that, also within the class of representation-tame finite dimensional algebras $\Lambda$, the big left finitistic dimension of $\Lambda$ may be strictly larger than the little. In fact, the discrepancies $Findim \Lambda - findim \Lambda$ need not even be bounded for special biserial algebras which constitute one of the (otherwise) most thoroughly understood classes of tame algebras. More precisely: For every positive integer $r$, we construct a special biserial algebra $\Lambda$ with the property that $findim \Lambda = r + 1$, while $Findim \Lambda = 2r + 1$. In particular, there are infinite dimensional representations of $\Lambda$ which have finite projective dimension, while not being direct limits of {\it finitely generated\/} representations of finite projective dimension.


title: A Revisit of Infinite Population Models for Evolutionary Algorithms on Continuous Optimization Problems

abstract: Infinite population models are important tools for studying population dynamics of evolutionary algorithms. They describe how the distributions of populations change between consecutive generations. In general, infinite population models are derived from Markov chains by exploiting symmetries between individuals in the population and analyzing the limit as the population size goes to infinity. In this paper, we study the theoretical foundations of infinite population models of evolutionary algorithms on continuous optimization problems. First, we show that the convergence proofs in a widely cited study were in fact problematic and incomplete. We further show that the modeling assumption of exchangeability of individuals cannot yield the transition equation. Then, in order to analyze infinite population models, we build an analytical framework based on convergence in distribution of random elements which take values in the metric space of infinite sequences. The framework is concise and mathematically rigorous. It also provides an infrastructure for studying the convergence of the stacking of operators and of iterating the algorithm which previous studies failed to address. Finally, we use the framework to prove the convergence of infinite population models for the mutation operator and the $k$-ary recombination operator. We show that these operators can provide accurate predictions for real population dynamics as the population size goes to infinity, provided that the initial population is identically and independently distributed.


title: A new result of the scaling law of weighted L1 minimization

abstract: This paper study recovery conditions of weighted L1 minimization for signal reconstruction from compressed sensing measurements. A sufficient condition for exact recovery by using the general weighted L1 minimization is derived, which builds a direct relationship between the weights and the recoverability. Simulation results indicates that this sufficient condition provides a precise prediction of the scaling law for the weighted L1 minimization.


title: Coalescence of spectral singularities and phase diagrams for one-dimensional PT symmetric photonic crystals

abstract: Non-Hermitian systems with parity-time (PT) symmetric complex potentials can exhibit a phase transition when the degree of non-Hermiticity is increased. Two eigenstates coalesce at a transition point, which is known as the exceptional point (EP) for a discrete spectrum and spectral singularity for a continuous spectrum. The existence of an EP is known to give rise to a great variety of novel behaviors in various fields of physics. In this work, we study the complex band structures of one-dimensional photonic crystals with PT symmetric complex potentials by setting up a Hamiltonian using the Bloch states of the photonic crystal without loss or gain as a basis. As a function of the degree of non-Hermiticity, two types of PT symmetry transitions are found. One is that a PT-broken phase can re-enter into a PT-exact phase at a higher degree of non-Hermiticity. The other is that two spectral singularities, one originating from the Brillouin zone center and the other from the Brillouin zone boundary, can coalesce at some k-point in the interior of the Brillouin zone and create a singularity of higher order. Furthermore, we can induce a band inversion by tuning the filling ratio of the photonic crystal, and we find that the geometric phases of the bands before and after the inversion are independent of the amount of non-Hermiticity as long as the PT-exact phase is not broken. The standard concept of topological transition can hence be extended to non-Hermitian systems.


title: Spin Dynamics and Relaxation in Graphene Dictated by Electron-hole Puddles

abstract: The understanding of spin dynamics and relaxation mechanisms in clean graphene and the upper time and length scales on which spin devices can operate are prerequisites to realizing graphene spintronic technologies. Here we theoretically reveal the nature of fundamental spin relaxation mechanisms in clean graphene on different substrates with spin-orbit Rashba fields as low as a few tens of micron eV. Spin lifetimes ranging from 50 picoseconds up to several nanoseconds are found to be dictated by substrate-induced electron-hole characteristics. A crossover in the spin relaxation mechanism from a Dyakonov-Perel type for SiO2 substrates to a broadening-induced dephasing for hBN substrates is described. The energy dependence of spin lifetimes, their ratio for spins pointing out-of-plane and in-plane, and the scaling with disorder provide a global picture about spin dynamics and relaxation in ultraclean graphene in presence of electron-hole puddles.


title: Mixed-ADC Massive MIMO Detectors: Performance Analysis and Design Optimization

abstract: Using a very low-resolution analog-to-digital convertor (ADC) unit at each antenna can remarkably reduce the hardware cost and power consumption of a massive multiple-input multiple-output (MIMO) system. However, such a pure low-resolution ADC architecture also complicates parameter estimation problems such as time/frequency synchronization and channel estimation. A mixed-ADC architecture, where most of the antennas are equipped with low-precision ADCs while a few antennas have full-precision ADCs, can solve these issues and actualize the potential of the pure low-resolution ADC architecture. In this paper, we present a unified framework to develop a family of detectors over the massive MIMO uplink system with the mixed-ADC receiver architecture by exploiting probabilistic Bayesian inference. As a basic setup, an optimal detector is developed to provide a minimum mean-squared-error (MMSE) estimate on data symbols. Considering the highly nonlinear steps involved in the quantization process, we also investigate the potential for complexity reduction on the optimal detector by postulating the common \emph{pseudo-quantization noise} (PQN) model. In particular, we provide asymptotic performance expressions including the MSE and bit error rate for the optimal and suboptimal MIMO detectors. The asymptotic performance expressions can be evaluated quickly and efficiently; thus, they are useful in system design optimization. We show that in the low signal-to-noise ratio (SNR) regime, the distortion caused by the PQN model can be ignored, whereas in the high-SNR regime, such distortion may cause 1-bit detection performance loss. The performance gap resulting from the PQN model can be narrowed by a small fraction of high-precision ADCs in the mixed-ADC architecture.


title: Error Gradient-based Variable-Lp Norm Constraint LMS Algorithm for Sparse System Identification

abstract: Sparse adaptive filtering has gained much attention due to its wide applicability in the field of signal processing. Among the main algorithm families, sparse norm constraint adaptive filters develop rapidly in recent years. However, when applied for system identification, most priori work in sparse norm constraint adaptive filtering suffers from the difficulty of adaptability to the sparsity of the systems to be identified. To address this problem, we propose a novel variable p-norm constraint least mean square (LMS) algorithm, which serves as a variant of the conventional Lp-LMS algorithm established for sparse system identification. The parameter p is iteratively adjusted by the gradient descent method applied to the instantaneous square error. Numerical simulations show that this new approach achieves better performance than the traditional Lp-LMS and LMS algorithms in terms of steady-state error and convergence rate.


title: Inserting Multiple Edges into a Planar Graph

abstract: Let $G$ be a connected planar (but not yet embedded) graph and $F$ a set of additional edges not yet in $G$. The {multiple edge insertion} problem (MEI) asks for a drawing of $G+F$ with the minimum number of pairwise edge crossings, such that the subdrawing of $G$ is plane. An optimal solution to this problem approximates the crossing number of the graph $G+F$. Finding an exact solution to MEI is NP-hard for general $F$, but linear time solvable for the special case of $|F|=1$ (SODA01, Algorithmica) or when all of $F$ are incident to a new vertex (SODA09). The complexity for general $F$ but with constant $k=|F|$ was open, but algorithms both with relative and absolute approximation guarantees have been presented (SODA11, ICALP11). We show that the problem is fixed parameter tractable (FPT) in $k$ for biconnected $G$, or if the cut vertices of $G$ have degrees bounded by a constant. We give the first exact algorithm for this problem; it requires only $O(|V(G)|)$ time for any constant $k$.


title: Optimal trading strategies - a time series approach

abstract: Motivated by recent advances in the spectral theory of auto-covariance matrices, we are led to revisit a reformulation of Markowitz' mean-variance portfolio optimization approach in the time domain. In its simplest incarnation it applies to a single traded asset and allows to find an optimal trading strategy which - for a given return - is minimally exposed to market price fluctuations. The model is initially investigated for a range of synthetic price processes, taken to be either second order stationary, or to exhibit second order stationary increments. Attention is paid to consequences of estimating auto-covariance matrices from small finite samples, and auto-covariance matrix cleaning strategies to mitigate against these are investigated. Finally we apply our framework to real world data.


title: Naked singularity and black hole formation in self-similar Einstein-scalar fields with exponential potentials

abstract: Motivated by cosmic censorship in general relativity and string theory, we extend Christodoulou's celebrated examples of naked singularity formation in the Einstein-massless scalar field system to include a positive or negative scalar potential of exponential forms, i.e., $V(\phi)=\pm\exp(2\phi/\kappa)$ with a parameter $\kappa$. Under spherical symmetry and a self-similar ansatz depending on $\kappa$, we derive a 3-dimensional autonomous system of first-order ordinary differential equations, which incorporates the equations for massless scalar fields as a special case. Local behavior of the phase space is studied analytically with global solutions constructed numerically. Within the 3-dimensional solution manifold, we observe, for the negative potentials, naked singularity formation from nonsingular initial data for $\kappa^2<1$. Meanwhile, transitions between solutions containing naked singularities and black holes are also identified. However, when the potential is taken positive, numerical evolutions result in formation of black holes, but not naked singularities.


title: Hierarchy of Spin Operators, Quantum Gates, Entanglement, Tensor Product and Eigenvalues

abstract: {\bf Abstract.} We show that two hierarchies of spin Hamilton operators admit the same spectrum. Both Hamilton operators play a central role for quantum gates in particular for the case spin-$\frac12$ and the case spin-1. The spin-$\frac12$, spin-1, spin-$\frac32$ and spin-2 cases are studied in detail. Entanglement and mutually unbiased bases of the eigenvectors is discussed. Two triple Hamilton operators are also investigated. Both are also admitting the same spectrum.


title: Examples of naked singularity formation in higher-dimensional Einstein-vacuum spacetimes

abstract: The vacuum Einstein equations in 5+1 dimensions are shown to admit solutions describing naked singularity formation in gravitational collapse from nonsingular asymptotically locally flat initial data that contain no trapped surface. We present a class of specific examples with topology $\mathbb{R}^{3+1} \times S^2$. Thanks to the Kaluza-Klein dimensional reduction, these examples are constructed by lifting continuously self-similar solutions of the 4-dimensional Einstein-scalar field system with a negative exponential potential. The latter solutions are obtained by solving a 3-dimensional autonomous system of first-order ordinary differential equations with a combined analytic and numerical approach. Their existence provides a new test-bed for weak cosmic censorship in higher-dimensional gravity. In addition, we point out that a similar attempt of lifting Christodoulou's naked singularity solutions of massless scalar fields fails to capture formation of naked singularities in 4+1 dimensions, due to a diverging Kretschmann scalar in the initial data.


title: The wave impedance of an atomically thin crystal

abstract: I propose an expression for the electromagnetic wave impedance of a two-dimensional atomic crystal, and I deduce the Fresnel coefficients in terms of this quantity. It is widely known that a two-dimensional crystal can absorb light, if its conductivity is different from zero. It is less emphasized that they can also store a certain amount of electromagnetic energy. The concept of impedance is useful to quantify this point.


title: Wavebreaking-associated transmitted emission of attosecond extreme-ultraviolet pulses from laser-driven overdense plasmas

abstract: We present a new mechanism of attosecond extreme-ultraviolet (XUV) pulses generation from a relativistic laser-driven overdense plasma surfaces in the wavebreaking regime. Through particle-in-cell simulations and analysis, we demonstrate that the observed ultrashort XUV emission for the parameters we considered is predominantly due to a strong plasma-density oscillation subsequent to wavebreaking. The coupling of the strong density variation and the transverse fields in the front surface layer gives rise to the transmitted emission with frequencies mainly around the local plasma frequency. This mechanism provides new insights into the scenarios of XUV generation from solid surfaces and the dynamics of laser-plasma interactions.


title: Physical observables from boundary artifacts: scalar glueball in Yang-Mills theory

abstract: By relating the functional averages of a generic scalar operator in simulations with Open (O) and Periodic (P) boundary conditions (BCs) respectively for $SU(3)$ lattice gauge theory, we show that the scalar glueball mass and the glueball to vacuum matrix element can be extracted very efficiently from the former. Numerical results are compared with those extracted from the two point function of the time slice energy density (both PBC and OBC). The scaling properties of the mass and the matrix element are studied with the help of Wilson (gradient) flow.


title: Adaptive low-rank approximation and denoised Monte-Carlo approach for high-dimensional Lindblad equations

abstract: We present a twofold contribution to the numerical simulation of Lindblad equations. First, an adaptive numerical approach to approximate Lindblad equations using low-rank dynamics is described: a deterministic low-rank approximation of the density operator is computed, and its rank is adjusted dynamically, using an on-the-fly estimator of the error committed when reducing the dimension. On the other hand, when the intrinsic dimension of the Lindblad equation is too high to allow for such a deterministic approximation, we combine classical ensemble averages of quantum Monte Carlo trajectories and a denoising technique. Specifically, a variance reduction method based upon the consideration of a low-rank dynamics as a control variate is developed. Numerical tests for quantum collapse and revivals show the efficiency of each approach, along with the complementarity of the two approaches.


title: Energy dissipation in magnetic null points at kinetic scales

abstract: We use kinetic particle-in-cell and magnetohydrodynamic simulations supported by an observational dataset to investigate magnetic reconnection in clusters of null points in space plasma. The magnetic configuration under investigation is driven by fast adiabatic flux rope compression that dissipates almost half of the initial magnetic field energy. In this phase powerful currents are excited producing secondary instabilities, and the system is brought into a state of `intermittent turbulence' within a few ion gyro-periods. Reconnection events are distributed all over the simulation domain and energy dissipation is rather volume-filling. Numerous spiral null points interconnected via their spines form null lines embedded into magnetic flux ropes; null point pairs demonstrate the signatures of torsional spine reconnection. However, energy dissipation mainly happens in the shear layers formed by adjacent flux ropes with oppositely directed currents. In these regions radial null pairs are spontaneously emerging and vanishing, associated with electron streams and small-scale current sheets. The number of spiral nulls in the simulation outweighs the number of radial nulls by a factor of 5\---10, in accordance with Cluster observations in the Earth's magnetosheath. Twisted magnetic fields with embedded spiral null points might indicate the regions of major energy dissipation for future space missions such as Magnetospheric Multiscale Mission (MMS).

categories: astro-ph.EP astro-ph.SR physics.plasm-ph physics.space-ph


title: Role of Z-pinches in magnetic reconnection in space plasmas

abstract: A widely accepted scenario of magnetic reconnection in collisionless space plasmas is the breakage of magnetic field lines in X-points. In laboratory, reconnection is commonly studied in pinches, current channels embedded into twisted magnetic fields. No model of magnetic reconnection in space plasmas considers both null-points and pinches as peers. We have performed a particle-in-cell simulation of magnetic reconnection in a three-dimensional configuration where null-points are present initially, and Z-pinches are formed during the simulation along the lines of spiral null-points. The non-spiral null-points are more stable than spiral ones, and no substantial energy dissipation is associated with them. On the contrary, turbulent magnetic reconnection in the pinches causes the magnetic energy to decay at a rate of ~1.5% per ion gyro period. Dissipation in similar structures is a likely scenario in space plasmas with large fraction of spiral null-points.

categories: physics.space-ph astro-ph.EP astro-ph.SR physics.plasm-ph


title: Load Shifting in the Smart Grid: To Participate or Not?

abstract: Demand-side management (DSM) has emerged as an important smart grid feature that allows utility companies to maintain desirable grid loads. However, the success of DSM is contingent on active customer participation. Indeed, most existing DSM studies are based on game-theoretic models that assume customers will act rationally and will voluntarily participate in DSM. In contrast, in this paper, the impact of customers' subjective behavior on each other's DSM decisions is explicitly accounted for. In particular, a noncooperative game is formulated between grid customers in which each customer can decide on whether to participate in DSM or not. In this game, customers seek to minimize a cost function that reflects their total payment for electricity. Unlike classical game-theoretic DSM studies which assume that customers are rational in their decision-making, a novel approach is proposed, based on the framework of prospect theory (PT), to explicitly incorporate the impact of customer behavior on DSM decisions. To solve the proposed game under both conventional game theory and PT, a new algorithm based on fictitious player is proposed using which the game will reach an epsilon-mixed Nash equilibrium. Simulation results assess the impact of customer behavior on demand-side management. In particular, the overall participation level and grid load can depend significantly on the rationality level of the players and their risk aversion tendency.


title: Lower bound on the blow-up rate of the 3D Navier-Stokes equations in H^{5/2}

abstract: Under assumption that $T^{\ast}$ is the maximal time of existence of smooth solution of the 3D Navier-Stokes equations in the Sobolev space $H^{s}$, we establish lower bounds for the blow-up rate of the type$\ \left( T^{\ast }-t\right) ^{-\varphi\left( n\right) }$, where $n$ is a natural number independent of $s$ and $\varphi$ is a linear function. Using this new type in the 3D Navier-Stokes equations in the $H^{5/2}$, both on the whole space and in the periodic case, we give an answer to a question left open by James et al (2012, J. Math. Phys.). We also prove optimal lower bounds for the blow-up rate in $\dot{H}^{3/2}$ and in $H^{1}$.


title: Dynamics of a planar thin shell at a Taub-FRW junction

abstract: We address the problem of stitching together the vacuum, static, planar-symmetric Taub spacetime and the flat Friedmann-Robertson-Walker cosmology using the Israel thin-shell formalism. The joining of Taub and FRW spacetimes is reminiscent of the Oppenheimer-Snyder collapse used in modeling the formation of a singularity from a collapsing spherical ball of dust. A possible mechanism for the formation of a planar singularity is provided. It is hoped that tackling such example will improve our intuition on planar-symmetric systems in Einstein's general relativity.


title: An Efficient Local Strategy to Control Information Spreading in Network

abstract: In social networks, control of rumor spread is an active area of research. SIR model is generally used to study the rumor dynamics in network while considering the rumor as an epidemic. In disease spreading model, epidemic is controlled by removing central nodes in the network. Full network information is needed for such removal. To have the information of complete network is difficult proposition. As a consequence, the search of an algorithm that may control epidemic without needing global information is a matter of great interest. In this paper, an immunization strategy is proposed that uses only local information available at a node, viz. degree of the node and average degree of its neighbour nodes. Proposed algorithm has been evaluated for scale-free network using SIR model. Numerical results show that proposed method has less complexity and gives significantly better results in comparison with other strategies while using only local information.


title: Logical operation of one-dimensional photonic crystal based on series and parallel connection

abstract: In this paper, we have proposed the compound structure of one-dimensional photonic crystal (PC), which includes series connection and parallel connection PC. We have studied the transmission characteristics of series connection and parallel connection PC, and obtained some new results. In addition, we have proved the series connection one-dimensional PC can realize the logical AND operation, and the parallel connection one-dimensional PC can realize the logical OR operation. The compound structure of one-dimensional PC can design more new type structure optical devices, and will provide the basic for designing quantum computer.


title: Discrete-Valued Control by Sum-of-Absolute-Values Optimization

abstract: In this paper, we propose a new design method of discrete-valued control for continuous-time linear time-invariant systems based on sum-of-absolute-values (SOAV) optimization. We first formulate the discrete-valued control design as a finite-horizon SOAV optimal control, which is an extended version of L1 optimal control. We then give simple conditions that guarantee the existence, discreteness, and uniqueness of the SOAV optimal control. Also, we give the continuity property of the value function, by which we prove the stability of infinite-horizon model predictive SOAV control systems. We provide a fast algorithm for the SOAV optimization based on the alternating direction method of multipliers (ADMM), which has an important advantage in real-time control computation. A simulation result shows the effectiveness of the proposed method.


title: Energetics of kinetic reconnection in a three-dimensional null points cluster

abstract: We performed three-dimensional Particle-in-Cell simulations of magnetic reconnection with multiple magnetic null points. Magnetic field energy conversion into kinetic energy was about five times higher than in traditional Harris sheet configuration. More than 85% of initial magnetic field energy was transferred to particle energy during 25 reversed ion cyclofrequencies. Magnetic reconnection in the cluster of null points evolved in three phases. During the first phase, ion beams were excited, that then gave part of their energy back to magnetic field in the second phase. In the third phase, magnetic reconnection occurs in many small patches around the current channels formed along the stripes of low magnetic field. Magnetic reconnection in null points presents essentially three-dimensional features, with no two dimensional symmetries or current sheets.