title: Dynamics of decoherence: universal scaling of the decoherence factor

abstract: We study the time dependence of the decoherence factor (DF) of a qubit globally coupled to an environmental spin system (ESS) which is driven across the quantum critical point (QCP) by varying a parameter of its Hamiltonian in time $t$ as $1 -t/\tau$ or $-t/\tau$, to which the qubit is coupled starting at the time $t \to -\infty$; here, $\tau$ denotes the inverse quenching rate. In the limit of weak coupling, we analyze the time evolution of the DF in the vicinity of the QCP (chosen to be at $t=0$) and define three quantities, namely, the generalized fidelity susceptibility $\chi_F(\tau)$ (defined right at the QCP), and the decay constants $\alpha_1 (\tau)$ and $\alpha_2 (\tau)$ which dictate the decay of the DF at a small but finite $t$($>0$). Using a dimensional analysis argument based on the Kibble-Zurek healing length, we show that $\chi_F(\tau)$ as well as $\alpha_1 (\tau)$ and $\alpha_2(\tau)$ indeed satisfy universal power-law scaling relations with $\tau$ and the exponents are solely determined by the spatial dimensionality of the ESS and the exponents associated with its QCP. Remarkably, using the numerical t-DMRG method, these scaling relations are shown to be valid in both the situations when the ESS is integrable and non-integrable and also for both linear and non-linear variation of the parameter. Furthermore, when an integrable ESS is quenched far away from the QCP, there is a predominant Gaussian decay of the DF with a decay constant which also satisfies a universal scaling relation.


title: Large Communities in a scale-free network

abstract: We prove the existence of a large complete subgraph w.h.p. in a preferential attachment random graph process with an edge-step. That is, we prove that the random graph $G_{t}$ produced by the so-called GLP model at time $t$ contains a complete subgraph of order $t^\alpha$, where $\alpha = (1-\varepsilon)\frac{1-p}{2-p}$, $\varepsilon$ is any number such that $0<\varepsilon<1$, and $0<p<1$ is a parameter of the model.


title: Optoelectronic switch based on intrinsic dual Schottky diodes in ambipolar MoSe$_2$ field-effect transistors

abstract: Here, we report the observation of a hitherto unreported optoelectronic effect, namely a light-induced diode-like response in multi-layered MoSe$_2$ field-effect transistors whose sense of current rectification is controllable through a gate voltage. We argue, through numerical simulations, that this behavior results from the difference in the size of the Schottky barriers between drain and source metal contacts. Each barrier can be modeled as a Schottky diode but with opposite senses of current rectification between them, with the diode response resulting from the light induced promotion of photo-generated carriers across the smaller barrier. The back gate voltage controls the sense of current rectification by modulating the relative amplitude between them. This effect, which gives rise to a novel type of optoelectronic switch, also yields a photovoltaic response. Hence, it could provide an alternative to PN-junctions when harvesting photovoltaic currents from transition metal dichalcogenides. We argue that the photovoltaic efficiency associated to this effect could be increased by just increasing the relative asymmetry between both Schottky barriers. We also suggest that this new electro-optical effect has potential for technological applications.


title: Wei-Norman-Kolokolov approach for Landau-Zener problems

abstract: The tunneling between the $2S+1$ Zeeman multiplets of an arbitrary spin $S$ using the $SU(2)$ Lie group exponential ordering Wei-Norman technique is investigated. A spin subjected to a regular time-dependent magnetic field and a transverse colored noise is considered and the regimes of fast and slow noise examined.


title: Non-linear Fano Interferences in Open Quantum Systems: an Exactly Solvable Model

abstract: We obtain an explicit solution for the stationary state populations of a dissipative Fano model, where a discrete excited state is coupled to a continumm set of states; both excited set of states are reachable by photo-excitation from the ground state. The dissipative dynamic is described by a Liouville equation in Lindblad form and the field intensity can take arbitrary values within the model. We show that the continuum states population as a function of laser frequency can always be expressed as a Fano profile plus a Lorentzian function with effective parameters whose explicit expressions are given in the case of a closed system coupled to a bath as well as for the original Fano scattering framework. Although the solution is intricate, it can be elegantly expressed as a linear transformation of the kernel of a $4\times 4$ matrix which has the meaning of an effective Liouvillian. We unveil key notable processes related to the optical nonlinearity and which had not been reported to date: electromagnetic induced-transparency, population inversions, power narrowing and broadening, as well as an effective reduction of the Fano asymmetry parameter.


title: Self-dual configurations in Abelian Higgs models with $k$-generalized gauge field dynamics

abstract: We have shown the existence of self-dual solutions in new Maxwell-Higgs scenarios where the gauge field possesses a $k$-generalized dynamic, i.e., the kinetic term of gauge field is a highly nonlinear function of $F_{\mu\nu}F^{\mu\nu}$. We have implemented our proposal by means of a $k$-generalized model displaying the spontaneous symmetry breaking phenomenon. We implement consistently the Bogomol'nyi-Prasad-Sommerfield formalism providing highly nonlinear self-dual equations whose solutions are electrically neutral possessing total energy proportional to the magnetic flux. Among the infinite set of possible configurations, we have found families of $k$-generalized models whose self-dual equations have a form mathematically similar to the ones arising in the Maxwell-Higgs or Chern-Simons-Higgs models. Furthermore, we have verified that our proposal also supports infinite twinlike models with $|\phi|^4$-potential or $|\phi |^6$-potential. With the aim to show explicitly that the BPS equations are able to provide well-behaved configurations, we have considered a test model in order to study axially symmetric vortices. By depending of the self-dual potential, we have shown that the $k$-generalized model is able to produce solutions that for long distances have a exponential decay (as Abrikosov-Nielsen-Olesen vortices) or have a power-law decay (characterizing delocalized vortices). In all cases, we observe that the generalization modifies the vortex core size, the magnetic field amplitude and the bosonic masses but the total energy remains proportional to the quantized magnetic flux.


title: Frequency shifts and relaxation rates for spin 1/2 particles moving in electromagnetic fields

abstract: We discuss the behaviour of the Larmor frequency shift and the longitudinal relaxation rate due to non-uniform electromagnetic fields on an assembly of spin 1/2 particles, in adiabatic and nonadiabatic regimes. We also show some general relations between the various frequency shifts and between the frequency shifts and relaxation rates. The remarkable feature of all our results is that they were obtained without any specific assumptions on the explicit form of the correlation functions of the fields. Hence, we expect that our results are valid both for diffusive and ballistic regime of motion and arbitrary cell shapes and surface scattering. These results can then be applied to a wide variety of realistic systems.


title: Simple model of the slingshot effect

abstract: We present a detailed quantitative description of the recently proposed "slingshot effect" [Fiore, Fedele, De angelis 2014]. Namely, we determine a broad range of conditions under which the impact of a very short and intense laser pulse normally onto a low-density plasma (or matter to be locally completely ionized into a plasma by the pulse) causes the expulsion of a bunch of surface electrons in the direction opposite to the one of propagation of the pulse, and the detailed, ready-for-experiments features of the expelled electrons (energy spectrum, collimation, etc). The effect is due to the combined actions of the ponderomotive force and the huge longitudinal field arising from charge separation. Our predictions are based on estimating 3D corrections to a simple, yet powerful plane magnetohydrodynamic (MHD) model where the equations to be solved are reduced to a system of Hamilton equations in one dimension (or a collection of) which become autonomous after the pulse has overcome the electrons. Experimental tests seem to be at hand. If confirmed by the latter, the effect would provide a new extraction and acceleration mechanism for electrons, alternative to traditional radio-frequency-based or Laser-Wake-Field ones.


title: Centrality-dependent modification of jet-production rates in deuteron-gold collisions at $\sqrt{s_{NN}}$=200 GeV

abstract: Jet production rates are measured in $p$$+$$p$ and $d$$+$Au collisions at $\sqrt{s_{NN}}$=200 GeV recorded in 2008 with the PHENIX detector at the Relativistic Heavy Ion Collider. Jets are reconstructed using the $R=0.3$ anti-$k_{t}$ algorithm from energy deposits in the electromagnetic calorimeter and charged tracks in multi-wire proportional chambers, and the jet transverse momentum ($p_T$) spectra are corrected for the detector response. Spectra are reported for jets with $12<p_T<50$ GeV/$c$, within a pseudorapidity acceptance of $\left|\eta\right|<0.3$. The nuclear-modification factor ($R_{d{\rm Au}}$) values for 0\%--100\% $d$$+$Au events are found to be consistent with unity, constraining the role of initial state effects on jet production. However, the centrality-selected $R_{d{\rm Au}}$ values and central-to-peripheral ratios ($R_{\rm CP}$) show large, $p_T$-dependent deviations from unity, which challenge the conventional models that relate hard-process rates and soft-particle production in collisions involving nuclei.


title: WSPEC: A waveguide filter-bank focal plane array spectrometer for millimeter wave astronomy and cosmology

abstract: Imaging and spectroscopy at (sub-)millimeter wavelengths are key frontiers in astronomy and cosmology. Large area spectral surveys with moderate spectral resolution (R=50-200) will be used to characterize large scale structure and star formation through intensity mapping surveys in emission lines such as the CO rotational transitions. Such surveys will also be used to study the SZ effect, and will detect the emission lines and continuum spectrum of individual objects. WSPEC is an instrument proposed to target these science goals. It is a channelizing spectrometer realized in rectangular waveguide, fabricated using conventional high-precision metal machining. Each spectrometer is coupled to free space with a machined feed horn, and the devices are tiled into a 2D array to fill the focal plane of the telescope. The detectors will be aluminum Lumped-Element Kinetic Inductance Detectors (LEKIDs). To target the CO lines and SZ effect, we will have bands at 135-175 GHz and 190-250 GHz, each Nyquist-sampled at R~200 resolution. Here we discuss the instrument concept and design, and successful initial testing of a WR10 (i.e. 90 GHz) prototype spectrometer. We recently tested a WR5 (180 GHz) prototype to verify that the concept works at higher frequencies, and also designed a resonant backshort structure that may further increase the optical efficiency. We are making progress towards integrating a spectrometer with a LEKID array and deploying a prototype device to a telescope for first light.


