title: A numerical method for the solution of relaxed one-sided Lipschitz algebraic inclusions

abstract: An existing solvability result for relaxed one-sided Lipschitz algebraic inclusions is substantially improved. This enhanced solvability result allows the design of a very robust numerical method for the approximation of a solution of the algebraic inclusion. Sharp error estimates for this method, illustrative analytic examples and a numerical example are provided.


title: Isomorphic universality and the number of pairwise non-isomorphic models in the class of Banach spaces

abstract: We study isomorphic universality of Banach spaces of a given density and a number of pairwise non-isomorphic models in the same class. We show that in the Cohen model the isomorphic universality number for Banach spaces of density $\aleph_1$ is $\aleph_2$, and analogous results are true for other cardinals (Theorem 1.2(1)) and that adding just one Cohen real to any model destroys the old universality of Banach spaces of density $\aleph_1$ (Theorem 1.5). Moreover, adding one Cohen real adds a weakly compactly generated Banach space of density $\aleph_1$ which does not embed into any Banach space of density $\aleph_1$ in the ground model. We develop the framework of natural spaces to study isomorphic embeddings of Banach spaces and use it to show that a sufficient failure of the generalized continuum hypothesis implies that the universality number of Banach spaces of a given density under a certain kind of positive embeddings (very positive embeddings), is high (Theorem 4.8(1)), and similarly for the number of pairwise non-isomorphic models (Theorem 4.8(2)).


title: Semi-implicit Euler schemes for ordinary differential inclusions

abstract: Two semi-implicit Euler schemes for differential inclusions are proposed and analyzed in depth. An error analysis shows that both semi-implicit schemes inherit favorable stability properties from the differential inclusion. Their performance is considerably better than that of the implicit Euler scheme, because instead of implicit inclusions only implicit equations have to be solved for computing their images. In addition, they are more robust with respect to spatial discretization than the implicit Euler scheme.


title: Local existence and uniqueness for exterior static vacuum Einstein metrics

abstract: We study solutions to the static vacuum Einstein equations on exterior domains with prescribed metric and mean curvature on the inner boundary. It is proved that for any such boundary data near the standard round boundary data in Euclidean space, there exists a unique AF solution to the static vacuum equations realizing the boundary data, which is close to the standard flat solution.


title: Robust boundary tracking for reachable sets of nonlinear differential inclusions

abstract: The Euler scheme is up to date the most important numerical method for ordinary differential inclusions, because the use of the available higher-order methods is prohibited by their enormous complexity after spatial discretization. Therefore, it makes sense to reassess the Euler scheme and optimize its performance. In the present paper, a considerable reduction of the computational cost is achieved by setting up a numerical method that computes the boundaries instead of the complete reachable sets of the fully discretized Euler scheme from lower-dimensional data only. Rigorous proofs for the propriety of this method are given, and numerical examples illustrate the gain of computational efficiency as well as the robustness of the scheme against changes of topology of the reachable sets.


title: Astrophysically robust systematics removal using variational inference: application to the first month of Kepler data

abstract: Space-based transit search missions such as Kepler are collecting large numbers of stellar light curves of unprecedented photometric precision and time coverage. However, before this scientific goldmine can be exploited fully, the data must be cleaned of instrumental artefacts. We present a new method to correct common-mode systematics in large ensembles of very high precision light curves. It is based on a Bayesian linear basis model and uses shrinkage priors for robustness, variational inference for speed, and a de-noising step based on empirical mode decomposition to prevent the introduction of spurious noise into the corrected light curves. After demonstrating the performance of our method on a synthetic dataset, we apply it to the first month of Kepler data. We compare the results, which are publicly available, to the output of the Kepler pipeline's pre-search data conditioning, and show that the two generally give similar results, but the light curves corrected using our approach have lower scatter, on average, on both long and short timescales. We finish by discussing some limitations of our method and outlining some avenues for further development. The trend-corrected data produced by our approach are publicly available.


title: The Suborbital Particle Aggregation and Collision Experiment (SPACE): Studying the Collision Behavior of Submillimeter-Sized Dust Aggregates on the Suborbital Rocket Flight REXUS 12

abstract: The Suborbital Particle Aggregation and Collision Experiment (SPACE) is a novel approach to study the collision properties of submillimeter-sized, highly porous dust aggregates. The experiment was designed, built and carried out to increase our knowledge about the processes dominating the first phase of planet formation. During this phase, the growth of planetary precursors occurs by agglomeration of micrometer-sized dust grains into aggregates of at least millimeters to centimeters in size. However, the formation of larger bodies from the so-formed building blocks is not yet fully understood. Recent numerical models on dust growth lack a particular support by experimental studies in the size range of submillimeters, because these particles are predicted to collide at very gentle relative velocities of below 1 cm/s that can only be achieved in a reduced-gravity environment. The SPACE experiment investigates the collision behavior of an ensemble of silicate-dust aggregates inside several evacuated glass containers which are being agitated by a shaker to induce the desired collisions at chosen velocities. The dust aggregates are being observed by a high-speed camera, allowing for the determination of the collision properties of the protoplanetary dust analog material. The data obtained from the suborbital flight with the REXUS (Rocket Experiments for University Students) 12 rocket will be directly implemented into a state-of-the-art dust growth and collision model.


title: Interpreting Cross-correlations of One-bit Filtered Seismic Noise

abstract: Seismic noise, generated by oceanic microseisms and other sources, illuminates the crust in a manner different from tectonic sources, and therefore provides independent information. The primary measurable is the two-point cross-correlation, evaluated using traces recorded at a pair of seismometers over a finite-time interval. However, raw seismic traces contain intermittent large-amplitude perturbations arising from tectonic activity and instrumental errors, which may corrupt the estimated cross-correlations of microseismic fluctuations. In order to diminish the impact of these perturbations, the recorded traces are filtered using the nonlinear one-bit digitizer, which replaces the measurement by its sign. Previous theory shows that for stationary Gaussian-distributed seismic noise fluctuations one-bit and raw correlation functions are related by a simple invertible transformation. Here we extend this to show that the simple correspondence between these two correlation techniques remains valid for {\it non-stationary} Gaussian and a very broad range of {\it non-Gaussian} processes as well. For a limited range of stationary and non-stationary Gaussian fluctuations, we find that one-bit filtering performs at least as well as spectral whitening. We therefore recommend using one-bit filtering when processing terrestrial seismic noise, with the substantial benefit that the measurements are fully compatible with current theoretical interpretation (e.g., adjoint theory). Given that seismic records are non-stationary and comprise small-amplitude fluctuations and intermittent, large-amplitude tectonic/other perturbations, we outline an algorithm to accurately retrieve the correlation function of the small-amplitude signals.


title: Bifurcation analysis of a smoothed model of a forced impacting beam and comparison with an experiment

abstract: A piecewise-linear model with a single degree of freedom is derived from first principles for a driven vertical cantilever beam with a localized mass and symmetric stops. The resulting piecewise-linear dynamical system is smoothed by a switching function (nonlinear homotopy). For the chosen smoothing function it is shown that the smoothing can induce bifurcations in certain parameter regimes. These induced bifurcations disappear when the transition of the switching is sufficiently and increasingly localized as the impact becomes harder. The bifurcation structure of the impact oscillator response is investigated via the one- and two-parameter continuation of periodic orbits in the driving frequency and/or forcing amplitude. The results are in good agreement with experimental measurements.


title: On the State and Importance of Reproducible Experimental Research in Parallel Computing

abstract: Computer science is also an experimental science. This is particularly the case for parallel computing, which is in a total state of flux, and where experiments are necessary to substantiate, complement, and challenge theoretical modeling and analysis. Here, experimental work is as important as are advances in theory, that are indeed often driven by the experimental findings. In parallel computing, scientific contributions presented in research articles are therefore often based on experimental data, with a substantial part devoted to presenting and discussing the experimental findings. As in all of experimental science, experiments must be presented in a way that makes reproduction by other researchers possible, in principle. Despite appearance to the contrary, we contend that reproducibility plays a small role, and is typically not achieved. As can be found, articles often do not have a sufficiently detailed description of their experiments, and do not make available the software used to obtain the claimed results. As a consequence, parallel computational results are most often impossible to reproduce, often questionable, and therefore of little or no scientific value. We believe that the description of how to reproduce findings should play an important part in every serious, experiment-based parallel computing research article. We aim to initiate a discussion of the reproducibility issue in parallel computing, and elaborate on the importance of reproducible research for (1) better and sounder technical/scientific papers, (2) a sounder and more efficient review process and (3) more effective collective work. This paper expresses our current view on the subject and should be read as a position statement for discussion and future work. We do not consider the related (but no less important) issue of the quality of the experimental design.


title: Logarithmic Sobolev inequalities and spectral concentration for the cubic Schr\"odinger equation

abstract: The nonlinear Schr\"odinger equation NLSE(p, \beta), -iu_t=-u_{xx}+\beta | u|^{p-2} u=0, arises from a Hamiltonian on infinite-dimensional phase space \Lp^2(\mT). For p\leq 6, Bourgain (Comm. Math. Phys. 166 (1994), 1--26) has shown that there exists a Gibbs measure \mu^{\beta}_N on balls \Omega_N= {\phi \in \Lp^2(\mT) : | \phi |^2_{\Lp^2} \leq N} in phase space such that the Cauchy problem for NLSE(p,\beta) is well posed on the support of \mu^{\beta}_N, and that \mu^{\beta}_N is invariant under the flow. This paper shows that \mu^{\beta}_N satisfies a logarithmic Sobolev inequality for the focussing case \beta <0 and 2\leq p\leq 4 on \Omega_N for all N>0; also \mu^{\beta} satisfies a restricted LSI for 4\leq p\leq 6 on compact subsets of \Omega_N determined by H\"older norms. Hence for p=4, the spectral data of the periodic Dirac operator in \Lp^2(\mT; \mC^2) with random potential \phi subject to \mu^{\beta}_N are concentrated near to their mean values. The paper concludes with a similar result for the spectral data of Hill's equation when the potential is random and subject to the Gibbs measure of KdV.