title: Top Local Cohomology Modules Are Almost Always Non-Artinian

abstract: Let $(R,\mathfrak m)$ be a Noetherian local ring, $I$ an ideal of $R$ and $M$ a weakly finite or a coatomic $R$-module of dimension $n$. In this article, we resolve the Artinianness and non-Artinianness of top local cohomology modules, $H^{cd(I,M)}_I(M)$, in all cases except in the case $cd(I,M)=n-1$ and $dim(R/{I+Ann(M)})>1$ for which we have some sorter results under certain conditions.


title: Solar Wind at 33 AU: Setting Bounds on the Pluto Interaction for New Horizons

abstract: The NASA New Horizons spacecraft flies past Pluto on July 14, 2015, carrying two instruments that detect charged particles. Pluto has a tenuous, extended atmosphere that is escaping the weak gravity of the planet. The interaction of the solar wind with the escaping atmosphere of Pluto depends on solar wind conditions as well as the vertical structure of the atmosphere. We have analyzed Voyager 2 particles and fields measurements between 25 and 39 AU and present their statistical variations. We have adjusted these predictions to allow for the declining activity of the Sun and solar wind output. We summarize the range of SW conditions that can be expected at 33 AU and survey the range of scales of interaction that New Horizons might experience. Model estimates for the solar wind stand-off distance vary from approximately 7 to 1000 RP with our best estimate being around 40 RP (where we take the radius of Pluto to be RP=1184 km).


title: Reheating of the Universe as holographic thermalization

abstract: Assuming gauge/gravity correspondence we study reheating of the Universe using its holographic dual. Inflaton decay and thermalisation of the decay products correspond to collapse of a spherical shell and formation of a blackhole in the dual anti-de Sitter (AdS) spacetime. The reheating temperature is computed as the Hawking temperature of the developed blackhole probed by a dynamical boundary, and is determined by the inflaton energy density and the AdS radius, with corrections from the dynamics of the shell collapse. For given initial energy density of the inflaton field the holographic model gives significantly lower reheating temperature than the instant reheating scenario, while it is shown to be safely within phenomenological bounds.


title: Single electron yields from semileptonic charm and bottom hadron decays in Au$+$Au collisions at $\sqrt{s_{NN}}=200$ GeV

abstract: The PHENIX Collaboration at the Relativistic Heavy Ion Collider has measured open heavy-flavor production in minimum bias Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV via the yields of electrons from semileptonic decays of charm and bottom hadrons. Previous heavy-flavor electron measurements indicated substantial modification in the momentum distribution of the parent heavy quarks due to the quark-gluon plasma created in these collisions. For the first time, using the PHENIX silicon vertex detector to measure precision displaced tracking, the relative contributions from charm and bottom hadrons to these electrons as a function of transverse momentum are measured in Au$+$Au collisions. We compare the fraction of electrons from bottom hadrons to previously published results extracted from electron-hadron correlations in $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV and find the fractions to be similar within the large uncertainties on both measurements for $p_T>4$ GeV/$c$. We use the bottom electron fractions in Au$+$Au and $p$$+$$p$ along with the previously measured heavy flavor electron $R_{AA}$ to calculate the $R_{AA}$ for electrons from charm and bottom hadron decays separately. We find that electrons from bottom hadron decays are less suppressed than those from charm for the region $3<p_T<4$ GeV/$c$.


title: Supernova Driving. I. The Origin of Molecular Cloud Turbulence

abstract: Turbulence is ubiquitous in molecular clouds (MCs), but its origin is still unclear because MCs are usually assumed to live longer than the turbulence dissipation time. Interstellar medium (ISM) turbulence is likely driven by SN explosions, but it has never been demonstrated that SN explosions can establish and maintain a turbulent cascade inside MCs consistent with the observations. In this work, we carry out a simulation of SN-driven turbulence in a volume of (250 pc)$^3$, specifically designed to test if SN driving alone can be responsible for the observed turbulence inside MCs. We find that SN driving establishes a velocity scaling consistent with the usual scaling laws of supersonic turbulence, suggesting that previous idealized simulations of MC turbulence, driven with a random, large-scale volume force, were correctly adopted as appropriate models for MC turbulence, despite the artificial driving. We also find that the same scaling laws extend to the interior of MCs, and that the velocity-size relation of the MCs selected from our simulation is consistent with that of MCs from the Outer-Galaxy Survey, the largest MC sample available. The mass-size relation and the mass and size probability distributions also compare successfully with those of the Outer Galaxy Survey. Finally, we show that MC turbulence is super-Alfv\'{e}nic with respect to both the mean and rms magnetic-field strength. We conclude that MC structure and dynamics are the natural result of SN-driven turbulence.


title: Self-Configuring and Evolving Fuzzy Image Thresholding