title: Continuum Contributions to the SDO/AIA Passbands During Solar Flares

abstract: Data from the Multiple EUV Grating Spectrograph (MEGS-A) component of the Extreme Ultraviolet Experiment (EVE) onboard the Solar Dynamics Observatory (SDO) were used to quantify the contribution of continuum emission to each of the EUV channels of the Atmospheric Imaging Assembly (AIA), also on SDO, during an X-class solar flare that occurred on 2011 February 15. Both the pre-flare subtracted EVE spectra and fits to the associated free-free continuum were convolved with the AIA response functions of the seven EUV passbands at 10 s cadence throughout the course of the flare. It was found that 10-25% of the total emission in the 94A, 131A, 193A, and 335A passbands throughout the main phase of the flare was due to free-free emission. Reliable measurements could not be made for the 171A channel, while the continuum contribution to the 304A channel was negligible due to the presence of the strong He II emission line. Up to 50% of the emission in the 211A channel was found to be due to free-free emission around the peak of the flare, while an additional 20% was due to the recombination continuum of He II. The analysis was extended to a number of M- and X-class flares and it was found that the level of free-free emission contributing to the 171A and 211A passbands increased with increasing GOES class. These results suggest that the amount of continuum emission that contributes to AIA observations during flares is more significant than that stated in previous studies which used synthetic, rather than observed, spectra. These findings highlight the importance of spectroscopic observations carried out in conjunction with those from imaging instruments so that the data are interpreted correctly.


title: Polyhedrons and PBIBDs from hyperbolic manifolds

abstract: By taking quotients of a certain tiling of hyperbolic plane / space by certain group actions, we obtain geometric polyhedra / cellulations with interesting symmetries and incidence structure.


title: Top quark anomalous tensor couplings in the two-Higgs-doublet models

abstract: We compute the one loop right and left anomalous tensor couplings ($g_R$ and $g_L$, respectively) for the top quark, in the aligned two-Higgs-doublet model. They are the magnetic-like couplings in the most general parameterization of the $tbW$ vertex. We find that the aligned two-Higgs doublet model, that includes as particular cases some of the most studied extensions of the Higgs sector, introduces new electroweak contributions and provides theoretical predictions that are very sensitive to both new scalar masses and the neutral scalar mixing angle. For a large area in the parameters space we obtain significant deviations in both the real and the imaginary parts of the couplings $g_R$ and $g_L$, compared to the predictions given by the electroweak sector of the Standard Model. The most important ones are those involving the imaginary part of the left coupling $g_{L}$ and the real part of the right coupling $g_R$. The real part of $g_L$ and the imaginary part of $g_R$ also show an important sensitivity to new physics scenarios. The model can also account for new CP violation effects via the introduction of complex alignment parameters that have important consequences on the values for the imaginary parts of the couplings. The top anomalous tensor couplings will be measured at the LHC and at future colliders providing a complementary insight on new physics, independent from the bounds in top decays coming from B physics and $b \rightarrow s \gamma$.


title: Low Energy Probes of PeV Scale Sfermions

abstract: We derive bounds on squark and slepton masses in mini-split supersymmetry scenario using low energy experiments. In this setup gauginos are at the TeV scale, while sfermions are heavier by a loop factor. We cover the most sensitive low energy probes including electric dipole moments (EDMs), meson oscillations and charged lepton flavor violation (LFV) transitions. A leading log resummation of the large logs of gluino to sfermion mass ratio is performed. A sensitivity to PeV squark masses is obtained at present from kaon mixing measurements. A number of observables, including neutron EDMs, mu->e transitions and charmed meson mixing, will start probing sfermion masses in the 100 TeV-1000 TeV range with the projected improvements in the experimental sensitivities. We also discuss the implications of our results for a variety of models that address the flavor hierarchy of quarks and leptons. We find that EDM searches will be a robust probe of models in which fermion masses are generated radiatively, while LFV searches remain sensitive to simple-texture based flavor models.


title: Connection Matrices and the Definability of Graph Parameters

abstract: In this paper we extend and prove in detail the Finite Rank Theorem for connection matrices of graph parameters definable in Monadic Second Order Logic with counting (CMSOL) from B. Godlin, T. Kotek and J.A. Makowsky (2008) and J.A. Makowsky (2009). We demonstrate its vast applicability in simplifying known and new non-definability results of graph properties and finding new non-definability results for graph parameters. We also prove a Feferman-Vaught Theorem for the logic CFOL, First Order Logic with the modular counting quantifiers.


title: Models with radiative neutrino masses and viable dark matter candidates

abstract: We provide a list of particle physics models at the TeV-scale that are compatible with neutrino masses and dark matter. In these models, the Standard Model particle content is extended with a small number (\leq 4) of scalar and fermion fields transforming as singlets, doublets or triplets under SU(2), and neutrino masses are generated radiatively via 1-loop diagrams. The dark matter candidates are stabilized by a Z_2 symmetry and are in general mixtures of the neutral components of such new multiplets. We describe the particle content of each of these models and determine the conditions under which they are consistent with current data. We find a total of 35 viable models, most of which have not been previously studied in the literature. There is a great potential to test these models at the LHC not only due to the TeV-scale masses of the new fields but also because about half of the viable models contain particles with exotic electric charges, which give rise to background-free signals. Our results should serve as a first step for detailed analysis of models that can simultaneously account for dark matter and neutrino masses.


title: Brownian Motion of Boomerang Colloidal Particles

abstract: We investigate the Brownian motion of boomerang colloidal particles confined between two glass plates. Our experimental observations show that the mean displacements are biased towards the center of hydrodynamic stress (CoH), and that the mean-square displacements exhibit a crossover from short time faster to long time slower diffusion with the short-time diffusion coefficients dependent on the points used for tracking. A model based on Langevin theory elucidates that these behaviors are ascribed to a superposition of two diffusive modes: the ellipsoidal motion of the CoH and the rotational motion of the tracking point with respect to the CoH.


title: Hoodsquare: Modeling and Recommending Neighborhoods in Location-based Social Networks

abstract: Information garnered from activity on location-based social networks can be harnessed to characterize urban spaces and organize them into neighborhoods. In this work, we adopt a data-driven approach to the identification and modeling of urban neighborhoods using location-based social networks. We represent geographic points in the city using spatio-temporal information about Foursquare user check-ins and semantic information about places, with the goal of developing features to input into a novel neighborhood detection algorithm. The algorithm first employs a similarity metric that assesses the homogeneity of a geographic area, and then with a simple mechanism of geographic navigation, it detects the boundaries of a city's neighborhoods. The models and algorithms devised are subsequently integrated into a publicly available, map-based tool named Hoodsquare that allows users to explore activities and neighborhoods in cities around the world. Finally, we evaluate Hoodsquare in the context of a recommendation application where user profiles are matched to urban neighborhoods. By comparing with a number of baselines, we demonstrate how Hoodsquare can be used to accurately predict the home neighborhood of Twitter users. We also show that we are able to suggest neighborhoods geographically constrained in size, a desirable property in mobile recommendation scenarios for which geographical precision is key.


title: The Effects of Dissolved Methane upon Liquid Argon Scintillation Light

abstract: In this paper we report on measurements of the effects of dissolved methane upon argon scintillation light. We monitor the light yield from an alpha source held 20 cm from a cryogenic photomultiplier tube (PMT) assembly as methane is injected into a high-purity liquid argon volume. We observe significant suppression of the scintillation light yield by dissolved methane at the 10 part per billion (ppb) level. By examining the late scintillation light time constant, we determine that this loss is caused by an absorption process and also see some evidence of methane-induced scintillation quenching at higher concentrations (50-100 ppb). Using a second PMT assembly we look for visible re-emission features from the dissolved methane which have been reported in gas-phase argon methane mixtures, and we find no evidence of visible re-emission from liquid-phase argon methane mixtures at concentrations between 10 ppb and 0.1%.


title: Area inequalities for stable marginally outer trapped surfaces in Einstein-Maxwell-dilaton theory

abstract: We prove area inequalities for stable marginally outer trapped surfaces in Einstein-Maxwell-dilaton theory. Our inspiration comes on the one hand from a corresponding upper bound for the area in terms of the charges obtained recently by Dain, Jaramillo and Reiris [1] in the pure Einstein-Maxwell case without symmetries, and on the other hand from Yazadjiev's inequality [2] in the axially symmetric Einstein-Maxwell-dilaton case. The common issue in these proofs and in the present one is a functional ${\mathscr W}$ of the matter fields for which the stability condition readily yields an {\it upper} bound. On the other hand, the step which crucially depends on whether or not a dilaton field is present is to obtain a {\it lower} bound for ${\mathscr W}$ as well. We obtain the latter by first setting up a variational principle for ${\mathscr W}$ with respect to the dilaton field $\phi$, then by proving existence of a minimizer $\psi$ as solution of the corresponding Euler-Lagrange equations and finally by estimating ${\mathscr W}(\psi)$. In the special case that the normal components of the electric and magnetic fields are proportional we obtain the area bound $A \ge 8\pi P Q$ in terms of the electric and magnetic charges. In the generic case our results are less explicit but imply rigorous `perturbation' results for the above inequality. All our inequalities are saturated for a 2-parameter family of static, extreme solutions found by Gibbons [3]. Via the Bekenstein-Hawking relation $A = 4S$ our results give positive lower bounds for the entropy $S$ which are particularly interesting in the Einstein-Maxwell-dilaton case.


title: Precise QCD predictions for the production of a photon pair in association with two jets

abstract: We compute the cross section for the production of a high-mass photon pair in association with two hadronic jets to next-to-leading order (NLO) in quantum chromodynamics (QCD). Our results allow for the first time to reliably predict the absolute normalisation of this process, and demonstrate that the shape of important kinematical distributions is modified by higher-order effects. The perturbative corrections will be an important ingredient to precision studies of Higgs boson properties from its production in association with two jets.


title: Giant vortex states in type I superconductors simulated by Ginzburg-Landau equations