abstract: Every segmentation algorithm has parameters that need to be adjusted in order to achieve good results. Evolving fuzzy systems for adjustment of segmentation parameters have been proposed recently (Evolving fuzzy image segmentation -- EFIS [1]. However, similar to any other algorithm, EFIS too suffers from a few limitations when used in practice. As a major drawback, EFIS depends on detection of the object of interest for feature calculation, a task that is highly application-dependent. In this paper, a new version of EFIS is proposed to overcome these limitations. The new EFIS, called self-configuring EFIS (SC-EFIS), uses available training data to auto-configure the parameters that are fixed in EFIS. As well, the proposed SC-EFIS relies on a feature selection process that does not require the detection of a region of interest (ROI).


title: Rotating cylindrical wormholes and energy conditions

abstract: We seek wormholes among rotating cylindrically symmetric configurations in general relativity. Exact wormhole solutions are presented with such sources of gravity as a massless scalar field, a cosmological constant, and a scalar field with an exponential potential. However, none of these solutions are asymptotically flat, which excludes the existence of wormhole entrances as local objects in our Universe. To overcome this difficulty, we try to build configurations with flat asymptotic regions using the cut-and-paste procedure: on both sides of the throat, a wormhole solution is matched to a properly chosen region of flat space-time at some surfaces $\Sigma_-$ and $\Sigma_+$. It is shown, however, that if the source of gravity in the throat region is a scalar field with an arbitrary potential, then one or both thin shells appearing on $\Sigma_-$ and $\Sigma_+$ inevitably violate the null energy condition. Thus, although rotating wormhole solutions are easily found without exotic matter, such matter is still necessary for obtaining asymptotic flatness.


title: Overlap distributions for quantum quenches in the anisotropic Heisenberg chain

abstract: The dynamics after a quantum quench is determined by the weights of the initial state in the eigenspectrum of the final Hamiltonian, i.e., by the distribution of overlaps in the energy spectrum. We present an analysis of such overlap distributions for quenches of the anisotropy parameter in the one-dimensional anisotropic spin-1/2 Heisenberg model (XXZ chain). We provide an overview of the form of the overlap distribution for quenches from various initial anisotropies to various final ones, using numerical exact diagonalization. We show that if the system is prepared in the antiferromagnetic N\'eel state (infinite anisotropy) and released into a non-interacting setup (zero anisotropy, XX point) only a small fraction of the final eigenstates gives contributions to the post-quench dynamics, and that these eigenstates have identical overlap magnitudes. We derive expressions for the overlaps, and present the selection rules that determine the final eigenstates having nonzero overlap. We use these results to derive concise expressions for time-dependent quantities (Loschmidt echo, longitudinal and transverse correlators) after the quench. We use perturbative analyses to understand the overlap distribution for quenches from infinite to small nonzero anisotropies, and for quenches from large to zero anisotropy.


title: Dielectron production in Au$+$Au collisions at $\sqrt{s_{NN}}$=200 GeV

abstract: We present measurements of $e^+e^-$ production at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of mass ($m_{ee} <$ 5 GeV/$c^2$) and pair transverse momentum ($p_T$ $<$ 5 GeV/$c$), for minimum bias and for five centrality classes. The \ee yield is compared to the expectations from known sources. In the low-mass region ($m_{ee}=0.30$--0.76 GeV/$c^2$) there is an enhancement that increases with centrality and is distributed over the entire pair \pt range measured. It is significantly smaller than previously reported by the PHENIX experiment and amounts to $2.3\pm0.4({\rm stat})\pm0.4({\rm syst})\pm0.2^{\rm model}$ or to $1.7\pm0.3({\rm stat})\pm0.3({\rm syst})\pm0.2^{\rm model}$ for minimum bias collisions when the open-heavy-flavor contribution is calculated with {\sc pythia} or {\sc mc@nlo}, respectively. The inclusive mass and $p_T$ distributions as well as the centrality dependence are well reproduced by model calculations where the enhancement mainly originates from the melting of the $\rho$ meson resonance as the system approaches chiral symmetry restoration. In the intermediate-mass region ($m_{ee}$ = 1.2--2.8 GeV/$c^2$), the data hint at a significant contribution in addition to the yield from the semileptonic decays of heavy-flavor mesons.


title: Non-Abelian cosmic strings in de Sitter and anti-de Sitter space

abstract: In this paper we investigate the non-Abelian cosmic string in de Sitter and anti-de Sitter spacetimes. In order to do that we construct the complete set of equations of motion considering the presence of a cosmological constant. By using numerical analysis we provide the behavior of the Higgs and gauge fields and also for the metric tensor for specific values of the physical parameters of the theory. For de Sitter case, we find the appearance of horizons that although being consequence of the presence of the cosmological constant it strongly depends on the value of the gravitational coupling. In the anti-de Sitter case, we find that the system does not present horizons. In fact the new feature of this system is related with the behavior of the $(00)$ and $(zz)$ components of the metric tensor. They present a strongly increasing for large distance from the string.


title: Warped cones over profinite completions

abstract: We construct metric spaces that do not have property A yet are coarsely embeddable into the Hilbert space. Our examples are so called warped cones, which were introduced by J. Roe to serve as examples of spaces non-embeddable into a Hilbert space and with or without property A. The construction provides the first examples of warped cones combining coarse embeddability and lack of property A. We also construct warped cones over manifolds with isometrically embedded expanders and generalise Roe's criteria for the lack of property A or coarse embeddability of a warped cone. Along the way, it is proven that property A of the warped cone over a profinite completion is equivalent to amenability of the group. In the appendix we solve a problem of Nowak regarding his examples of spaces with similar properties.


title: Searches for Higgs boson pair production in the $hh\to bb\tau\tau, \gamma\gamma WW*, \gamma\gamma bb, bbbb$ channels with the ATLAS detector

abstract: Searches for both resonant and non-resonant Higgs boson pair production are performed in the $hh\to bb\tau\tau, \gamma\gamma WW^*$ final states using 20.3 fb$^{-1}$ of $pp$ collision data at a center-of-mass energy of 8 TeV recorded with the ATLAS detector at the Large Hadron Collider. No evidence of their production is observed and 95\% confidence level upper limits on the production cross sections are set. These results are then combined with the published results of the $hh\to \gamma\gamma bb, bbbb$ analyses. An upper limit of 0.69 (0.47) pb on the non-resonant Standard Model like $hh$ production is observed (expected), corresponding to 70 (48) times of the SM $gg\to hh$ cross section. For production via narrow resonances, cross section limits of $hh$ production from a heavy Higgs boson decay are set as a function of the heavy Higgs boson mass. The observed (expected) limits range from 2.1 (1.1) pb at 260 GeV to 0.011 (0.018) pb at 1000 GeV. These results are interpreted in the context of two simplified scenarios of the Minimal Supersymmetric Standard Model.


title: Unravelling Medium Effects in Heavy Ion Collisions with Zeal

abstract: We propose a new observable, called zeal, to analyze events with jets in heavy ion collisions. The observable measures how a thermal medium affects the multiplicity and distribution of energetic particles in a jet. Using few known models for energy loss and jet quenching, we demonstrate its capability to distinguish the physics of these models.


title: Lambda Boo Abundance Patterns: Accretion from Orbiting Sources

abstract: The abundance anomalies in lambda Boo stars are popularly explained by element-specific mass inflows at rates that are much greater than empirically-inferred bounds for interstellar accretion. Therefore, a lambda Boo star's thin outer envelope must derive from a companion star, planet, analogs to Kuiper Belt Objects or a circumstellar disk. Because radiation pressure on gas-phase ions might selectively allow the accretion of carbon, nitrogen, and oxygen and inhibit the inflow of elements such as iron, the source of the acquired matter need not contain dust. We propose that at least some lambda Boo stars accrete from the winds of hot Jupiters.


title: Comment on "On Phase Selective Quantum Eraser" (arXiv:1501.00817 [quant-ph])

abstract: In a recent interesting article A. Heuer, G. Pieplow, and R. Menzel discuss a quantum-eraser experiment to investigate the complementarity of wave-like and particle-like behavior of photons. I would like to draw your attention to the fact that the very same experimental setup was suggested in a previous paper, and take advantage of this opportunity to examine some aspects of this controversial subject.


title: Impact of Transceiver Impairments on the Capacity of Dual-Hop Relay Massive MIMO Systems

abstract: Despite the deleterious effect of hardware impairments on communication systems, most prior works have not investigated their impact on widely used relay systems. Most importantly, the application of inexpensive transceivers, being prone to hardware impairments, is the most cost-efficient way for the implementation of massive multiple-input multiple-output (MIMO) systems. Consequently, the direction of this paper is towards the investigation of the impact of hardware impairments on MIMO relay networks with large number of antennas. Specifically, we obtain the general expression for the ergodic capacity of dual-hop (DH) amplify-and-forward (AF) relay systems. Next, given the advantages of the free probability (FP) theory with comparison to other known techniques in the area of large random matrix theory, we pursue a large limit analysis in terms of number of antennas and users by shedding light to the behavior of relay systems inflicted by hardware impairments.


title: No Keplerian Disk >10 AU around the Protostar B335: Magnetic Braking or Young Age?

abstract: We have conducted ALMA cycle 2 observations in the 1.3 mm continuum and in the C18O (2-1) and SO (5_6-4_5) lines at a resolution of ~0.3" toward the Class 0 protostar B335. The 1.3 mm continuum, C18O, and SO emission all show central compact components with sizes of ~40-180 AU within more extended components. The C18O component shows signs of infalling and rotational motion. By fitting simple kinematic models to the C18O data, the protostellar mass is estimated to be 0.05 Msun. The specific angular momentum, on a 100 AU scale, is ~4.3E-5 km/s*pc. A similar specific angular momentum, ~3E-5 to 5E-5 km/s*pc, is measured on a 10 AU scale from the velocity gradient observed in the central SO component, and there is no clear sign of an infalling motion in the SO emission. By comparing the infalling and rotational motion, our ALMA results suggest that the observed rotational motion has not yet reached Keplerian velocity neither on a 100 AU nor even on a 10 AU scale. Consequently, the radius of the Keplerian disk in B335 (if present) is expected to be 1-3 AU. The expected disk radius in B335 is one to two orders of magnitude smaller than those of observed Keplerian disks around other Class 0 protostars. Based on the observed infalling and rotational motion from 0.1 pc to inner 100 AU scales, there are two possible scenarios to explain the presence of such a small Keplerian disk in B335: magnetic braking and young age. If our finding is the consequence of magnetic braking, ~50% of the angular momentum of the infalling material within a 1000 AU scale might have been removed, and the magnetic field strength on a 1000 AU scale is estimated to be ~200 uG. If it is young age, the infalling radius in B335 is estimated to be ~2700 AU, corresponding to a collapsing time scale of ~5E4 yr.


title: Band spectrum is D-brane

abstract: We show that band spectrum of topological insulators can be identified as the shape of D-branes in string theory. The identification is based on a relation between the Berry connection associated with the band structure and the ADHM/Nahm construction of solitons whose geometric realization is available with D-branes. We also show that chiral and helical edge states are identified as D-branes representing a noncommutative monopole.


title: Note: "Lock-in accelerometry" to follow sink dynamics in shaken granular matter

abstract: Understanding the penetration dynamics of intruders in granular beds is relevant not only for fundamental Physics, but also for geophysical processes and construction on sediments or granular soils in areas potentially affected by earthquakes. While the penetration of intruders in two dimensional (2D) laboratory granular beds can be followed using video recording, it is useless in three dimensional (3D) beds of non-transparent materials such as common sand. Here we propose a method to quantify the sink dynamics of an intruder into laterally shaken granular beds based on the temporal correlations between the signals from a reference accelerometer fixed to the shaken granular bed, and a probe accelerometer deployed inside the intruder. Due to its analogy with the working principle of a lock in amplifier, we call this technique Lock in accelerometry (LIA). During Earthquakes, some soils can lose their ability to sustain shear and deform, causing subsidence and sometimes substantial building damage due to deformation or tumbling


title: Application of bond valence method in the non-isovalent semiconductor alloy (GaN)$_{1-x}$(ZnO)$_x$

abstract: This paper studies the bond valence method (BVM) and its application in the non-isovalent semiconductor alloy (GaN)$_{\rm{1-x}}$(ZnO)$_{\rm{x}}$. Particular attention is paid to the role of short-range order (SRO). A physical interpretation based on atomic orbital interaction is proposed and examined by density-functional theory (DFT) calculations. Combining BVM with Monte-Carlo simulations and a DFT-based cluster expansion model, bond-length distributions and bond-angle variations are predicted. The correlation between bond valence and bond stiffness is also revealed. Finally the concept of bond valence is extended into the modelling of an atomistic potential.


title: $1$-cohomology of simplicial amalgams of groups

abstract: We develop a cohomological method to classify amalgams of groups. We generalize this to simplicial amalgams in any concrete category. We compute the non-commutative 1-cohomology for several examples of amalgams defined over small simplices.


title: Analytic approximation of transmutation operators and related systems of functions

abstract: In arXiv:1306.2914 a method for approximate solution of Sturm-Liouville equations and related spectral problems was presented based on the construction of the Delsarte transmutation operators. The problem of numerical approximation of solutions was reduced to approximation of a primitive of the potential by a finite linear combination of certain specially constructed functions obtained from the generalized wave polynomials introduced in arXiv:1208.5984 and arXiv:1208.6166. The method allows one to compute both lower and higher eigendata with an extreme accuracy. Since the solution of the approximation problem is the main step in the application of the method, the properties of the system of functions involved are of primary interest. In arXiv:1306.2914 two basic properties were established: the completeness in appropriate functional spaces and the linear independence. In this paper we present a considerably more complete study of the systems of functions. We establish their relation with another linear differential second-order equation, find out certain operations (in a sense, generalized derivatives and antiderivatives) which allow us to generate the next such function from a previous one. We obtain the uniqueness of the coefficients of expansions in terms of such functions and a corresponding generalized Taylor theorem. We construct the invertible integral operators transforming powers of the independent variable into the functions under consideration and establish their commutation relations with differential operators. We present some error bounds for the solution of the approximation problem depending on the smoothness of the potential and show that these error bounds are close to optimal in order. Also, we provide a rigorous justification of the alternative formulation of the proposed method allowing one to make use of the known initial values of the solutions at an endpoint.


title: Large-Scale Optimization Algorithms for Sparse Conditional Gaussian Graphical Models

abstract: This paper addresses the problem of scalable optimization for L1-regularized conditional Gaussian graphical models. Conditional Gaussian graphical models generalize the well-known Gaussian graphical models to conditional distributions to model the output network influenced by conditioning input variables. While highly scalable optimization methods exist for sparse Gaussian graphical model estimation, state-of-the-art methods for conditional Gaussian graphical models are not efficient enough and more importantly, fail due to memory constraints for very large problems. In this paper, we propose a new optimization procedure based on a Newton method that efficiently iterates over two sub-problems, leading to drastic improvement in computation time compared to the previous methods. We then extend our method to scale to large problems under memory constraints, using block coordinate descent to limit memory usage while achieving fast convergence. Using synthetic and genomic data, we show that our methods can solve one million dimensional problems to high accuracy in a little over a day on a single machine.


title: Robust Sensitivity Analysis of the Optimal Value of Linear Programming

abstract: We propose a framework for sensitivity analysis of linear programs (LPs) in minimization form, allowing for simultaneous perturbations in the objective coefficients and right-hand sides, where the perturbations are modeled in a compact, convex uncertainty set. This framework unifies and extends multiple approaches for LP sensitivity analysis in the literature and has close ties to worst-case linear optimization and two-stage adaptive optimization. We define the minimum (best-case) and maximum (worst-case) LP optimal values, p- and p+, over the uncertainty set, and we discuss issues of finiteness, attainability, and computational complexity. While p- and p+ are difficult to compute in general, we prove that they equal the optimal values of two separate, but related, copositive programs. We then develop tight, tractable conic relaxations to provide lower and upper bounds on p- and p+, respectively. We also develop techniques to assess the quality of the bounds, and we validate our approach computationally on several examples from--and inspired by--the literature. We find that the bounds on p- and p+ are very strong in practice and, in particular, are at least as strong as known results for specific cases from the literature.


title: Effect of cross-coupling on the phase behavior in a biaxial nematic: Insights from Monte Carlo studies

abstract: Phase sequences of the biaxial nematic liquid crystal in the interior of the essential triangle are studied with Wang Landau sampling. The evidence points to the existence of an intermediate unixial phase with low biaxiality in the isotropic to biaxial nematic phase sequence.


title: Universal energy diffusion in a quivering billiard

abstract: We introduce and study a model of time-dependent billiard systems with billiard boundaries undergoing infinitesimal wiggling motions. The so-called quivering billiard is simple to simulate, straightforward to analyze, and is a faithful representation of time-dependent billiards in the limit of small boundary displacements. We assert that when a billiard's wall motion approaches the quivering motion, deterministic particle dynamics become inherently stochastic. Particle ensembles in a quivering billiard are shown to evolve to a universal energy distribution through an energy diffusion process, regardless of the billiard's shape or dimensionality, and as a consequence universally display Fermi acceleration. Our model resolves a known discrepancy between the one-dimensional Fermi-Ulam model and the simplified static wall approximation. We argue that the quivering limit is the true fixed wall limit of the Fermi-Ulam model.


title: ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method

abstract: Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using $N$-body simulations. COLA is a novel method designed to make this feasible by following an approximate dynamics but with up to 3 orders of magnitude speed-ups when compared to an exact $N$-body. In this paper we investigate the optimization of the code parameters in the compromise between computational cost and recovered accuracy in observables such as two-point clustering and halo abundance. We benchmark those observables with a state-of-the-art $N$-body run, the MICE Grand Challenge simulation (MICE-GC). We find that using 40 time steps linearly spaced since $z_i \sim 20$, and a force mesh resolution three times finer than that of the number of particles, yields a matter power spectrum within $1\%$ for $k \lesssim 1\,h {\rm Mpc}^{-1}$ and a halo mass function within $5\%$ of those in the $N$-body. In turn the halo bias is accurate within $2\%$ for $k \lesssim 0.7\,h {\rm Mpc}^{-1}$ whereas, in redshift space, the halo monopole and quadrupole are within $4\%$ for $k \lesssim 0.4\,h {\rm Mpc}^{-1}$. These results hold for a broad range in redshift ($0 < z < 1$) and for all halo mass bins investigated ($M > 10^{12.5} \, h^{-1} \, {\rm M_{\odot}}$). To bring accuracy in clustering to one percent level we study various methods that re-calibrate halo masses and/or velocities. We thus propose an optimized choice of COLA code parameters as a powerful tool to optimally exploit future galaxy surveys.


title: The fitness of the strongest individual in the subcritical GMS model

abstract: We derive the strongest individual fitness distribution on a variation for a species survival model proposed by Guiol, Machado and Schinazi \cite{GMS11}. We point out to the fact that this distribution relies on the Gauss hypergeometric function and when $p=\frac{1}{2}$ on the Hypergeometric function type I distribution


title: New expansion rate measurements of the Crab Nebula in radio and optical

abstract: We present new radio measurements of the expansion rate of the Crab nebula's synchrotron nebula over a ~30-yr period. We find a convergence date for the radio synchrotron nebula of CE 1255 +- 27. We also re-evaluated the expansion rate of the optical line emitting filaments, and we show that the traditional estimates of their convergence dates are slightly biased. Using an un-biased Bayesian analysis, we find a convergence date for the filaments of CE 1091 +- 34 (~40 yr earlier than previous estimates). Our results show that both the synchrotron nebula and the optical line-emitting filaments have been accelerated since the explosion in CE 1054, but that the synchrotron nebula has been relatively strongly accelerated, while the optical filaments have been only slightly accelerated. The finding that the synchrotron emission expands more rapidly than the filaments supports the picture that the latter are the result of the Rayleigh-Taylor instability at the interface between the pulsar-wind nebula and the surrounding freely-expanding supernova ejecta, and rules out models where the pulsar wind bubble is interacting directly with the pre-supernova wind of the Crab's progenitor.


title: Ultrafast terahertz spectroscopy study of Kondo insulating thin film SmB$_{6}$: evidence for an emergent surface state

abstract: We utilize terahertz time domain spectroscopy to investigate thin films of the heavy fermion compound SmB6, a prototype Kondo insulator. Temperature dependent terahertz (THz) conductivity measurements reveal a rapid decrease in the Drude weight and carrier scattering rate at ~T*=20 K, well below the hybridization gap onset temperature (100 K). Moreover, a low-temperature conductivity plateau (below 20K) indicates the emergence of a surface state with an effective electron mass of 0.1me. Conductivity dynamics following optical excitation are also measured and interpreted using Rothwarf-Taylor (R-T) phenomenology, yielding a hybridization gap energy of 17 meV. However, R-T modeling of the conductivity dynamics reveals a deviation from the expected thermally excited quasiparticle density at temperatures below 20K, indicative of another channel opening up in the low energy electrodynamics. Taken together, these results suggest the onset of a surface state well below the crossover temperature (100K) after long-range coherence of the f-electron Kondo lattice is established.


title: Random matrix techniques in quantum information theory

abstract: The purpose of this review article is to present some of the latest developments using random techniques, and in particular, random matrix techniques in quantum information theory. Our review is a blend of a rather exhaustive review, combined with more detailed examples -- coming from research projects in which the authors were involved. We focus on two main topics, random quantum states and random quantum channels. We present results related to entropic quantities, entanglement of typical states, entanglement thresholds, the output set of quantum channels, and violations of the minimum output entropy of random channels.


title: Properties of the Nearby Brown Dwarf WISEP J180026.60+013453.1

abstract: We present new spectroscopy and astrometry to characterize the nearby brown dwarf WISEP J180026.60+013453.1. The optical spectral type, L7.5, is in agreement with the previously reported near-infrared spectral type. The preliminary trigonometric parallax places it at a distance of $8.01 \pm 0.21$ pc, confirming that it is the fourth closest known late-L (L7-L9) dwarf. The measured luminosity, our detection of lithium, and the lack of low surface gravity indicators indicates that WISEP J180026.60+013453.1 has a mass $0.03 < M < 0.06 M_\odot$ and an age between 300 million and 1.5 billion years according to theoretical substellar evolution models. The low space motion is consistent with this young age. We have measured the rotational broadening ($v \sin i = 13.5 \pm 0.5$ km/s), and use it to estimate a maximum rotation period of 9.3 hr.


title: On the functoriality of Khovanov-Floer theories

abstract: We introduce the notion of a Khovanov-Floer theory. Roughly, such a theory assigns a filtered chain complex over Z/2 to a link diagram such that (1) the E_2 page of the resulting spectral sequence is naturally isomorphic to the Khovanov homology of the link; (2) this filtered complex behaves nicely under planar isotopy, disjoint union, and 1-handle addition; and (3) the spectral sequence collapses at the E_2 page for any diagram of the unlink. We prove that a Khovanov-Floer theory naturally yields a functor from the link cobordism category to the category of spectral sequences. In particular, every page (after E_1) of the spectral sequence accompanying a Khovanov-Floer theory is a link invariant, and an oriented cobordism in R^3 \times [0,1] between links in R^3 induces a map between each page of their spectral sequences, invariant up to smooth isotopy of the cobordism rel boundary. We then show that the spectral sequences relating Khovanov homology to Heegaard Floer homology and singular instanton knot homology are induced by Khovanov-Floer theories and are therefore functorial in the manner described above, as has been conjectured for some time. We further show that Szabo's geometric spectral sequence comes from a Khovanov-Floer theory, and is thus functorial as well. In addition, we illustrate how our framework can be used to give another proof that Lee's spectral sequence is functorial and that Rasmussen's invariant is a knot invariant. Finally, we use this machinery to define some potentially new knot invariants.


title: Scalable Metropolis Monte Carlo for simulation of hard shapes

abstract: We design and implement HPMC, a scalable hard particle Monte Carlo simulation toolkit, and release it open source as part of HOOMD-blue. HPMC runs in parallel on many CPUs and many GPUs using domain decomposition. We employ BVH trees instead of cell lists on the CPU for fast performance, especially with large particle size disparity, and optimize inner loops with SIMD vector intrinsics on the CPU. Our GPU kernel proposes many trial moves in parallel on a checkerboard and uses a block-level queue to redistribute work among threads and avoid divergence. HPMC supports a wide variety of shape classes, including spheres / disks, unions of spheres, convex polygons, convex spheropolygons, concave polygons, ellipsoids / ellipses, convex polyhedra, convex spheropolyhedra, spheres cut by planes, and concave polyhedra. NVT and NPT ensembles can be run in 2D or 3D triclinic boxes. Additional integration schemes permit Frenkel-Ladd free energy computations and implicit depletant simulations. In a benchmark system of a fluid of 4096 pentagons, HPMC performs 10 million sweeps in 10 minutes on 96 CPU cores on XSEDE Comet. The same simulation would take 7.6 hours in serial. HPMC also scales to large system sizes, and the same benchmark with 16.8 million particles runs in 1.4 hours on 2048 GPUs on OLCF Titan.


title: Well Control Optimization using Derivative-Free Algorithms and a Multiscale Approach

abstract: In this paper, we use numerical optimization algorithms and a multiscale approach in order to find an optimal well management strategy over the life of the reservoir. The large number of well rates for each control step make the optimization problem more difficult and at a high risk of achieving a suboptimal solution. Moreover, the optimal number of adjustments is not known a priori. Adjusting well controls too frequently will increase unnecessary well management and operation cost, and an excessively low number of control adjustments may not be enough to obtain a good yield. We investigate three derivative-free optimization algorithms, chosen for their robust and parallel nature, to determine optimal well control strategies. The algorithms chosen include generalized pattern search (GPS), particle swarm optimization (PSO) and covariance matrix adaptation evolution strategy (CMA-ES). These three algorithms encompass the breadth of available black-box optimization strategies: deterministic local search, stochastic global search and stochastic local search. In addition, we hybridize the three derivative-free algorithms with a multiscale regularization approach. Starting with a reasonably small number of control steps, the control intervals are subsequently refined during the optimization. Results for experiments studied indicate that CMA-ES performs best among the three algorithms in solving both small and large scale problems. When hybridized with a multiscale regularization approach, the ability to find the optimal solution is further enhanced, with the performance of GPS improving the most. Topics affecting the performance of the multiscale approach are discussed in this paper, including the effect of control frequency on the well control problem. The parameter settings for GPS, PSO, and CMA-ES, within the multiscale approach are considered.


title: Optimization of Spatial Dose Distribution for Controlling Sidewall Shape in Electron-beam Lithography

abstract: Electron-beam (e-beam) lithography is widely employed in fabrication of 2-D patterns and 3-D structures. A certain type or shape of the sidewall in the remaining resist profile may be desired in an application, e.g., an undercut for lift-off and a vertical sidewall for etching, or required for a device. Also, as the feature size is decreased well below a micron, a small variation of the sidewall slope can lead to a significant (relative) CD error in certain layers of resist. Therefore, it is important to understand effects of spatial dose distribution on sidewall shape and be able to achieve the desired shape. In this study, via simulation, the relationship among the total dose, spatial distribution of dose, developing time and sidewall shape, and performance of the method developed to optimize the dose distribution for a target sidewall shape have been analyzed. The simulation results have been verified through experiments.


title: A Multivariate Cure Model for Left- and Right-Censored Data with Application to Colorectal Cancer Screening Patterns

abstract: We develop a multivariate cure survival model to estimate lifetime patterns of colorectal cancer screening. Screening data cover long periods of time, with sparse observations for each person. Some events may occur before the study begins or after the study ends, so the data are both left- and right-censored, and some individuals are never screened (the "cured" population). We propose a multivariate parametric cure model that can be used with left- and right-censored data. Our model allows for the estimation of the time to screening as well as the average number of times individuals will be screened. We calculate likelihood functions based on the observations for each subject using a distribution that accounts for within-subject correlation, and estimate parameters using Markov Chain Monte Carlo methods. We apply our methods to the estimation of lifetime colorectal cancer screening behavior in the SEER-Medicare data set.


title: On the cop number of generalized Petersen graphs

abstract: We show that the cop number of every generalized Petersen graph is at most 4. The strategy is to play a modified game of cops and robbers on an infinite cyclic covering space where the objective is to capture the robber or force the robber towards an end of the infinite graph. We prove that finite isometric subtrees are 1-guardable and apply this to determine the exact cop number of some families of generalized Petersen graphs. We also extend these ideas to prove that the cop number of any connected I-graph is at most 5.


title: Sensitivity to initial conditions of a $d$-dimensional long-range-interacting quartic Fermi-Pasta-Ulam model: Universal scaling

abstract: We introduce a generalized $d$-dimensional Fermi-Pasta-Ulam (FPU) model in presence of long-range interactions, and perform a first-principle study of its chaos for $d=1,2,3$ through large-scale numerical simulations. The nonlinear interaction is assumed to decay algebraically as $d_{ij}^{-\alpha}$ ($\alpha \ge 0$), $\{d_{ij}\}$ being the distances between $N$ oscillator sites. Starting from random initial conditions we compute the maximal Lyapunov exponent $\lambda_{max}$ as a function of $N$. Our $N>>1$ results strongly indicate that $\lambda_{max}$ remains constant and positive for $\alpha/d>1$ (implying strong chaos, mixing and ergodicity), and that it vanishes like $N^{-\kappa}$ for $0 \le \alpha/d < 1$ (thus approaching weak chaos and opening the possibility of breakdown of ergodicity). The suitably rescaled exponent $\kappa$ exhibits universal scaling, namely that $(d+2) \kappa$ depends only on $\alpha/d$ and, when $\alpha/d$ increases from zero to unity, it monotonically decreases from unity to zero, remaining so for all $\alpha/d >1$. The value $\alpha/d=1$ can therefore be seen as a critical point separating the ergodic regime from the anomalous one, $\kappa$ playing a role analogous to that of an order parameter. This scaling law is consistent with Boltzmann-Gibbs statistics for $\alpha/d > 1$, and possibly with $q$-statistics for $0 \le \alpha/d < 1$.


title: Optimal Policies for Wireless Networks with Energy Harvesting Transmitters and Receivers: Effects of Decoding Costs

abstract: We consider the effects of decoding costs in energy harvesting communication systems. In our setting, receivers, in addition to transmitters, rely solely on energy harvested from nature, and need to spend some energy in order to decode their intended packets. We model the decoding energy as an increasing convex function of the rate of the incoming data. In this setting, in addition to the traditional energy causality constraints at the transmitters, we have the decoding causality constraints at the receivers, where energy spent by the receiver for decoding cannot exceed its harvested energy. We first consider the point-to-point single-user problem where the goal is to maximize the total throughput by a given deadline subject to both energy and decoding causality constraints. We show that decoding costs at the receiver can be represented as generalized data arrivals at the transmitter, and thereby moving all system constraints to the transmitter side. Then, we consider several multi-user settings. We start with a two-hop network where the relay and the destination have decoding costs, and show that separable policies, where the transmitter's throughput is maximized irrespective of the relay's transmission energy profile, are optimal. Next, we consider the multiple access channel (MAC) and the broadcast channel (BC) where the transmitters and the receivers harvest energy from nature, and characterize the maximum departure region. In all multi-user settings considered, we decompose our problems into inner and outer problems. We solve the inner problems by exploiting the structure of the particular model, and solve the outer problems by water-filling algorithms.


title: Confluence of Layered Rewrite Systems

abstract: We investigate the new, Turing-complete class of layered systems, whose lefthand sides of rules can only be overlapped at a multiset of disjoint or equal positions. Layered systems define a natural notion of rank for terms: the maximal number of non-overlapping redexes along a path from the root to a leaf. Overlappings are allowed in finite or infinite trees. Rules may be non-terminating, non-left-linear, or non-right-linear. Using a novel unification technique, cyclic unification, and the so-alled subrewriting relation, we show that rank non-increasing layered systems are confluent provided their cyclic critical pairs have cyclic-joinable decreasing diagrams.


title: One-nucleon transfer reactions and the optical potential

abstract: We provide a summary of new developments in the area of direct reaction theory with a particular focus on one-nucleon transfer reactions. We provide a status of the methods available for describing (d,p) reactions. We discuss the effects of nonlocality in the optical potential in transfer reactions. The results of a purely phenomenological potential and the optical potential obtained from the dispersive optical model are compared; both point toward the importance of including nonlocality in transfer reactions explicitly. Given the large ambiguities associated with optical potentials, we discuss some new developments toward the quantification of this uncertainty. We conclude with some general comments and a brief account of new advances that are in the pipeline.


title: Subdivisions in the Robber Locating Game

abstract: We consider a game in which a cop searches for a moving robber on a graph using distance probes, which is a slight variation on one introduced by Seager. Carragher, Choi, Delcourt, Erickson and West showed that for any n-vertex graph $G$ there is a winning strategy for the cop on the graph $G^{1/m}$ obtained by replacing each edge of $G$ by a path of length $m$, if $m \geqslant n$. They conjectured that this bound was best possible for complete graphs, but the present authors showed that in fact the cop wins on $K^{1/m}$ if and only if $m \geqslant n/2$, for all but a few small values of $n$. In this paper we extend this result to general graphs by proving that the cop has a winning strategy on $G^{1/m}$ provided $m \geqslant n/2$ for all but a few small values of $n$; this bound is best possible. We also consider replacing the edges of $G$ with paths of varying lengths.


title: High-order Tail in Schwarzschild Space-time

abstract: We present an analysis of the behaviour at late-times of linear field perturbations of a Schwarzschild black hole space-time. In particular, we give explicit analytic expressions for the field perturbations (for a specific multipole) of general spin up to the first four orders at late times. These expressions are valid at arbitrary radius and include, apart from the well-known power-law tail decay at leading order ($\sim t^{-2\ell-3}$), a new logarithmic behaviour at third leading order ($\sim t^{-2\ell-5}\ln t$). We obtain these late-time results by developing the so-called MST formalism and by expanding the various MST Fourier-mode quantities for small frequency. While we give explicit expansions up to the first four leading orders (for small-frequency for the Fourier modes, for late-time for the field perturbation), we give a prescription for obtaining expressions to arbitrary order within a `perturbative regime'.


title: Newtonian, post-Newtonian and Relativistic Cosmological Perturbation Theory

abstract: Newtonian cosmological perturbation equations valid to full nonlinear order are well known in the literature. Assuming the absence of the transverse-tracefree part of the metric, we present the general relativistic counterpart valid to full nonlinear order. The relativistic equations are presented without taking the slicing (temporal gauge) condition. The equations do have the proper Newtonian and first post-Newtonian limits. We also present the relativistic pressure correction terms in the Newtonian hydrodynamic equations.


title: Central limit theorems for network driven sampling

abstract: Respondent-Driven Sampling is a popular technique for sampling hidden populations. This paper models Respondent-Driven Sampling as a Markov process indexed by a tree. Our main results show that the Volz-Heckathorn estimator is asymptotically normal below a critical threshold. The key technical difficulties stem from (i) the dependence between samples and (ii) the tree structure which characterizes the dependence. The theorems allow the growth rate of the tree to exceed one and suggest that this growth rate should not be too large. To illustrate the usefulness of these results beyond their obvious use, an example shows that in certain cases the sample average is preferable to inverse probability weighting. We provide a test statistic to distinguish between these two cases.


title: Forecasting Method for Grouped Time Series with the Use of k-Means Algorithm

abstract: The paper is focused on the forecasting method for time series groups with the use of algorithms for cluster analysis. $K$-means algorithm is suggested to be a basic one for clustering. The coordinates of the centers of clusters have been put in correspondence with summarizing time series data the centroids of the clusters. A description of time series, the centroids of the clusters, is implemented with the use of forecasting models. They are based on strict binary trees and a modified clonal selection algorithm. With the help of such forecasting models, the possibility of forming analytic dependences is shown. It is suggested to use a common forecasting model, which is constructed for time series the centroid of the cluster, in forecasting the private (individual) time series in the cluster. The promising application of the suggested method for grouped time series forecasting is demonstrated.


title: Direct high-order edge-preserving regularization for tomographic image reconstruction

abstract: In this paper we present a new two-level iterative algorithm for tomographic image reconstruction. The algorithm uses a regularization technique, which we call edge-preserving Laplacian, that preserves sharp edges between objects while damping spurious oscillations in the areas where the reconstructed image is smooth. Our numerical simulations demonstrate that the proposed method outperforms total variation (TV) regularization and it is competitive with the combined TV-L2 penalty. Obtained reconstructed images show increased signal-to-noise ratio and visually appealing structural features. Computer implementation and parameter control of the proposed technique is straightforward, which increases the feasibility of it across many tomographic applications. In this paper, we applied our method to the under-sampled computed tomography (CT) projection data and also considered a case of reconstruction in emission tomography The MATLAB code is provided to support obtained results.


title: KMOS view of the Galactic Centre I. Young stars are centrally concentrated

abstract: The Galactic centre hosts a crowded, dense nuclear star cluster with a half-light radius of 4 pc. Most of the stars in the Galactic centre are cool late-type stars, but there are also >100 hot early-type stars in the central parsec of the Milky Way. These stars are only 3-8 Myr old. Our knowledge of the number and distribution of early-type stars in the Galactic centre is incomplete. Only a few spectroscopic observations have been made beyond a projected distance of 0.5 pc of the Galactic centre. The distribution and kinematics of early-type stars are essential to understand the formation and growth of the nuclear star cluster. We cover the central >4pc^2 of the Galactic centre using the integral-field spectrograph KMOS. We extracted more than 1,000 spectra from individual stars and identified early-type stars based on their spectra. Our data set contains 114 bright early-type stars: 6 have narrow emission lines, 23 are Wolf-Rayet stars, 9 stars have featureless spectra, and 76 are O/B type stars. Our wide-field spectroscopic data confirm that the distribution of young stars is compact, with 90% of the young stars identified within 0.5 pc of the nucleus. We identify 24 new O/B stars primarily at large radii. We estimate photometric masses of the O/B stars and show that the total mass in the young population is >12,000M_sun. The O/B stars all appear to be bound to the Milky Way nuclear star cluster, while less than 30% belong to the clockwise rotating disk. The central concentration of the early-type stars is a strong argument that they have formed in situ. A large part of the young O/B stars is not on the disk, which either means that the early-type stars did not all form on the same disk or that the disk is dissolving rapidly. [abridged]


title: No evidence for bulk velocity from type Ia supernovae

abstract: We revisit the effect of peculiar velocities on low-redshift type Ia supernovae. Velocities introduce an additional guaranteed source of correlations between supernova magnitudes that should be considered in all analyses of nearby supernova samples but has largely been neglected in the past. Applying a likelihood analysis to the latest compilation of nearby supernovae, we find no evidence for the presence of these correlations, although, given the significant noise, the data is also consistent with the correlations predicted for the standard $\Lambda$CDM model. We then consider the dipolar component of the velocity correlations -- the frequently studied "bulk velocity" -- and explicitly demonstrate that including the velocity correlations in the data covariance matrix is crucial for drawing correct and unambiguous conclusions about the bulk flow. In particular, current supernova data is consistent with no excess bulk flow on top of what is expected in $\Lambda$CDM and effectively captured by the covariance. We further clarify the nature of the apparent bulk flow that is inferred when the velocity covariance is ignored. We show that a significant fraction of this quantity is expected to be noise bias due to uncertainties in supernova magnitudes and not any physical peculiar motion.


title: Pair Correlations in Doped Hubbard Ladders

abstract: Hubbard ladders are an important stepping stone to the physics of the two-dimensional Hubbard model. While many of their properties are accessible to numerical and analytical techniques, the question of whether weakly hole-doped Hubbard ladders are dominated by superconducting or charge-density-wave correlations has so far eluded a definitive answer. In particular, previous numerical simulations of Hubbard ladders have seen a much faster decay of superconducting correlations than expected based on analytical arguments. We revisit this question using a state-of-the-art implementation of the density matrix renormalization group algorithm that allows us to simulate larger system sizes with higher accuracy than before. Performing careful extrapolations of the results, we obtain improved estimates for the Luttinger liquid parameter and the correlation functions at long distances. Our results confirm that, as suggested by analytical considerations, superconducting correlations become dominant in the limit of very small doping.


title: NGC 3628-UCD1: A possible $\omega$ Cen Analog Embedded in a Stellar Stream

abstract: Using Subaru/Suprime-Cam wide-field imaging and both Keck/ESI and LBT/MODS spectroscopy, we identify and characterize a compact star cluster, which we term NGC 3628-UCD1, embedded in a stellar stream around the spiral galaxy NGC 3628. The size and luminosity of UCD1 are similar to $\omega$ Cen, the most luminous Milky Way globular cluster, which has long been suspected to be the stripped remnant of an accreted dwarf galaxy. The object has a magnitude of $i=19.3$ mag (${\rm L}_{\rm i}=1.4\times10^{6}~{\rm L}_{\odot}$). UCD1 is marginally resolved in our ground-based imaging, with a half-light radius of $\sim10$ pc. We measure an integrated brightness for the stellar stream of $i=13.1$ mag, with $(g-i)=1.0$. This would correspond to an accreted dwarf galaxy with an approximate luminosity of ${\rm L}_i\sim4.1\times10^{8}~{\rm L}_{\odot}$. Spectral analysis reveals that UCD1 has an age of $6.6$ Gyr , $[\rm{Z}/\rm{H}]=-0.75$, an $[{\alpha}/\rm{Fe}]=-0.10$. We propose that UCD1 is an example of an $\omega$ Cen-like star cluster possibly forming from the nucleus of an infalling dwarf galaxy, demonstrating that at least some of the massive star cluster population may be created through tidal stripping.


title: Interpretations of quantum theory: A map of madness

abstract: Motivated by some recent news, a journalist asks a group of physicists: "What's the meaning of the violation of Bell's inequality?" One physicist answers: "It means that non-locality is an established fact". Another says: "There is no non-locality; the message is that measurement outcomes are irreducibly random". A third one says: "It cannot be answered simply on purely physical grounds, the answer requires an act of metaphysical judgement". Puzzled by the answers, the journalist keeps asking questions about quantum theory: "What is teleported in quantum teleportation?" "How does a quantum computer really work?" Shockingly, for each of these questions, the journalist obtains a variety of answers which, in many cases, are mutually exclusive. At the end of the day, the journalist asks: "How do you plan to make progress if, after 90 years of quantum theory, you still don't know what it means? How can you possibly identify the physical principles of quantum theory or expand quantum theory into gravity if you don't agree on what quantum theory is about?" Here we argue that it is becoming urgent to solve this too long lasting problem. For that, we point out that the interpretations of quantum theory are, essentially, of two types and that these two types are so radically different that there must be experiments that, when analyzed outside the framework of quantum theory, lead to different empirically testable predictions. Arguably, even if these experiments do not end the discussion, they will add new elements to the list of strange properties that some interpretations must have, therefore they will indirectly support those interpretations that do not need to have all these strange properties.


title: A Pan-STARRS1 study of the relationship between wide binarity and planet occurrence in the Kepler field

abstract: The NASA Kepler mission has revolutionised time-domain astronomy and has massively expanded the number of known extrasolar planets. However, the effect of wide multiplicity on exoplanet occurrence has not been tested with this dataset. We present a sample of 401 wide multiple systems containing at least one Kepler target star. Our method uses Pan-STARRS1 and archival data to produce an accurate proper motion catalogue of the Kepler field. Combined with Pan-STARRS1 SED fits and archival proper motions for bright stars, we use a newly developed probabilistic algorithm to identify likely wide binary pairs which are not chance associations. As by-products of this we present stellar SED templates in the Pan-STARRS1 photometric system and conversions from this system to Kepler magnitudes. We find that Kepler target stars in our binary sample with separations above 6 arcseconds are no more or less likely to be identified as confirmed or candidate planet hosts than a weighted comparison sample of Kepler targets of similar brightness and spectral type. Therefore we find no evidence that binaries with projected separations greater than 3,000AU affect the occurrence rate of planets with P<300days around FGK stars.


title: Observation of the Dirac fluid and the breakdown of the Wiedemann-Franz law in graphene

abstract: Interactions between particles in quantum many-body systems can lead to collective behavior described by hydrodynamics. One such system is the electron-hole plasma in graphene near the charge neutrality point which can form a strongly coupled Dirac fluid. This charge neutral plasma of quasi-relativistic fermions is expected to exhibit a substantial enhancement of the thermal conductivity, due to decoupling of charge and heat currents within hydrodynamics. Employing high sensitivity Johnson noise thermometry, we report the breakdown of the Wiedemann-Franz law in graphene, with a thermal conductivity an order of magnitude larger than the value predicted by Fermi liquid theory. This result is a signature of the Dirac fluid, and constitutes direct evidence of collective motion in a quantum electronic fluid.


title: Sensitivity and foreground modelling for large-scale CMB B-mode polarization satellite missions

abstract: The measurement of the large-scale B-mode polarization in the cosmic microwave background (CMB) is a fundamental goal of future CMB experiments. However, because of unprecedented sensitivity, future CMB experiments will be much more sensitive to any imperfect modelling of the Galactic foreground polarization in the reconstruction of the primordial B-mode signal. We compare the sensitivity to B-modes of different concepts of CMB satellite missions (LiteBIRD, COrE, COrE+, PRISM, EPIC, PIXIE) in the presence of Galactic foregrounds. In particular, we quantify the impact on the tensor-to-scalar parameter of incorrect foreground modelling in the component separation process. Using Bayesian fitting and Gibbs sampling, we perform the separation of the CMB and Galactic foreground B-modes. The recovered CMB B-mode power spectrum is used to compute the likelihood distribution of the tensor-to-scalar ratio. We focus the analysis to the very large angular scales that can be probed only by CMB space missions, i.e. the Reionization bump, where primordial B-modes dominate over spurious B-modes induced by gravitational lensing. We find that fitting a single modified blackbody component for thermal dust where the "real" sky consists of two dust components strongly bias the estimation of the tensor-to-scalar ratio by more than 5{\sigma} for the most sensitive experiments. Neglecting in the parametric model the curvature of the synchrotron spectral index may bias the estimated tensor-to-scalar ratio by more than 1{\sigma}. For sensitive CMB experiments, omitting in the foreground modelling a 1% polarized spinning dust component may induce a non-negligible bias in the estimated tensor-to-scalar ratio.


title: Successive magnetic field-induced transitions and colossal magnetoelectric effect in Ni$_{3}$TeO$_{6}$

abstract: We report the discovery of a metamagnetic phase transition in a polar antiferromagnet Ni$_3$TeO$_6$ that occurs at 52 T. The new phase transition accompanies a colossal magnetoelectric effect, with a magnetic-field-induced polarization change of 0.3 $\mu$C/cm$^2$, a value that is 4 times larger than for the spin-flop transition at 9 T in the same material, and also comparable to the largest magnetically-induced polarization changes observed to date. Via density-functional calculations we construct a full microscopic model that describes the data. We model the spin structures in all fields and clarify the physics behind the 52 T transition. The high-field transition involves a competition between multiple different exchange interactions which drives the polarization change through the exchange-striction mechanism. The resultant spin structure is rather counter-intuitive and complex, thus providing new insights on design principles for materials with strong magnetoelectric coupling.


title: The Arches Cluster: Extended Structure and Tidal Radius

abstract: At a projected distance of ~26 pc from Sgr A*, the Arches cluster provides insight to star formation in the extreme Galactic Center (GC) environment. Despite its importance, many key properties such as the cluster's internal structure and orbital history are not well known. We present an astrometric and photometric study of the outer region of the Arches cluster (R > 6.25") using HST WFC3IR. Using proper motions we calculate membership probabilities for stars down to F153M = 20 mag (~2.5 M_sun) over a 120" x 120" field of view, an area 144 times larger than previous astrometric studies of the cluster. We construct the radial profile of the Arches to a radius of 75" (~3 pc at 8 kpc), which can be well described by a single power law. From this profile we place a 3-sigma lower limit of 2.8 pc on the observed tidal radius, which is larger than the predicted tidal radius (1 - 2.5 pc). Evidence of mass segregation is observed throughout the cluster and no tidal tail structures are apparent along the orbital path. The absence of breaks in the profile suggests that the Arches has not likely experienced its closest approach to the GC between ~0.2 - 1 Myr ago. If accurate, this constraint indicates that the cluster is on a prograde orbit and is located front of the sky plane that intersects Sgr A*. However, further simulations of clusters in the GC potential are required to interpret the observed profile with more confidence.


title: Electron Thermodynamics in GRMHD Simulations of Low-Luminosity Black Hole Accretion

abstract: Simple assumptions made regarding electron thermodynamics often limit the extent to which general relativistic magnetohydrodynamic (GRMHD) simulations can be applied to observations of low-luminosity accreting black holes. We present, implement, and test a model that self-consistently evolves an electron entropy equation and takes into account the effects of spatially varying electron heating and relativistic anisotropic thermal conduction along magnetic field lines. We neglect the back-reaction of electron pressure on the dynamics of the accretion flow. Our model is appropriate for systems accreting at $\ll 10^{-5}$ of the Eddington rate, so radiative cooling by electrons can be neglected. It can be extended to higher accretion rates in the future by including electron cooling and proton-electron Coulomb collisions. We present a suite of tests showing that our method recovers the correct solution for electron heating under a range of circumstances, including strong shocks and driven turbulence. Our initial applications to axisymmetric simulations of accreting black holes show that (1)~physically-motivated electron heating rates yield electron temperature distributions significantly different from the constant electron to proton temperature ratios assumed in previous work, with higher electron temperatures concentrated in the coronal region between the disc and the jet; (2)~electron thermal conduction significantly modifies the electron temperature in the inner regions of black hole accretion flows if the effective electron mean free path is larger than the local scale-height of the disc (at least for the initial conditions and magnetic field configurations we study). The methods developed in this work are important for producing more realistic predictions for the emission from accreting black holes such as Sagittarius A* and M87; these applications will be explored in future work.


title: Null Fluids - A New Viewpoint of Galilean Fluids

abstract: This article is a detailed version of our short letter `On equilibrium partition function for non-relativistic fluid' [arXiv:1505.05677] extended to include an anomalous $U(1)$ symmetry. We construct a relativistic system, which we call null fluid and show that it is in one-to-one correspondence with a Galilean fluid living in one lower dimension. The correspondence is based on light cone reduction, which is known to reduce the Poincare symmetry of a theory to Galilean in one lower dimension. We show that the proposed null fluid and the corresponding Galilean fluid have exactly same symmetries, thermodynamics, constitutive relations, and equilibrium partition to all orders in derivative expansion. We also devise a mechanism to introduce $U(1)$ anomaly in even dimensional Galilean theories using light cone reduction, and study its effect on the constitutive relations of a Galilean Fluid.


title: Spectral properties for a family of two-dimensional quantum antiferromagnets

abstract: We study the spectral properties of a family of quantum antiferromagnets on two-dimensional (2D) lattices. This family of models is obtained by a deformation of the well-studied 2D quantum antiferromagnetic model of Affleck, Kennedy, Lieb and Tasaki (AKLT); they are described by two-body, frustration-free Hamiltonians on a three-colourable lattice of spins. Although the existence of a spectral gap in the 2D AKLT model remains an open question, we rigorously prove the existence of a gap for a subset of this family of quantum antiferromagnets. Along with providing new progress for the gap problem in AKLT-type antiferromagnets in 2D, this result has implications for the theory of quantum computation as it provides a family of two-body Hamiltonians for which the ground state is a resource for universal quantum computation and for which a spectral gap is proven to exist.


title: The inferred evolution of the cold gas properties of CANDELS galaxies at 0.5 < z < 3.0

abstract: We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.


title: Ground state on the dumbbell graph

abstract: We consider standing waves in the focusing nonlinear Schr\"odinger (NLS) equation on a dumbbell graph (two rings attached to a central line segment subject to the Kirchhoff boundary conditions at the junctions). In the limit of small $L^2$ norm, the ground state (the orbitally stable standing wave of the smallest energy at a fixed $L^2$ norm) is represented by a constant solution. However, when the $L^2$ norm is increased, this constant solution undertakes two bifurcations, where the first is the pitchfork (symmetry breaking) bifurcation and the second one is the symmetry preserving bifurcation. As a result of the first symmetry breaking bifurcation, the standing wave becomes more localized in one of the two rings. As a result of the second symmetry preserving bifurcation, the standing wave becomes localized in the central line segment. In the limit of large norm solutions, both standing waves are represented by a truncated solitary wave localized in either the ring or the central line segment. Although the asymmetric wave supported in the ring is a ground state near the symmetry breaking bifurcation of the constant solution, it is the symmetric wave supported in the central line segment which becomes the ground state in the limit of large $L^2$ norm. The analytical results are confirmed by numerical approximations of the ground state on the dumbbell graph.


title: Nef cones of Hilbert schemes of points on surfaces

abstract: Let X be a smooth projective surface of irregularity 0. The Hilbert scheme of n points on X parameterizes zero-dimensional subschemes of X of length n. In this paper, we discuss general methods for studying the cone of ample divisors on the Hilbert scheme. We then use these techniques to compute the cone of ample divisors on the Hilbert scheme for several surfaces where the cone was previously unknown. Our examples include families of surfaces of general type and del Pezzo surfaces of degree 1. The methods rely on Bridgeland stability and the Positivity Lemma of Bayer and Macri.


title: Crystal structure and magnetism in $\alpha$-RuCl3: an ab-initio study

abstract: $\alpha$-RuCl$_3$ has been proposed recently as an excellent playground for exploring Kitaev physics on a two-dimensional (2D) honeycomb lattice. However, structural clarification of the compound has not been completed, which is crucial in understanding the physics of this system. Here, using {\it ab-initio} electronic structure calculations, we study a full three dimensional (3D) structure of $\alpha$-RuCl$_3$ including the effects of spin-orbit coupling (SOC) and electronic correlations. Three major results are as follows; i) SOC suppresses dimerization of Ru atoms, which exists in other Ru compounds such as isostructural Li$_2$RuO$_3$, and making the honeycomb closer to an ideal one. ii) The nearest-neighbor Kitaev exchange interaction between the $j_{\rm eff}$=1/2 pseudospin depends strongly on the Ru-Ru distance and the Cl position, originating from the nature of the edge-sharing geometry. iii) The optimized 3D structure without electronic correlations has $P{\bar 3}1m$ space group symmetry independent of SOC, but including electronic correlation changes the optimized 3D structure to either $C2/m$ or $Cmc2_1$ within 0.1 meV per formula unit (f.u.) energy difference. The reported $P3_112$ structure is also close in energy. The interlayer spin exchange coupling is a few percent of in-plane spin exchange terms, confirming $\alpha$-RuCl$_3$ is close to a 2D system. We further suggest how to increase the Kitaev term via tensile strain, which sheds new light in realizing Kitaev spin liquid phase in this system.