abstract: The quantization of magnetic flux in superconductors is usually seen as vortices penetrating the sample. While vortices are unstable in bulk type I superconductors, restricting the superconductor causes a variety of vortex structures to appear. We present a systematic study of giant vortex states in type I superconductors obtained by numerically solving the Ginzburg-Landau equations. The size of the vortices is seen to increase with decreasing film thickness. In type I superconductors, giant vortices appear at intermediate thicknesses but they do not form a well-defined vortex lattice. In the thinnest type I films, singly quantized vortices seem to be stabilized by the geometry of the sample instead of an increase in the effective Ginzburg-Landau parameter.


title: A Convex Framework for Optimal Investment on Disease Awareness in Social Networks

abstract: We consider the problem of controlling the propagation of an epidemic outbreak in an arbitrary network of contacts by investing on disease awareness throughout the network. We model the effect of agent awareness on the dynamics of an epidemic using the SAIS epidemic model, an extension of the SIS epidemic model that includes a state of "awareness". This model allows to derive a condition to control the spread of an epidemic outbreak in terms of the eigenvalues of a matrix that depends on the network structure and the parameters of the model. We study the problem of finding the cost-optimal investment on disease awareness throughout the network when the cost function presents some realistic properties. We propose a convex framework to find cost-optimal allocation of resources. We validate our results with numerical simulations in a real online social network.


title: The role of blood circulatory system in thermal regulation of animals explained by entropy production analysis

abstract: A novel model of thermal regulation of homoeothermic animals has been implemented. The model is based on a non-equilibrium thermodynamic approach which introduces entropy balance and the rate of entropy generation as a formulation of The Second Law. The model explains linear correlation between an animals skin and environment temperatures using the first principles and demonstrates the role of blood circulation in the thermoregulation of homoeothermic animals.


title: Constraints on the Skyrme Equations of State from Properties of Doubly Magic Nuclei

abstract: I use properties of doubly-magic nuclei to constrain nuclear matter and neutron matter equations of state. I conclude that the data determined the value of the neutron equation of state and the symmetry energy near a density of $\rho_{on}$ = 0.10 nucleons/fm$^{3}$. The slope at that point is constrained by the value of the neutron skin.


title: On Sparsification for Computing Treewidth

abstract: We investigate whether an n-vertex instance (G,k) of Treewidth, asking whether the graph G has treewidth at most k, can efficiently be made sparse without changing its answer. By giving a special form of OR-cross-composition, we prove that this is unlikely: if there is an e > 0 and a polynomial-time algorithm that reduces n-vertex Treewidth instances to equivalent instances, of an arbitrary problem, with O(n^{2-e}) bits, then NP is in coNP/poly and the polynomial hierarchy collapses to its third level. Our sparsification lower bound has implications for structural parameterizations of Treewidth: parameterizations by measures that do not exceed the vertex count, cannot have kernels with O(k^{2-e}) bits for any e > 0, unless NP is in coNP/poly. Motivated by the question of determining the optimal kernel size for Treewidth parameterized by vertex cover, we improve the O(k^3)-vertex kernel from Bodlaender et al. (STACS 2011) to a kernel with O(k^2) vertices. Our improved kernel is based on a novel form of treewidth-invariant set. We use the q-expansion lemma of Fomin et al. (STACS 2011) to find such sets efficiently in graphs whose vertex count is superquadratic in their vertex cover number.


title: The Nuclear Pairing Gap -- How Low Can It Go?

abstract: The pairing gap for $^{53}$Ca obtained from new experimental data on the masses of $^{52-54}$Ca has the smallest value yet observed. This is explained in the framework of the nuclear shell model with schematic and realistic Hamiltonians as being due to shell gaps around the low-$ j $ orbital $ 1p_{1/2} $. Minima in the pairing gaps for all nuclei are shown and discussed


title: Chemical energy in an introductory physics course for the life sciences

abstract: Energy is a complex idea that cuts across scientific disciplines. For life science students, an approach to energy that incorporates chemical bonds and chemical reactions is better equipped to meet the needs of life sciences students than a traditional introductory physics approach that focuses primarily on mechanical energy. We present a curricular sequence, or thread, designed to build up students' understanding of chemical energy in an introductory physics course for the life sciences. This thread is designed to connect ideas about energy from physics, biology, and chemistry. We describe the kinds of connections among energetic concepts that we intended to develop to build interdisciplinary coherence, and present some examples of curriculum materials and student data that illustrate our approach.


title: Kinetic properties in inhomogeneous self-aware media

abstract: The new framework for finance is proposed. This framework based on three known approaches in econophysics. Assumptions of the framework are the following: 1. For the majority of situations market follows non-arbitrage condition. 2. For the small number of situations market influenced by the actions of big firms. 3. If actions of big players lead to the arbitrage opportunity, small players could self-organize to take advantage of this opportunity. The framework is an attempt to combine approaches of Bouchaud, Gabaix, Sornette, Stanley and coauthors. Suggested framework is applied for the analysis of market impact models, behavior of big players, self-organization of market firm and volatility description.


title: A generalisation of core partitions

abstract: Suppose $s$ and $t$ are coprime natural numbers. A theorem of Olsson says that the $t$-core of an $s$-core partition is again an $s$-core. We generalise this theorem, showing that the $s$-weight of the $t$-core of a partition $\lambda$ is at most the $s$-weight of $\lambda$. Then we consider the set $\mathcal C_{s:t}$ of partitions for which equality holds, which we call $[s{:}t]$-cores; this set has interesting structure, and we expect that it will be the subject of future study. We show that the set of $[s{:}t]$-cores is a union of finitely many orbits for an action of a Coxeter group of type $\tilde A_{s-1}\times\tilde A_{t-1}$ on the set of partitions. We also consider the problem of constructing an $[s{:}t]$-core with specified $s$-core and $t$-core.


title: Tight Lagrangian homology spheres in compact homogeneous K\"ahler manifolds

abstract: For any irreducible compact homogeneous K\"ahler manifold, we classify the compact tight Lagrangian submanifolds which have the Z_2-homology of a sphere.


title: THz generation by optical Cherenkov emission from ionizing two-color laser pulses

abstract: Two-color photoionization produces a cycle-averaged current driving broadband, conically emitted THz radiation. We investigate, through simulation, the processes determining the angle of conical emission. We find that the emission angle is determined by an optical Cherenkov effect, where the front velocity of the current source is faster than the THz phase velocity.


title: Emergent spin

abstract: Quantum mechanics and relativity in the continuum imply the well known spin-statistics connection. However for particles hopping on a lattice, there is no such constraint. If a lattice model yields a relativistic field theory in a continuum limit, this constraint must "emerge" for physical excitations. We discuss a few models where a spin-less fermion hopping on a lattice gives excitations which satisfy the continuum Dirac equation. This includes such well known systems such as graphene and staggered fermions.


title: Substance Abuse via Legally Prescribed Drugs: The Case of Vicodin in the United States

abstract: Vicodin is the most commonly prescribed pain reliever in the United States. Research indicates that there are two million people who are currently abusing Vicodin, and the majority of those who abuse Vicodin were initially exposed to it via prescription. Our goal is to determine the most effective strategies for reducing the overall population of Vicodin abusers. More specifically, we focus on whether prevention methods aimed at educating doctors and patients on the potential for drug abuse or treatment methods implemented after a person abuses Vicodin will have a greater overall impact. We consider one linear and two non-linear compartmental models in which medical users of Vicodin can transition into the abuser compartment or leave the population by no longer taking the drug. Once Vicodin abusers, people can transition into a treatment compartment, with the possibility of leaving the population through successful completion of treatment or of relapsing and re-entering the abusive compartment. The linear model assumes no social interaction, while both non-linear models consider interaction. One considers interaction with abusers affecting the relapse rate, while the other assumes both this and an additional interaction between the number of abusers and the number of new prescriptions. Sensitivity analyses are conducted varying the rates of success of these intervention methods measured by the parameters to determine which strategy has the greatest impact on controlling the population of Vicodin abusers. From these models and analyses, we determine that manipulating parameters tied to prevention measures has a greater impact on reducing the population of abusers than manipulating parameters associated with treatment. We also note that increasing the rate at which abusers seek treatment affects the population of abusers more than the success rate of treatment itself.


title: Black hole triple dynamics: breakdown of the orbit average approximation and implications for gravitational wave detections

abstract: Coalescing black hole (BH) binaries forming in the dense core of globular clusters (GCs) are expected to be one the brightest sources of gravitational wave (GW) radiation for the next generation of ground-based laser interferometers. Favorable conditions for merger are initiated by the Kozai resonance in which the gravitational interaction with a third distant object, typically another BH, induces quasi-periodic variations of the inner BH binary eccentricity. In this paper we perform high precision N-body simulations of the long term evolution of hierarchical BH triples and investigate the conditions that lead to the merging of the BH binary and the way it might become an observable source of GW radiation. We find that the secular orbit average treatment, adopted in previous works, does not reliably describe the dynamics of these systems if the binary is orbited by the outer BH on a highly inclined orbit at a moderate distance. During the high eccentricity phase of a Kozai cycle the torque due to the outer BH can drive the binary to extremely large eccentricities in a fraction of the binary's orbital period. This occurs before relativistic terms become important to the evolution and allows the binary GW signal to reach large GW frequencies (>~10 Hz) at high eccentricities. We show that ~50 % of coalescing BH binaries driven by the Kozai mechanism in GCs will have eccentricities larger than 0.1, with 10 % of them being extremely eccentric, (1-e)<~10^-5, when they first chirp in the frequency band of ground based laser interferometers. This implies that a large fraction of such GW sources could be missed if conventional quasi-circular templates are used for analysis of GW detectors data. The efficient detection of all coalescing BH binaries in GCs will therefore require template banks of eccentric inspiral waveforms for matched-filtering and dedicated search strategies.


title: Dynamics and Control of an Invasive Species: The Case of the Rasberry Crazy Ant Colonies