title: Heavy Resonances in the Electroweak Effective Lagrangian

abstract: As a first step towards the construction of a general electroweak effective Lagrangian incorporating heavy states, we present here a simplified version where only vector and axial-vector spin-1 triplets are involved. We adopt an effective field theory formalism, implementing the electroweak chiral symmetry breaking SU(2)L x SU(2)R -> SU(2)L+R, which couples the heavy states to the SM fields. At low energies, the heavy degrees of freedom are integrated out from the action and their effects are hidden in the low-energy couplings of the Electroweak Effective Theory, which can be tested experimentally. Short-distance constraints are also implemented, requiring a proper behaviour in the high-energy regime. We analyze the phenomenological constraints from the oblique parameters S and T , at the next- to-leading order. Our results show that present data allow for strongly-coupled scenarios with massive bosons above the TeV scale.


title: Resonances in W_L W_L, Z_L Z_L and hh scattering from dispersive analysis of the non-linear Electroweak+Higgs Effective Theory

abstract: If new resonances of the electroweak symmetry breaking sector (longitudinal-gauge and Higgs) bosons are found in the 1-3 TeV region, the right tool to assess their properties and confront experimental data in a largely model-independent yet simple manner is Unitarized Effective Theory. Its ingredients are: 1) custodial symmetry and the Equivalence Theorem, that allow to approximate W_L and Z_L by an isospin-triplet of Goldstone bosons omega^a in the 1-TeV region. 2) The effective coupling of a generic, approximately massless scalar-isoscalar h to those Goldstone bosons, and the chiral Lagrangian describing them, valid up to about 3 TeV. 3) The Inverse Amplitude or other unitarization techniques that allow to extend the reach of perturbation theory to the first resonance in each partial wave. We highlight some of the parameter space that can give rise to 2-TeV resonances, for example a simultaneous scalar-isoscalar and a vector-isovector ones (motivated by the ATLAS excess) and also the potential importance of coupled-channel dynamics between hh and omega omega.