abstract: This project is motivated by the costs related with the documented risks of the introduction of non-native invasive species of plants, animals, or pathogens associated with travel and international trade. Such invasive species often have no natural enemies in their new regions. The spatiotemporal dynamics related to the invasion/spread of Nylanderia fulva, commonly known as the Rasberry crazy ant, are explored via the use of models that focus on the reproduction of ant colonies. A Cellular Automaton (CA) simulates the spatially explicit spread of ants on a grid. The impact of local spatial correlations on the dynamics of invasion is investigated numerically and analytically with the aid of a Mean Field (MF) model and a Pair Approximation (PA) model, the latter of which accounts for adjacent cell level efects. The PA model approach considers the limited mobility range of N. fulva, that is, the grid cell dynamics are not strongly in uenced by non-adjacent cells. The model determines the rate of growth of colonies of N. fulva under distinct cell spatial architecture. Numerical results and qualitative conclusions on the spread and control of this invasive ant species are discussed.


title: Thermodynamic anomalies and structural fluctuations in aqueous solutions of tertiary butyl alcohol

abstract: In this work, we discuss the connection between the anomalies of the thermodynamic properties, experimentally observed in tertiary butyl alcohol (TBA) and water solutions, and the molecular clustering in these solutions, as revealed by molecular dynamics (MD) simulations. These anomalies are observed in relatively dilute solutions of about 0.03 to 0.08 mole fraction of TBA and become more pronounced at low temperatures. MD simulations show that these solutions exhibit shortranged (order of 1 nm), shortlived (tens of picoseconds) "micelle-like" structural fluctuations in the same concentration range. We attribute the anomalies in the thermodynamic properties of aqueous TBA solutions to these structural fluctuations on the molecular scale.


title: Neutron-hole states in 45Ar from p(46Ar,d)45Ar reactions

abstract: To improve the effective interactions in the pf shell, it is important to measure the single particle- and hole- states near the N=28 shell gap. In this paper, the neutron spectroscopic factors of hole-states from the unstable neutron-rich 45Ar (Z=18, N=27) nucleus have been studied using 1H(46Ar, 2H)45Ar transfer reaction in inverse kinematics. Comparison of our results with the particle-states of 45Ar produced in 2H(44Ar, H)45Ar reaction shows that the two reactions populate states with different angular momentum. Using the angular distributions, we are able to confirm the spin assignments of four low-lying states of 45Ar. These are the ground state (f7/2), the first-excited (p3/2), the s1/2 and the d3/2 states. While large basis shell model predictions describe spectroscopic properties of the ground and p3/2 states very well, they fail to describe the s1/2 and d3/2 hole-states.


title: Preventing Exceptions to Robins InEquality

abstract: For sufficiently large n Ramanujan gave a sufficient condition for the truth Robin's InEquality $X(n):=\frac{\sigma(n)}{n\ln\ln n}<e^{\gamma}$ (RIE). The largest known violation of RIE is $n_8=5040$. In this paper Robin's multipliers are split into logarithmic terms $\mathcal{L}$ and relative divisor sums $\mathcal{G}$. A violation of RIE above $n_{8}$ is proposed to imply oscillations that cause $\mathcal{G}$ to exceed $\mathcal{L}$. To this aim Alaoglu and Erd\H{o}s's conjecture for the CA numbers algorithm is used and the paper's key points are in section 4.2


title: Just In Time Indexing

abstract: One of the major challenges being faced by Database managers today is to manage the performance of complex SQL queries which are dynamic in nature. Since it is not possible to tune each and every query because of its dynamic nature, there is a definite possibility that these queries may cause serious database performance issues if left alone. Conventional indexes are useful only for those queries which are frequently executed or those columns which are frequently joined in SQL queries. This proposal is regarding a method, a query optimizer for optimizing database queries in a database management system. Just In Time(JIT) indexes are On Demand, temporary indexes created on the fly based on current needs so that they would be able to satisfy any kind of queries. JIT indexes are created only when the configured threshold values for resource consumption are exceeded for a query. JIT indexes will be stored in a temporary basis and will get replaced by new JIT indexes in course of time. The proposal is substantiated with the help of experimental programs and with various test cases. The idea of parallel programming is also brought into picture as it can be effectively used in a multiprocessor system. Multiple threads are employed while one set of threads proceed in the conventional way and the other set of threads proceed with the proposed way. A live switch over is made when a suitable stage is reached and from then onwards the proposed method will only come into picture.


title: Tilted Lema\^itre model and the dark flow

abstract: In the last years, the peculiar velocities of many X-ray galaxies clusters with respect to the distance have been measured directly in the rest frame of the cosmic microwave background radiation (CBR), using the kinematic Sunyaev-Zeldovich (kSZ) effect. These measures prove that exists a highly coherent motion, extending out to at least to ~1 Gpc, of the matter rest frame with respect to the CBR rest frame. This global motion was named "dark flow". By using an inhomogeneous spherically symmetric "tilted" Lema\^itre model, we could explain the dark flow if we assume a linear increase with distance of the peculiar velocities, which is in principle allowed by these observations. This linear increase of the dark flow with the distance has the same behavior that the intrinsic dipole, due to the kinematic acceleration, which appears in the Hubble law of the Lema\^itre model. In the "tilted" Lema\^itre model considered, we consider that the radiation orthogonal congruence is a perfect fluid and the matter "tilted" congruence is an imperfect fluid with heat flux.


title: Prediction Capabilities of VLF/LF Emission as the Main Precursor of Earthquake

abstract: Recent satellite and ground-based observations proved that in earthquake preparation period in the seismogenic area we have VLF/LF and ULF electromagnetic emissions. According to the opinion of the authors of the present paper this phenomenon is more universal and reliable than other earthquake indicators. Hypothetically, in case of availability of adequate methodological grounds, in the nearest future, earth VLF/LF electromagnetic emission might be declared as the main precursor of earthquake. In particular, permanent monitoring of frequency spectrum of earth electromagnetic emission generated in the earthquake preparation period might turn out very useful with the view of prediction of large (M 5) inland earthquakes. The present paper offers a scheme of the methodology according to which the reality of the above given hypothesis can be checked up. To prove the prediction capabilities of earth electromagnetic emission we have used avalanche-like unstable model of fault formation and an analogous model of electromagnetic contour, synthesis of which, according to our opinion, is rather harmonious.


title: A Closer Look at the Fluctuations in Brightness of SN 2009ip During Its Late 2012 Eruption

abstract: The supernova impostor SN 2009ip has re-brightened several times since its initial discovery in August 2009. During its last outburst in late September 2012 it reached a peak brightness of m$_v$ $\sim$ 13.5 (M$_v$ brighter than -18) causing some to speculate that it had undergone a terminal core-collapse supernova. Relatively high-cadence multi-wavelength photometry of the post-peak decline revealed bumps in brightness infrequently observed in other Type IIn supernovae. These bumps occurred synchronously in all UV and optical bands with amplitudes of 0.1 -- 0.4 mag at intervals of 10 -- 30 days. Episodic continuum brightening and dimming in the UV and optical with these characteristics is not easilly explained within the context of models that have been proposed for the late September 2012 outburst of SN 2009ip. We also present evidence that the post peak fluctuations in brightness occur at regular intervals and raise more questions about their origin.


title: Non-Zeeman Circular Polarization of CO rotational lines in SNR IC 443

abstract: Context: We investigate non-Zeeman circular polarization and linear polarization levels of up to 1% of $^{12}$CO spectral line emission detected in a shocked molecular clump around the supernova remnant (SNR) IC 443, with the goal of understanding the magnetic field structure in this source. Aims: We examine our polarization results to confirm that the circular polarization signal in CO lines is caused by a conversion of linear to circular polarization, consistent with anisotropic resonant scattering. In this process background linearly polarized CO emission interacts with similar foreground molecules aligned with the ambient magnetic field and scatters at a transition frequency. The difference in phase shift between the orthogonally polarized components of this scattered emission can cause a transformation of linear to circular polarization. Methods: We compared linear polarization maps from dust continuum, obtained with PolKa at APEX, and $^{12}$CO ($J=2\rightarrow1$) and ($J=1\rightarrow0$) from the IRAM 30-m telescope and found no consistency between the two sets of polarization maps. We then reinserted the measured circular polarization signal in the CO lines across the source to the corresponding linear polarization signal to test whether before this linear to circular polarization conversion the linear polarization vectors of the CO maps were aligned with those of the dust. Results: After the flux correction for the two transitions of the CO spectral lines, the new polarization vectors for both CO transitions aligned with the dust polarization vectors, establishing that the non-Zeeman CO circular polarization is due to a linear to circular polarization conversion.


title: Dielectric relaxation of water below the melting point: the effect of inner pressure

abstract: Despite water is the most studied substance in the Earth, it is not completely understood why its structural and dynamical properties give rise to some anomalous behaviors. Interesting properties emerge when experiments at low temperatures and/or high pressures, are performed. Here we report dielectric measurements of cold water under constrained conditions, i.e. water that below the melting point can not freeze. The inner pressure shifts the {\alpha} relaxation peak to similar frequencies as seen in ice Ih. Also, when we reach the triple point at 251 K, ice III seems to form. As far as we know, this via to obtain such crystalline phase has not been observed.


title: The effect of magnetic fields on the r-modes of slowly rotating relativistic neutron stars

abstract: We study here the r-modes in the Cowling approximation of a slowly rotating and magnetized neutron star with a poloidal magnetic field, where we neglect any deformations of the spherical symmetry of the star. We were able to quantify the influence of the magnetic field in both the oscillation frequency $\sigma_r$ of the r-modes and the growth time $t_{gw}$ of the gravitational radiation emission. We conclude that magnetic fields of the order $10^{15}$ G at the center of the star are necessary to produce any changes. Our results for $\sigma_r$ show a decrease of up to $\sim$ $5\%$ in the frequency with increasing magnetic field, with a $B^2$ dependence for rotation rates $\Omega/\Omega_K \gtrsim 0.07$ and $B^4$ for $\Omega/\Omega_K \lesssim 0.07$. (These results should be trusted only within slow rotation approximation and we kept $\Omega/\Omega_{K}< 0.3$.) For $t_{gw}$, we find that it is approximately $30\%$ smaller than previous Newtonian results for non-magnetized stars, which would mean a faster growth of the emission of gravitational radiation. The effect of the magnetic field in $t_{gw}$ causes a non-monotonic effect, that first slightly increases $t_{gw}$ and then decreases it further by another $\sim$ $5\%$. (The value of magnetic field for which $t_{gw}$ starts to decrease depends on the rotational frequency, but it is generally around $10^{15}$G.) Future work should be dedicated to the study of the effect of viscosity in the presence of magnetic fields, in order to establish the magnetic correction to the instability window.


title: A new filtration of the Magnus kernel of the Torelli group

abstract: For a oriented genus g surface with one boundary component, S, the Torelli group is the group of orientation preserving homeomorphisms of S that induce the identity on homology. The Magnus representation of the Torelli group represents the action on F/F" where F=pi_1(S) and F" is the second term of the derived series. We show that the kernel of the Magnus representation, Mag(S), is highly non-trivial and has a rich structure as a group. Specifically, we define an infinite filtration of Mag(S) by subgroups, called the higher order Magnus subgroups, M_k(S). We develop methods for generating nontrivial mapping classes in M_k(S) for all k and g>1. We show that for each k the quotient M_k(S)/M_{k+1}(S) contains a subgroup isomorphic to a lower central series quotient of free groups E(g-1)_k/E(g-1)_{k+1}. Finally We show that for g>2 the quotient M_k(S)/M_{k+1}(S) surjects onto an infinite rank torsion free abelian group. To do this, we define a Johnson-type homomorphism on each higher order Magnus subgroup quotient and show it has a highly non-trivial image.


title: A semiclassical collective response of heated, asymmetric and rotating nuclei

abstract: The Landau Fermi-liquid and extended Gutzwiller periodic-orbit theories are presented for the semiclassical description of collective excitations in nuclei, which are close to main topics of the fruitful activity of S.T. Belyaev. Static susceptibilities show the ergodicity of Fermi liquids. Transport coefficients (nuclear friction and inertia) as functions of the temperature for the hydrodynamic and zero-sound modes are derived within the response theory by using the Fermi-liquid droplet model, in agreement with the shell model for large temperatures. The surface symmetry binding-energy constants are obtained as functions of Skyrme force parameters in the approximation of a sharp edged proton-neutron asymmetric nucleus.The energies and sum rules of the isovector dipole giant resonances are in fairly good agreement with the experimental data. An analysis of the specific structure of these resonances in terms of a main, and satellite peaks, in comparison with the experimental data and microscopic theoretical models, might turn out to be of importance for a better understanding of the values of the surface symmetry-energy constant. The semiclassical collective moment of inertia is derived analytically beyond the quantum perturbation approximation of the cranking model for any potential well as a mean field. It is shown that this moment of inertia can be approximated by its rigid-body value for the rotation with a given frequency within the ETF and more general periodic orbit theories in the nearly local long-length approximation. Its semiclassical shell-structure components are derived in terms of the periodic-orbit free-energy shell corrections. We obtained good agreement between the semiclassical and quantum shell-structure components of the moment of inertia for several critical bifurcation deformations for the harmonic oscillator mean field.


title: Stable reflexive sheaves and localization

abstract: We study moduli spaces $\mathcal{N}$ of rank 2 stable reflexive sheaves on $\mathbb{P}^3$. Fixing Chern classes $c_1$, $c_2$, and summing over $c_3$, we consider the generating function $\mathsf{Z}^{\mathrm{refl}}(q)$ of Euler characteristics of such moduli spaces. The action of the torus $T$ on $\mathbb{P}^3$ lifts to $\mathcal{N}$ and we classify all sheaves in $\mathcal{N}^T$. This leads to an explicit expression for $\mathsf{Z}^{\mathrm{refl}}(q)$. Since $c_3$ is bounded below and above, $\mathsf{Z}^{\mathrm{refl}}(q)$ is a polynomial. We find a simple formula for its leading term when $c_1=-1$. Next, we study moduli spaces of rank 2 stable torsion free sheaves on $\mathbb{P}^3$ and consider the generating function of Euler characteristics of such moduli spaces. We give an expression for this generating function in terms of $\mathsf{Z}^{\mathrm{refl}}(q)$ and Euler characteristics of Quot schemes of certain $T$-equivariant reflexive sheaves, which are studied elsewhere. Many techniques of this paper apply to any toric 3-fold. In general, $\mathsf{Z}^{\mathrm{refl}}(q)$ depends on the choice of polarization which leads to wall-crossing phenomena. We briefly illustrate this in the case of $\mathbb{P}^2 \times \mathbb{P}^1$.


title: Evolving a Behavioral Repertoire for a Walking Robot

abstract: Numerous algorithms have been proposed to allow legged robots to learn to walk. However, the vast majority of these algorithms is devised to learn to walk in a straight line, which is not sufficient to accomplish any real-world mission. Here we introduce the Transferability-based Behavioral Repertoire Evolution algorithm (TBR-Evolution), a novel evolutionary algorithm that simultaneously discovers several hundreds of simple walking controllers, one for each possible direction. By taking advantage of solutions that are usually discarded by evolutionary processes, TBR-Evolution is substantially faster than independently evolving each controller. Our technique relies on two methods: (1) novelty search with local competition, which searches for both high-performing and diverse solutions, and (2) the transferability approach, which com-bines simulations and real tests to evolve controllers for a physical robot. We evaluate this new technique on a hexapod robot. Results show that with only a few dozen short experiments performed on the robot, the algorithm learns a repertoire of con-trollers that allows the robot to reach every point in its reachable space. Overall, TBR-Evolution opens a new kind of learning algorithm that simultaneously optimizes all the achievable behaviors of a robot.


title: Reduction of breast cancer relapses with perioperative non-steroidal anti-inflammatory drugs: new findings and a review

abstract: To explain a bimodal pattern of hazard of relapse among early stage breast cancer patients identified in multiple databases, we proposed that late relapses result from steady stochastic progressions from single dormant malignant cells to avascular micrometastases and then on to growing deposits. However in order to explain early relapses, we had to postulate that something happens at about the time of surgery to provoke sudden exits from dormant phases to active growth and then to detection. Most relapses in breast cancer are in the early category. Recent data from Forget et al suggests an unexpected mechanism. They retrospectively studied results from 327 consecutive breast cancer patients comparing various perioperative analgesics and anesthetics in one Belgian hospital and one surgeon. Patients were treated with mastectomy and conventional adjuvant therapy. Relapse hazard updated Sept 2011 are presented. A common Non-Steroidal Anti-Inflammatory Drug (NSAID) analgesic used in surgery produced far superior disease-free survival in the first 5 years after surgery. The expected prominent early relapse events in months 9-18 are reduced 5-fold. If this observation holds up to further scrutiny, it could mean that the simple use of this safe, inexpensive and effective anti-inflammatory agent at surgery might eliminate early relapses. Transient systemic inflammation accompanying surgery could facilitate angiogenesis of dormant micrometastases, proliferation of dormant single cells, and seeding of circulating cancer stem cells (perhaps in part released from bone marrow) resulting in early relapse and could have been effectively blocked by the perioperative anti-inflammatory agent.


title: Surprisingly Long Length Scales for Semiclassical Loop Quantum Gravity and Their Physical Consequences

abstract: When gauge field theory coherent states for loop quantum gravity (LQG) were introduced, an optimized semiclassical proper length emerged, corresponding to the edge length $\epsilon$ of a graph embedded in a given classical geometry. Here $\epsilon$ is explored in more detail. $\epsilon$ at the Earth's surface is found to lie between 100 $\mu $m and 0.7 m. The implied quantum fluctuating space-time strain amplitude and noise spectrum are estimated to be $4\; 1/2$ orders smaller than the current experimental detectability. However, such a macroscopic $\epsilon $ makes regularization of the semiclassical electromagnetic Hamiltonian problematic for photon wavelengths shorter than $\epsilon$. The origin of a large $\epsilon$ is traced to an edge-wise tensor product of independent edge-based coherent states for the whole graph state. This provides physical grounds for recently proposed collective coherent states, where $\epsilon$ acquires the interpretation of a sliding scale. A new proper distance $\xi$ emerges as the characteristic length of semiclassical LQG. $\xi$ will affect the LQG photon vacuum dispersion relations, and is also accessible to current measurements of space-time strain. Matter interactions may also affect $\xi$.


title: Monte Carlo study of the XY-model on Sierpnski carpet

abstract: We have performed a Monte Carlo study of the classical XY-model on a Sierpi\' nski carpet, which is a planar fractal structure with infinite order of ramification and fractal dimension 1.8928. We employed the Wolff cluster algorithm in our simulations and our results, in particular those for the susceptibility and the helicity modulus, indicate the absence of finite-temperature Berezinskii-Kosterlitz-Thouless (BKT) transition in this system.


title: Business and social evaluation of denial of service attacks of communications networks in view of scaling economic counter-measures

abstract: This paper gives an analytical method to determine the economic and indirect implications of denial of service and distributed denial of service attacks. It is based on time preference dynamics applied to the monetary mass for the restoration of capabilities, on long term investments to rebuild capabilities, and of the usability level of the capabilities after an attack. A simple illustrative example is provided for a denial of service on a corporate data centre. The needed data collection methodologies are categorized by classes of targets. The use of the method is explained in the context of legal or policy driven dissuasive, retaliation or compensation/ restoration actions. A concrete set of deployment cases in mobile communications services is discussed. The conclusion includes policy recommendations as well as information exchange requirements.


title: Predicting the $B -> \rho$ form factors using AdS/QCD Distribution Amplitudes for the $\rho$ meson

abstract: We use QCD light-cone sum rules (LCSR) with holographic anti de Sitter/Chromodynamics (AdS/QCD) Distribution Amplitudes (DAs) for the $\rho$ meson to predict the form factors that govern the leading amplitudes of the rare radiative $B -> \rho \gamma$ decay and the semileptonic decay $B -> \rho l \nu$. We test the Isgur-Wise relation between the radiative and semileptonic form factors. We also compute the total width (in units of $|V_{ub}|^2$) for the semileptonic decay $B -> \rho l \nu$ as well as ratios of partial decay widths which are independent of $|V_{ub}|$.