title: Beautiful mirrors for a pNGB Higgs

abstract: We consider one of the most significant deviations from the Standard Model: the forward-backward asymmetry of the b-quark measured at leptonic colliders. We investigate the possibility to solve this discrepancy by introducing new physics at the TeV scale. We focus on models where the Higgs is a pseudo Nambu-Goldstone boson of a new strongly coupled sector with a global SO(5) symmetry broken spontaneously to SO(4). Besides the usual top partners, we introduce bottom partners in the representations 16 and 4 of SO(5) and show that they can improve significantly the fit by correcting the Zbb couplings. We also estimate the corrections to the couplings at one-loop and obtain that the tree-level ones dominate and can give a reliable estimation. We find that the large shift required for Zb_Rb_R leads to light custodians associated to the b-quark, similar to the top partners, as well as a rich phenomenology involving neutral interactions in the bottom-sector.


title: Laplace Symbols and Invariant Distributions

abstract: We introduce a new kind of symbol in the framework of It\^o processes which are bounded on one side. The connection between this symbol and the infinitesimal generator is analyzed. Based on this concept, an integral criterion for invariant distributions of the underlying process is derived. Some applications are mentioned.


title: Detecting topological sectors in continuum Yang-Mills theory and the fate of BRST symmetry

abstract: In this work, motivated by Laplacian type center gauges in the lattice, designed to avoid the Gribov problem, we introduce a new family of gauge fixings for pure Yang-Mills theories in the continuum. This procedure separates the partition function into partial contributions associated with different sectors, containing center vortices and correlated monopoles. We show that, on each sector, the gauge fixed path-integral displays a BRST symmetry, however, it cannot be globally extended due to sector dependent boundary conditions on the ghost fields. These are nice features as they would permit to discuss the independence of the partial contributions on gauge parameters,, while opening a window for the space of quantum states to be different from the perturbative one, which would be implied if topological configurations were removed.