title: Holographic EPR Pairs, Wormholes and Radiation

abstract: As evidence for the ER=EPR conjecture, it has recently been observed that the string that is holographically dual to an entangled quark-antiquark pair separating with (asymptotically) uniform acceleration has a wormhole on its worldsheet. We point out that a two-sided horizon and a wormhole actually appear for much more generic quark-antiquark trajectories, which is consistent with the fact that the members of an EPR pair need not be permanently out of causal contact. The feature that determines whether the causal structure of the string worldsheet is trivial or not turns out to be the emission of gluonic radiation by the dual quark and antiquark. In the strongly-coupled gauge theory, it is only when radiation is emitted that one obtains an unambiguous separation of the pair into entangled subsystems, and this is what is reflected on the gravity side by the existence of the worldsheet horizon.


title: Universal dynamics of a degenerate unitary Bose gas

abstract: Understanding the rich behavior that emerges from systems of interacting quantum particles, such as electrons in materials, nucleons in nuclei or neutron stars, the quark-gluon plasma, and superfluid liquid helium, requires investigation of systems that are clean, accessible, and have tunable parameters. Ultracold quantum gases offer tremendous promise for this application largely due to an unprecedented control over interactions. Specifically, $a$, the two-body scattering length that characterizes the interaction strength, can be tuned to any value. This offers prospects for experimental access to regimes where the behavior is not well understood because interactions are strong, atom-atom correlations are important, mean-field theory is inadequate, and equilibrium may not be reached or perhaps does not even exist. Of particular interest is the unitary gas, where $a$ is infinite, and where many aspects of the system are universal in that they depend only on the particle density and quantum statistics. While the unitary Fermi gas has been the subject of intense experimental and theoretical investigation, the degenerate unitary Bose gas has generally been deemed experimentally inaccessible because of three-body loss rates that increase dramatically with increasing $a$. Here, we investigate dynamics of a unitary Bose gas for timescales that are short compared to the loss. We find that the momentum distribution of the unitary Bose gas evolves on timescales fast compared to losses, and that both the timescale for this evolution and the limiting shape of the momentum distribution are consistent with universal scaling with density. This work demonstrates that a unitary Bose gas can be created and probed dynamically, and thus opens the door for further exploration of this novel strongly interacting quantum liquid.


title: Frozen Vacuum

abstract: Modes just outside the horizon of a typical old black hole are thermally entangled with distant Hawking radiation. This precludes their entangled purity with interior modes, leading to a firewall. Identifying the interior with the distant radiation ("A=R_B", "ER=EPR") can resolve the entanglement conflict. But the map must adjust for any interactions, or else the firewall will reappear if the Hawking radiation scatters off the CMB. With a self-correcting map, an infalling observer is unable to excite the vacuum near the horizon. This allows the horizon to be locally detected and so violates the equivalence principle.


title: Impact of Short-Range Scattering on the Metallic Transport of Strongly Correlated 2D Holes in GaAs Quantum Wells

abstract: Understanding the non-monotonic behavior in the temperature dependent resistance, R(T), of strongly correlated two-dimensional (2D) carriers in clean semiconductors has been a central issue in the studies of 2D metallic states and metal-insulator-transitions. We have studied the transport of high mobility 2D holes in 20nm wide GaAs quantum wells (QWs) with varying short-range disorder strength by changing the Al fraction x in the Al_xGa_{1-x}As barrier. Via varying the short range interface roughness and alloy scattering, it is observed that increasing x suppresses both the strength and characteristic temperature scale of the 2D metallicity, pointing to the distinct role of short-range versus long-range disorder in the 2D metallic transport in this correlated 2D hole system with interaction parameter r_s~ 20.


title: Spherical collapse in Galileon gravity: fifth force solutions, halo mass function and halo bias

abstract: We study spherical collapse in the Quartic and Quintic Covariant Galileon gravity models within the framework of the excursion set formalism. We derive the nonlinear spherically symmetric equations in the quasi-static and weak-field limits, focusing on model parameters that fit current CMB, SNIa and BAO data. We demonstrate that the equations of the Quintic model do not admit physical solutions of the fifth force in high density regions, which prevents the study of structure formation in this model. For the Quartic model, we show that the effective gravitational strength deviates from unity at late times ($z \lesssim 1$), becoming larger if the density is low, but smaller if the density is high. This shows that the Vainshtein mechanism at high densities is not enough to screen all of the modifications of gravity. This makes halos that collapse at $z \lesssim 1$ feel an overall weaker gravity, which suppresses halo formation. However, the matter density in the Quartic model is higher than in standard $\Lambda$CDM, which boosts structure formation and dominates over the effect of the weaker gravity. In the Quartic model there is a significant overabundance of high-mass halos relative to $\Lambda$CDM. Dark matter halos are also less biased than in $\Lambda$CDM, with the difference increasing appreciably with halo mass. However, our results suggest that the bias may not be small enough to fully reconcile the predicted matter power spectrum with LRG clustering data.


title: Sailfish: Alignment-free Isoform Quantification from RNA-seq Reads using Lightweight Algorithms

abstract: RNA-seq has rapidly become the de facto technique to measure gene expression. However, the time required for analysis has not kept up with the pace of data generation. Here we introduce Sailfish, a novel computational method for quantifying the abundance of previously annotated RNA isoforms from RNA-seq data. Sailfish entirely avoids mapping reads, which is a time-consuming step in all current methods. Sailfish provides quantification estimates much faster than existing approaches (typically 20-times faster) without loss of accuracy.


title: Theory of metallic double perovskites with spin orbit coupling and strong correlations; application to ferrimagnetic Ba2FeReO6

abstract: We consider a model of the double perovskite Ba2FeReO6, a room temperature ferrimagnet with correlated and spin-orbit coupled Re t2g electrons moving in the background of Fe moments stabilized by Hund's coupling. We show that for such 3d/5d double perovskites, strong correlations on the 5d-element (Re) are essential in driving a half-metallic ground state. Incorporating both strong spin-orbit coupling and the Hubbard repulsion on Re leads to a band structure consistent with ab initio calculations. Using our model, we find a large spin polarization at the Fermi level, and obtain a semi-quantitative understanding of the saturation magnetization of Ba2FeReO6, as well as X-ray magnetic circular dichroism data indicating a significant orbital magnetization. Based on the orbital populations obtained in our theory, we predict a specific doping dependence to the tetragonal distortion accompanying ferrimagnetic order. Finally, the combination of a net magnetization and spin-orbit interactions is shown to induce Weyl nodes in the band structure, and we predict a significant intrinsic anomalous Hall effect in hole-doped Ba2FeReO6. The uncovered interplay of strong correlations and spin-orbit coupling lends partial support to our previous work, which used a local moment description to capture the spin wave dispersion found in neutron scattering measurements. Our work is of broad interest for understanding metallic 4d-based and 5d-based double perovskites which are of fundamental interest and of possible relevance to spintronic applications.


title: Regional variations in the dense gas heating and cooling in M51 from Herschel far-infrared spectroscopy

abstract: We present Herschel PACS and SPIRE spectroscopy of the most important far-infrared cooling lines in M51, [CII](158 \mu m), [NII](122 & 205 \mu m), [OI](63 and 145 \mu m) and [OIII](88 \mu m). We compare the observed flux of these lines with the predicted flux from a photon dominated region model to determine characteristics of the cold gas such as density, temperature and the far-ultraviolet radiation field, G_0, resolving details on physical scales of roughly 600 pc. We find an average [CII]/F_TIR of 4 x 10^{-3}, in agreement with previous studies of other galaxies. A pixel-by-pixel analysis of four distinct regions of M51 shows a radially decreasing trend in both the far-ultraviolet (FUV) radiation field, G_0 and the hydrogen density, n, peaking in the nucleus of the galaxy, then falling off out to the arm and interarm regions. We see for the first time that the FUV flux and gas density are similar in the differing environments of the arm and interarm regions, suggesting that the inherent physical properties of the molecular clouds in both regions are essentially the same.


title: ALMA observations of the variable 12CO/13CO ratio around the asymptotic giant branch star R Sculptoris

abstract: [abridged] The 12CO/13CO ratio is often used as a measure of the 12C/13C ratio in the circumstellar environment, carrying important information about the stellar nucleosynthesis. External processes can change the 12CO and 13CO abundances, and spatially resolved studies of the 12CO/13CO ratio are needed to quantify the effect of these processes on the globally determined values. Additionally, such studies provide important information on the conditions in the circumstellar environment. The detached-shell source R Scl, displaying CO emission from recent mass loss, in a binary-induced spiral structure as well as in a clumpy shell produced during a thermal pulse, provides a unique laboratory for studying the differences in CO isotope abundances throughout its recent evolution. We observed both the 12CO(J=3-2) and the 13CO(J=3-2) line using ALMA. We find significant variations in the 12CO/13CO intensity ratios and consequently in the abundance ratios. The average CO isotope abundance ratio is at least a factor three lower in the shell (~19) than that in the present-day (<300 years) mass loss (>60). Additionally, variations in the ratio of more than an order of magnitude are found in the shell itself. We attribute these variations to the competition between selective dissociation and isotope fractionation in the shell, of which large parts cannot be warmer than ~35 K. However, we also find that the 12CO/13CO ratio in the present-day mass loss is significantly higher than the 12C/13C ratio determined in the stellar photosphere from molecular tracers (~19). The origin of this discrepancy is still unclear, but we speculate that it is due to an embedded source of UV-radiation that is primarily photo-dissociating 13CO. This radiation source could be the hitherto hidden companion. Alternatively, the UV-radiation could originate from an active chromosphere of R Scl itself....


title: Searching for Oscillations in the Primordial Power Spectrum: Perturbative Approach (Paper I)