title: Integrating Research Data Management into Geographical Information Systems

abstract: Ocean modelling requires the production of high-fidelity computational meshes upon which to solve the equations of motion. The production of such meshes by hand is often infeasible, considering the complexity of the bathymetry and coastlines. The use of Geographical Information Systems (GIS) is therefore a key component to discretising the region of interest and producing a mesh appropriate to resolve the dynamics. However, all data associated with the production of a mesh must be provided in order to contribute to the overall recomputability of the subsequent simulation. This work presents the integration of research data management in QMesh, a tool for generating meshes using GIS. The tool uses the PyRDM library to provide a quick and easy way for scientists to publish meshes, and all data required to regenerate them, to persistent online repositories. These repositories are assigned unique identifiers to enable proper citation of the meshes in journal articles.


title: Rare semileptonic decay of $\chi_{c1}(1p)$ meson In QCD

abstract: The rare semileptonic $\chi_{c1}(1p)\rightarrow D_{s}^{+}e\bar{\nu }$ decay is analyzed, by using the three-point QCD sum rules. Taking into account the two-gluon condensate contributions, the transition form factors related to this decay are calculated and are used to determine the total decay width and branching fraction. Our findings may be approved by future experiments.


title: Beyond packing of hard spheres: The effects of core softness, non-additivity, intermediate-range repulsion, and many-body interactions on the glass-forming ability of bulk metallic glasses

abstract: When a liquid is cooled well below its melting temperature at a rate that exceeds the critical cooling rate $R_c$, the crystalline state is bypassed and an amorphous glassy state forms instead. $R_c$ (or the corresponding critical casting thickness $d_c$) characterizes the glass-forming ability (GFA) of each material. While silica is an excellent glass-former with small $R_c<10^{-2}$ K/s, pure metals and most alloys are poor glass-formers with large $R_c>10^{10}$ K/s. Only in the past thirty years have bulk metallic glasses (BMGs) been identified with $R_c$ approaching that for silica. Recent simulations have shown that hard-sphere models are able to identify the atomic size ratio and number fraction regime where BMGs exist with critical cooling rates more than 13 orders of magnitude smaller than those for pure metals. However, there are many other features of interatomic potentials beyond hard-core interactions. How do these other features affect the glass-forming ability of BMGs? We perform molecular dynamics simulations to determine how variations in the softness and non-additivity of the repulsive core and form of the interatomic pair potential at intermediate distances affect the GFA of binary alloys. These variations in the interatomic pair potential allow us to introduce geometric frustration and change the crystal phases that compete with glass formation. We also investigate the effect of tuning the strength of the many-body interactions from zero to the full embedded atom model on the GFA for pure metals. We then employ the full embedded atom model for binary BMGs and show that hard-core interactions play the dominant role in setting the GFA of alloys, while other features of the interatomic potential only change the GFA by one to two orders of magnitude.