abstract: In this first of two papers, we present a new method for searching for oscillatory features in the primordial power spectrum. A wide variety of models predict these features in one of two different flavors: logarithmically spaced oscillations and linearly spaced oscillations. The proposed method treats the oscillations as perturbations on top of the scale-invariant power spectrum, allowing us to vary all cosmological parameters. This perturbative approach reduces the computational requirements for the search as the transfer functions and their derivatives can be precomputed. We show that the most significant degeneracy in the analysis is between the distance to last scattering and the overall amplitude at low frequencies. For models with logarithmic oscillations, this degeneracy leads to an uncertainty in the phase. For linear spaced oscillations, it affects the frequency of the oscillations. In this first of two papers, we test our code on simulated Planck-like data, and show we are able to recover fiducial input oscillations with an amplitude of a few times order 10^{-2}. We apply the code to WMAP9-year data and confirm the existence of two intriguing resonant frequencies for log spaced oscillations. For linear spaced oscillations we find a single resonance peak. We use numerical simulations to assess the significance of these features and conclude that the data do not provide compelling evidence for the existence of oscillatory features in the primordial spectrum.


title: Searching for Oscillations in the Primordial Power Spectrum: Constraints from Planck (Paper II)

abstract: We apply our recently developed code to search for resonance features in the Planck CMB temperature data. We search both for log spaced oscillations or linear spaced oscillations and compare our findings with results of our WMAP9 analysis and the Planck team analysis. While there are hints of log spaced resonant features present in the WMAP9 data, the significance of these features weaken with more data. With more accurate small scale measurements, we also find that the best fit frequency has shifted and the amplitude has been reduced. We confirm the presence of a several low frequency peaks, earlier identified by the Planck team, but with a better improvement of fit (delta chi^2 ~ 12). We further investigate this improvement by allowing the lensing potential to vary as well, showing mild correlation between the amplitude of the oscillations and the lensing amplitude. We find that the improvement of the fit increases even more (delta chi^2 ~ 14) for the low frequencies that modify the spectrum in a way that mimics the lensing effect. Since these features were not present in the WMAP data, they are primarily due to better measurements of Planck at small angular scales. For linear spaced oscillations we find a maximum delta chi^2 ~ 13 scanning two orders of magnitude in frequency space, and the biggest improvements are at extremely high frequencies. We recover a best fit frequency very close to the one found in WMAP9, which confirms that the fit improvement is driven by low l. Further comparisons with WMAP9 show Planck contains many more features, both for linear and log space oscillations, but with a smaller improvement of fit. We discuss the improvement as a function of the number of modes and study the effect of the 217 GHz map, which appears to drive most of the improvement for log spaced oscillations. We conclude that none of the detected features are statistically significant.


title: Quantifying mass segregation and new core radii for 54 milky way globular clusters

abstract: We present core radii for 54 Milky Way globular clusters determined by fitting King-Michie models to cumulative projected star count distributions. We find that fitting star counts rather than surface brightness profiles produces results that differ significantly due to the presence of mass segregation. The sample in each cluster is further broken down into various mass groups, each of which is fit independently, allowing us to determine how the concentration of each cluster varies with mass. The majority of the clusters in our sample show general agreement with the standard picture that more massive stars will be more centrally concentrated. We find that core radius vs. stellar mass can be fit with a two parameter power-law. The slope of this power-law is a value that describes the amount of mass segregation present in the cluster, and is measured independently of our distance from the cluster. This value correlates strongly with the core relaxation time and physical size of each cluster. Supplementary figures are also included showing the best fits and likelihood contours of fit parameters for all 54 clusters.


title: Galaxy pairs in the Sloan Digital Sky Survey - VIII: The observational properties of post-merger galaxies

abstract: In order to investigate the effects of galaxy mergers throughout the interaction sequence, we present a study of 10,800 galaxies in close pairs and a smaller sample of 97 post-mergers identified in the Sloan Digital Sky Survey. We find that the average central star formation rate (SFR) enhancement (x 3.5) and the fraction of starbursts (20 per cent) peak in the post-merger sample. The post-mergers also show a stronger deficit in gas phase metallicity than the closest pairs, being more metal-poor than their control by -0.09 dex. Combined with the observed trends in SFR and the timescales predicted in merger simulations, we estimate that the post-mergers in our sample have undergone coalescence within the last few hundred Myr. In contrast with the incidence of star-forming galaxies, the frequency of active galactic nuclei (AGN) peaks in the post-mergers, outnumbering AGN in the control sample by a factor of 3.75. Moreover, amongst the galaxies that host an AGN, the black hole accretion rates in the closest pairs and post-mergers are higher by a factor of ~3 than AGN in the control sample. These results are consistent with a picture in which star formation is initiated early on in the encounter, with AGN activity peaking post-coalescence.


title: Multidimensional optomechanical cantilevers for high frequency atomic force microscopy

abstract: High-frequency atomic force microscopy has enabled extraordinary new science through large bandwidth, high speed measurements of atomic and molecular structures. However, traditional optical detection schemes restrict the dimensions, and therefore the frequency, of the cantilever - ultimately setting a limit to the time resolution of experiments. Here we demonstrate optomechanical detection of low-mass, high-frequency nanomechanical cantilevers (up to 20 MHz) that surpass these limits, anticipating their use for single-molecule force measurements. These cantilevers achieve 2 fm / sqrt(Hz) displacement noise floors, and force sensitivity down to 132 aN / sqrt(Hz). Furthermore, the ability to resolve both in-plane and out-of-plane motion of our cantilevers opens the door for ultrasensitive multidimensional force spectroscopy, and optomechanical interactions, such as tuning of the cantilever frequency in situ, provide new opportunities in high-speed, high-resolution experiments.


title: Multifrequency studies of the narrow-line Seyfert 1 galaxy SBS 0846+513

abstract: The narrow-line Seyfert 1 galaxy SBS 0846+513 was first detected by the Large Area Telescope (LAT) on-board Fermi in 2011 June-July when it underwent a period of flaring activity. Since then, as Fermi continues to accumulate data on this source, its flux has been monitored on a daily basis. Two further gamma-ray flaring episodes from SBS 0846+513 were observed in 2012 May and August, reaching a daily peak flux integrated above 100 MeV of (50+/-12)x10^-8 ph/cm^2/s, and (73+/-14)x10^-8 ph/cm^2/s on May 24 and August 7, respectively. Three outbursts were detected at 15 GHz by the Owens Valley Radio Observatory 40-m telescope in 2012 May, 2012 October, and 2013 January, suggesting a complex connection with the gamma-ray activity. The most likely scenario suggests that the 2012 May gamma-ray flare may not be directly related to the radio activity observed over the same period, while the two gamma-ray flaring episodes may be related to the radio activity observed at 15 GHz in 2012 October and 2013 January. The gamma-ray flare in 2012 May triggered Swift observations that confirmed that SBS 0846+513 was also exhibiting high activity in the optical, UV and X-ray bands, thus providing a firm identification between the gamma-ray source and the lower-energy counterpart. We compared the spectral energy distribution (SED) of the flaring state in 2012 May with that of a quiescent state. The two SEDs, modelled as an external Compton component of seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. No significant evidence of thermal emission from the accretion disc has been observed. Interestingly, in the 5 GHz radio luminosity vs. synchrotron peak frequency plot SBS 0846+513 seems to lie in the flat spectrum radio quasar part of the so-called `blazar sequence'.


title: The Second Main Theorem Vector for the modular regular representation of $C_2$

abstract: We study the ring of invariants for a finite dimensional representation $V$ of the group $C_2$ of order 2 in characteristic $2$. Let $\sigma$ denote a generator of $C_2$ and $\{x_1,y_1 \dots, x_m,y_m\}$ a basis of $V^*$. Then $\sigma(x_i) = x_i$, and $\sigma(y_i) = y_i + x_i$. To our knowledge, this ring (for any prime $p$) was first studied by David Richman in 1990. He gave a first main theorem for $(V_2, C_2)$, that is, he proved that the ring of invariants when $p=2$ is generated by $\{x_i, N_i = y_i^2 + x_iy_i, tr(A) | 2 \le |A| \le m\}$ where $A \subset \{0,1\}^m$, $y^A = y_1^{a_1} y_2^{a_2} \cdots y_m^{a_m}$ and $tr(A) = y^A + (y_1+x_1)^{a_1}(y_2+x_2)^{a_2} \cdots (y_m+x_m)^{a_m}.$ In this paper, we prove the second main theorem for $(V_2, C_2)$, that is, we show that all relations between these generators are generated by relations of type I: $\sum_{I \subset A } x^I tr(A-I)$ and of type II: $tr(A) tr(B) = \sum_{L < I} x^{I-L} N^L tr(I-L+J+K) + N^I \sum_{L < J} x^{J-L} tr(L+K)$ for all $m$. We also derive relations of type III which are simpler and can be used in place of the relations of type II.


title: Effective electrostatic interactions in mixtures of charged colloids

abstract: We present a theory of effective electrostatic interactions in polydisperse suspensions of charged macroions, generalizing to mixtures a theory previously developed for monodisperse suspensions. Combining linear response theory with a random phase approximation for microion correlations, we coarse-grain the microion degrees of freedom to derive general expressions for effective macroion-macroion pair potentials and a one-body volume energy. For model mixtures of charged hard-sphere colloids, we give explicit analytical expressions. The resulting effective pair potentials have the same general form as predicted by linearized Poisson-Boltzmann theory, but consistently incorporate dependence on macroion density and excluded volume via the Debye screening constant. The volume energy, which depends on the average macroion density, contributes to the free energy and so can influence thermodynamic properties of deionized suspensions. To validate the theory, we compute radial distribution functions of binary mixtures of oppositely charged colloidal macroions from molecular dynamics simulations of the coarse-grained model (with implicit microions), taking effective pair potentials as input. Our results agree closely with corresponding results from more computationally intensive Monte Carlo simulations of the primitive model (with explicit microions). Simulations of a mixture with large size and charge asymmetries indicate that charged nanoparticles can enhance electrostatic screening of charged colloids. The theory presented here lays a foundation for future large-scale modeling of complex mixtures of charged colloids, nanoparticles, and polyelectrolytes.


title: Tracking Down the Top Quark Forward-Backward Asymmetry with Monotops

abstract: We revisit the possibility that the top quark forward-backward asymmetry arises from the on-shell production and decay of scalar top partners to top-antitop pairs with missing transverse energy. Although the asymmetry is produced by t-channel exchange of a light mediator, the model remains unconstrained by low energy atomic parity violation tests. An interesting connection to the active neutrino sector through a Type-I seesaw operator helps to evade stringent monojet constraints and opens up a richer collider phenomenology. After performing a global fit to top data from both the Tevatron and the LHC, we obtain a viable region of parameter space consistent with all phenomenological and collider constraints. We also discuss the discovery potential of a predicted monotop signal and related lepton charge asymmetry at the LHC.


title: Heavy dust obscuration of z=7 galaxies in a cosmological hydrodynamic simulation

abstract: Hubble Space Telescope observations with Wide Field Camera 3/IR reveal that galaxies at z~7 have very blue ultraviolet (UV) colors, consistent with these systems being dominated by young stellar populations with moderate or little attenuation by dust. We investigate UV and optical properties of the high-z galaxies in the standard cold dark matter model using a high-resolution adaptive mesh refinement cosmological hydrodynamic simulation. For this purpose, we perform panchromatic three-dimensional dust radiative transfer calculations on 198 galaxies of stellar mass 5x10^8-3x10^{10} Msun with three parameters, the dust-to-metal ratio, the extinction curve, and the fraction of directly escaped light from stars (\fesc). Our stellar mass function is found to be in broad agreement with Gonzalez et al., independent of these parameters. We find that our heavily dust-attenuated galaxies (A_V~1.8) can also reasonably match modest UV-optical colors, blue UV slopes, as well as UV luminosity functions, provided that a significant fraction (~10%) of light directly escapes from them. The observed UV slope and scatter are better explained with an Small Magellanic Cloud-type extinction curve, whereas Milky Way-type dust predicts too blue UV colors due to the 2175A bump. We expect that upcoming observations by ALMA will be able to test this heavily obscured model.


title: Randomization and the Gross-Pitaevskii hierarchy

abstract: We study the Gross-Pitaevskii hierarchy on the spatial domain $\mathbb{T}^3$. By using an appropriate randomization of the Fourier coefficients in the collision operator, we prove an averaged form of the main estimate which is used in order to contract the Duhamel terms that occur in the study of the hierarchy. In the averaged estimate, we do not need to integrate in the time variable. An averaged spacetime estimate for this range of regularity exponents then follows as a direct corollary. The range of regularity exponents that we obtain is $\alpha>\frac{3}{4}$. It was shown in our previous joint work with Gressman that the range $\alpha>1$ is sharp in the corresponding deterministic spacetime estimate. This is in contrast to the non-periodic setting, which was studied by Klainerman and Machedon, in which the spacetime estimate is known to hold whenever $\alpha \geq 1$. The goal of our paper is to extend the range of $\alpha$ in this class of estimates in a \emph{probabilistic sense}. We use the new estimate and the ideas from its proof in order to study randomized forms of the Gross-Pitaevskii hierarchy. More precisely, we consider hierarchies similar to the Gross-Pitaevskii hierarchy, but in which the collision operator has been randomized. For these hierarchies, we show convergence to zero in low regularity Sobolev spaces of Duhamel expansions of fixed deterministic density matrices. We believe that the study of the randomized collision operators could be the first step in the understanding of a nonlinear form of randomization.


title: Quantum kinetics of spinning neutral particles: General theory and Spin wave dispersion

abstract: Plasma physics give an example of physical system of particles with the long range interaction. At small velocity of particles we can consider the plasma approximately as a system of particles with the Coulomb interaction. The Coulomb interaction is isotropic. Systems of spinning neutral particles have long-range anisotropic interparticle interaction. So, they can reveal more reach properties than plasma. Furthermore for studying of systems of spinning particles we can develop kinetic and hydrodynamic methods analogous to used for the plasma. We derive kinetic equations by a new method, which is the generalization of the many-particle quantum hydrodynamics. Obtained set of kinetic equations is truncated, so we have closed set of two equations. One of them is the kinetic equation for quantum distribution function. The second equation is the equation for the spin-distribution. Which describes the spin kinetic evolution and gives contribution in time evolution of the distribution function. Our method allows to obtain equations as for three dimensional system of particles and for low dimensional systems. So, we consider spin waves in three- and two dimensional systems of neutral spinning particles.

categories: physics.plasm-ph cond-mat.quant-gas cond-mat.stat-mech quant-ph


title: Gravitational Dynamics From Entanglement "Thermodynamics"

abstract: In a general conformal field theory, perturbations to the vacuum state obey the relation \delta S = \delta E, where \delta S is the change in entanglement entropy of an arbitrary ball-shaped region, and \delta E is the change in "hyperbolic" energy of this region. In this note, we show that for holographic conformal field theories, this relation, together with the holographic connection between entanglement entropies and areas of extremal surfaces and the standard connection between the field theory stress tensor and the boundary behavior of the metric, implies that geometry dual to the perturbed state satisfies Einstein's equations expanded to linear order about pure AdS.


title: Persistent spin excitations in doped antiferromagnets revealed by resonant inelastic light scattering

abstract: How coherent quasiparticles emerge by doping quantum antiferromagnets is a key question in correlated electron systems, whose resolution is needed to elucidate the phase diagram of copper oxides. Recent resonant inelastic X-ray scattering (RIXS) experiments in hole-doped cuprates have purported to measure high-energy collective spin excitations that persist well into the overdoped regime and bear a striking resemblance to those found in the parent compound, challenging the perception that spin excitations should weaken with doping and have a diminishing effect on superconductivity. Here we show that RIXS at the Cu L3-edge indeed provides access to the spin dynamical structure factor once one considers the full influence of light polarization. Further we demonstrate that high-energy spin excitations do not correlate with the doping dependence of Tc, while low-energy excitations depend sensitively on doping and show ferromagnetic correlations. This suggests that high-energy spin excitations are marginal to pairing in cuprate superconductors.


title: The dynamics of internal gravity waves in the ocean: theory and applications

abstract: In this paper we consider fundamental processes of the disturbance and propagation of internal gravity waves in the ocean modeled as a vertically stratified, horizontally non-uniform, and non-stationary medium. We develop asymptotic methods for describing the wave dynamics by generalizing the spatiotemporal ray-tracing method (a geometrical optics method). We present analytical and numerical algorithms for calculating the internal gravity wave fields using actual ocean parameters such as physical characteristics of the sea water, topography of its floor, etc. We demonstrate that our mathematical models can realistically describe the internal gravity wave dynamics in the ocean. Our numerical and analytical results show that the internal gravity waves have a significant impact on underwater objects in the ocean.

categories: physics.flu-dyn math-ph math.MP physics.ao-ph


title: Consequences of Mechanical and Radiative Feedback from Black Holes in Disc Galaxy Mergers

abstract: We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $\epsilon_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - \sigma$) as well as X-ray luminosity ($L_X - \sigma$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim \sigma \lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model.


title: Probing the Gamma-Ray Burst Rate with Trigger Simulations of the Swift Burst Alert Telescope

abstract: The long gamma-ray burst (GRB) rate is essential for revealing the connection between GRBs, supernovae and stellar evolution. Additionally, the GRB rate at high redshift provides a strong probe of star formation history in the early universe. While hundreds of GRBs are observed by Swift, it remains difficult to determine the intrinsic GRB rate due to the complex trigger algorithm of Swift. Current studies usually approximate the Swift trigger algorithm by a single detection threshold. However, unlike the previously flown GRB instruments, Swift has over 500 trigger criteria based on photon count rate and additional image threshold for localization. To investigate possible systematic biases and explore the intrinsic GRB properties, we developed a program that is capable of simulating all the rate trigger criteria and mimicking the image trigger threshold. We use this program to search for the intrinsic GRB rate. Our simulations show that adopting the complex trigger algorithm of Swift increases the detection rate of dim bursts. As a result, we find that either the GRB rate is much higher than previously expected at large redshift, or the luminosity evolution is non-negligible. We will discuss the best results of the GRB rate in our search, and their impact on the star-formation history.


title: Mechanical cell-matrix feedback explains pairwise and collective endothelial cell behavior in vitro

abstract: In vitro cultures of endothelial cells are a widely used model system of the collective behavior of endothelial cells during vasculogenesis and angiogenesis. When seeded in an extracellular matrix, endothelial cells can form blood vessel-like structures, including vascular networks and sprouts. Endothelial morphogenesis depends on a large number of chemical and mechanical factors, including the compliancy of the extracellular matrix, the available growth factors, the adhesion of cells to the extracellular matrix, cell-cell signaling, etc. Although various computational models have been proposed to explain the role of each of these biochemical and biomechanical effects, the understanding of the mechanisms underlying in vitro angiogenesis is still incomplete. Most explanations focus on predicting the whole vascular network or sprout from the underlying cell behavior, and do not check if the same model also correctly captures the intermediate scale: the pairwise cell-cell interactions or single cell responses to ECM mechanics. Here we show, using a hybrid cellular Potts and finite element computational model, that a single set of biologically plausible rules describing (a) the contractile forces that endothelial cells exert on the ECM, (b) the resulting strains in the extracellular matrix, and (c) the cellular response to the strains, suffices for reproducing the behavior of individual endothelial cells and the interactions of endothelial cell pairs in compliant matrices. With the same set of rules, the model also reproduces network formation from scattered cells, and sprouting from endothelial spheroids. Combining the present mechanical model with aspects of previously proposed mechanical and chemical models may lead to a more complete understanding of in vitro angiogenesis.


title: On the internal gravity waves in the stratified ocean with shear flows

abstract: In this paper, we consider a fundamental problem of describing the dynamics of internal gravity waves in the stratified ocean with shear flows. We develop an asymptotic representation of the wave fields in terms of the Green's functions. We explore the far field of the internal gravity waves generated by disturbing sources, and propose asymptotic algorithms for calculating its dynamics.