title: Background subtraction methods for precision measurements of di-hadron and jet-hadron correlations in heavy ion collisions

abstract: Di-hadron correlations and jet-hadron correlations are frequently used to study interactions of the Quark Gluon Plasma with the hard partons that form jets. The existing background subtraction methods for these studies depend on several assumptions and independent measurements of the Fourier coefficients of the combinatorial background. In this paper, we present two background subtraction methods which use the observation that the near-side correlation signal is negligible for large separations between the associated particle and the trigger. The combinatorial background is fit in this background-dominated region. Including the reaction plane dependence of the raw correlation improves these results. To test the accuracy of these methods, a simple model is used to simulate di-hadron and jet-hadron correlations with a combinatorial background similar to that observed in the data. The true signal is compared to the extracted signal. The results are compared to results from two variants of the Zero-Yield-At-Minimum (ZYAM) method. We test these methods for mid-peripheral and central collisions and for di-hadron and jet-hadron correlations. These methods are more precise than the ZYAM method with fewer assumptions about the shape and level of the combinatorial background, even in central collisions where the experimental resolution on the measurement of the reaction plane dependence is poor. These methods will allow more accurate studies of modifications of the away-side jet and will be particularly useful for studies of jet-hadron correlations, where the combinatorial background is poorly constrained from previous studies.


title: Factorization threshold models for scale-free networks generation

abstract: Many real networks such as the World Wide Web, financial, biological, citation and social networks have a power-law degree distribution. Networks with this feature are also called scale-free. Several models for producing scale-free networks have been obtained by now and most of them are based on the preferential attachment approach. We will offer the model with another scale-free property explanation. The main idea is to approximate the network's adjacency matrix by multiplication of the matrices $V$ and $V^T$, where $V$ is the matrix of vertices' latent features. This approach is called matrix factorization and is successfully used in the link prediction problem. To create a generative model of scale-free networks we will sample latent features $V$ from some probabilistic distribution and try to generate a network's adjacency matrix. Entries in the generated matrix are dot products of latent features which are real numbers. In order to create an adjacency matrix, we approximate entries with the Boolean domain $\{0, 1\}$. We have incorporated the threshold parameter $\theta$ into the model for discretization of a dot product. Actually, we have been influenced by the geographical threshold models which were recently proven to have good results in a scale-free networks generation. The overview of our results is the following. First, we will describe our model formally. Second, we will tune the threshold $\theta$ in order to generate sparse growing networks. Finally, we will show that our model produces scale-free networks with the fixed power-law exponent which equals two. In order to generate oriented networks with tunable power-law exponents and to obtain other model properties, we will offer different modifications of our model. Some of our results will be demonstrated using computer simulation.


title: Information transfer in community structured multiplex networks

abstract: The study of complex networks that account for different types of interactions has become a subject of interest in the last few years, specially because its representational power in the description of users interactions in diverse online social platforms (Facebook, Twitter, Instagram, etc.). The mathematical description of these interacting networks has been coined under the name of multilayer networks, where each layer accounts for a type of interaction. It has been shown that diffusive processes on top of these networks present a phenomenology that cannot be explained by the naive superposition of single layer diffusive phenomena but require the whole structure of interconnected layers. Nevertheless, the description of diffusive phenomena on multilayer networks has obviated the fact that social networks have strong mesoscopic structure represented by different communities of individuals driven by common interests, or any other social aspect. In this work, we study the transfer of information in multilayer networks with community structure. The final goal is to understand and quantify, if the existence of well-defined community structure at the level of individual layers, together with the multilayer structure of the whole network, enhances or deteriorates the diffusion of packets of information.